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SECTION 1 
INTRODUCTION 


Vertical takeoff and landing (VTOL) aircraft have considerable potential 
for use in a viable short-haul air transportation system. The VTOL aircraft used in this 
context would provide convenient, safe and reliable access to long-haul air transporta- 
tion by providing air service from major and smaller cities to regional airports. They 
would also contribute to the achievement of a more balanced total transportation system 
by providing direct links between smaller cities and major cities and between nearby 
major cities. 

In order for such a VTOL short-haul system to be economically feasible, the 
aircraft must provide schedule relability in all-weather conditions, must provide accept- 
able levels of ride quality, and must be operated directly into the city centers to provide 
the requisite convenience to passengers. Before a viable VTOL system can become a 
reality, technology developments are needed in a number of areas. During the past 
several years many advanced VTOL aircraft design prc^rams have been carried out by 
NASA, DOD and the aircraft industry to develop economical vehicles with improved 
ride qualities and controllability which would be suitable for a commercial VTOL trans- 
portation system. However, to effectively utilize these vehicles and to exploit their 
unique characteristics for minimizing noise and both air and ground space requirements, 
corresponding advances must be made in handling qualities, operating procedures, and 
all-weather avionics. 

The NASA Langley Research Center (LaRC) has undertaken a research progrom 
to develop the navigation, guidance, control, display and fl ight management tech- 
nology base needed by Government and industry in establishing systems design con- 
cepts and operating procedures for VTOL short-haul transportation systems in the 1980s 
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time period and beyond. The VALT (VTOL Automatic _Londing Technology) Program 
encompasses the investigation of operating systems and piloting techniques associated 
with VTOL operations under all-weather conditions from downtown vertiports; the 
investigation of terminal air traffic and airspace requirements; and the development 
of avionics including navigation, guidance, controls, and displays for automated 
takeoff, cruise, and landing operations. 

In support of the VALT Program, Aerospace Systems, Inc. (ASI) has conducted 
a number of research studies for LaRC . which provide a technology base for the present 
study. In the initial effort (Reference 1), ASI analyzed the navigation and guidance 
requirements for commercial VTOL operations in the takeoff, cruise, terminal area, 
and landing phases of flight in weather conditions up to and including Category III. 

A digital computer simulation (Program VALT) was developed to provide a means for 
evaluating the performance of candidate VTOL guidance and control systems. This 
program was used to conduct a sensitivity study of several VTOL guidance and control 
system concepts (Reference 2). 

One conclusion in Reference 1 was that curved decelerating approaches will 
be required for safe, efficient, and independent VTOL operations. To facilitate these 

' ' • ' ... ■ t 

maneuvers, a spiral descent technique was formulated as a possible standard VTOL 
approach procedure. The spiral descent investigated by ASI in Reference 3 uses minimal 
airspace, accommodates arrivals from any direction, and can service multipad landings. 

The spiral approach also provides the benefits of a vertical descent, but avoids the 
vortex ring state, maintains a stable airspeed, and uses less fuel. 

The guidance of a VTOL along this type of spiral descent trajectory is compli- 
cated by many control problems which are unique to this class of aircraft. To reduce the 
workload for the guidance and control tasks to a tolerable level for multiple daily land- 
ings, the aircraft controls will be partially or completely automated. As the level of 
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automation increases, the pilot's role shifts from primarily that of d controller towards 
that of a system monitor and manager. The purpose of the present study was to examine 
which tasks should be allocated to the pilot of an automated VTOL aircraft utilized as 
part of a short-haul air transportation system and to determine what displayed information 
will be required in performing these tasks. 


While the study was intended to provide insight into problems associated with 
VTOL pilot tasks in an automated VTOL aircraft in general, several guidelines were 
specified by NASA to provide a frame of reference and to ensure that the results could 
be readily used and evaluated in VALT and other existing Langley Research Center 
programs. These guidelines included the following; 

• Flight Profile - The emphasis of the study was to be on the approach 
and landing phase of flight; however, sufficient general consideration 
to the takeoff and enroute phases of flight was to be included to ensure 
that the study results would be compatible with the overall task of 
operating the vehicle as a commercial transport. 

• Vehicle Dynamics - The study was to utilize the CH-46C and CH-47 
helicopters used in the LaRC flight research programs. 

• Crew - Crew tasks were to be configured to permit operation by one 
pilot. Routine calls, communication channel selection, or other 
tasks which might be handled by a second crew member in an opera- 
tional context were not included in the scope of work. 

• Pilot Involvement - The levels of automation considered were to be 
varied over a range extending from a fully automatic system with the 
pilot in a passive mode with respect to control activity to a system 
with full manual control. 

• Technology Date - In defining a level of system automation, allocating 
tasks to automatic systems, and in conceiving displays for the control/ 
display concept, decisions were to be based on the relevant technology 
projected as being available in the mid 1980s. 

• Pilot/Hardware Experiments - Hardware tests, flight tests, and pilot/ 
hardware interaction experiments were to be specifically excluded 
from the scope of the work . 
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Some preliminary results obtained during the present study were' included in 
a paper presented by the principal investigators in Reference 4. This report includes d 
complete description of the work performed and the results obtained during the investi- 
gation. Section 2 discusses the VTOL display/control methodology used, including a 
description of the model developed for simultaneous monitoring and control. Section 3 
includes a discussion of the logic used for the choice of control systems, a description 
of the different levels of automation investigated, and the models that were formulated. 
The different approaches for pilot monitoring performance prediction and choices of 
monitoring workload performance metrics are examined in Section 4. Section 5 pre- 
sents the model validation analysis results which utilize data obtained by NASA during 
CH-46C helicopter flight tests. Results obtained for eight CH-47 levels of control 
automation are presented in Section 6. . In Section 7 a CH-47 display concept and format 
is presented . Conclusions and recommendations are given in Sections 8 and 9 respectively, 
followed by a List of References and a Bibliography. The Appendices contain detailed 
technical information concerning the equations for the .optimal control model for the 
human operator, attenfional allocation using the optimal control model, and CH-46C 
and CH-47 aircraft characteristics used in the study. 
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SECTION 2 


VTOL DISPLAY/CONTROL DESIGN METHODOLOGY 

The pilot's functions throughout the entire flight of a VTOL aircraft can be 
broadly classified in the two areas of (1) control and (2) monitoring; both of these func- 
tions are necessary for safe and efficient operations. This study has addressed itself to 
the development of a pilot/vehicle model for simultaneous monitoring and control in 
order to explore the dispiqy/control tradeoffs inherent in the selection of a display and 
control system for an automated VTOL aircraft. The advantage of using the model in 
the preliminary stages is that a wide variety of candidate systems can be explored with 
a minimal amount of effort and cost. The detailed man-in-the-loop simulations can and 
should be reserved to resolve the minor details between competing display/control sys- 
tems and to confirm the predictions of the model . 

Regardless of the form of the model used to carry out the display/control 
design, it is imperative that this model have realistic and quantitative metrics for the 
following: 

• System control performance 

• Workload for control 

• Monitoring performance 

• Workload for monitoring 

I , 

Without measures for these four quantities, one would be unable to explore the many 
and varied tradeoffs between control, monitoring, augmentation systems and displays. 

2. 1 A MODEL FOR SIMULTANEOUS MONITORING AND CONTROL 

This section briefly outlines a model for simultaneous monitoring and control 
which is used in the design methodology outlined in Subsection 2.2. This model is based 
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on fhe optimal control model of the human operator (References Sand 6)/ which is 
described in more detail in Appendix A. 

The performance metric for the control task is given by 


j =T, /./ 

C X. I 

i=l ' 


max 


(2-1) 


where = performance pietric, for control 

o = - standard deviation for X. 

X. ■ I 

I 

ith 

X. = i component of state vector x 


Equation (2-1) is a relative weighting of the variances of the individual components of 
the state vector normalized by their maximum allowable excursions. When evaluating 
competing systems, one may use a subset of the state vector components; this is , 
equivalent to assuming that some components have infinite allowable deviations. 


The workload metric for control is based on the task -interference model of 
Reference 7. This metric states that for a single display element, the covariance of 
the observation noise at the input to the pilot model (see Appendix A) is given by 




( 2 - 2 ) 


th 

where the noise-to-signal ratio pj is approximately 0.01, and y. is the i*^" displayed 
variable. When the pilot devotes only a part of his attention to any one display 
element, then the observation noise can be modelled as 


V = 


p. CT 

' Xt 


i f ■ 

c. 


(2-3) 
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th 

where is the fraction of attention allocated to the i display. This observation 
i 

noise has the appropriate limiting values: if f^. ~ 1/ then-Equation (2-3) reduces to 

the full-attention results. Equation (2-2); when the controller is not using the i*^^ 

display, then.fg^ = 0, and the observation noise is infinite, corresponding to no 

observations at all. The fractions of attention fg must all be positive and sum to the 

i 

total fraction of attention being used for control . Thus the control workload metric is 
given by 


f 

c 



fc. ^ “• 


(2-4) 


The pilot allocates his attention between the displays, spending the larger 
fractions of attentions on displays which are most useful for control. This behavior is 
formulated in the model by the assumption that the pilot minimizes a quadratic per- 
formance index with respect to the f^. subject to the constraints of Equation (2-4), i.e.. 


min J = min E 


(U 

I 


<'c> 


( 

T . . ) 

lim 1 

j jx' Q^x+ u' Q^uldtJ 

)T-»oo T 

V 

0 ^ 


(2-5) 


where and are positive semidefinite weighting matrices. Since the performance 
index in Equation (2-5) can be evaluated for any specific values of f^^ in Equation 
(2-3), it is conceptually a simple problem to further minimize this cost function subject 
to the constraints of Equation (2-4). Appendix B contains details of the associated 
equations and the techniques for performing this minimization. 


The basic form of the model for simultaneous monitoring and control can now 

be described. For a given level of control task workload, f , the pilot will allocate 

c 

his attention as described by Equation (2-5) to minimize the performance index. This 
is first in the hierarchy of control and monitoring, i.e., the pilot will first attend to 
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the control task, and with any available attention remaining, will then attend to the 
monitoring task. Thus the control task workload metric of Equation (2-4) and the 
control task performance metric of Equation (2-1). can be used to draw the control 
performance/workload tradeoff curve. At any given level of system performance (or 
workload), the pilot will have some residual capacity or fraction of attention after 
attending to the control task. We assume a total fraction of attention of 0.8 should 
be allocated to both the control and monitoring tasks: e.g., if the control task 
requires a fraction of attention of 0.3, a fraction of attention of 0.5 |s available for 
monitoring; similarly, if the control task requires f = 0.7, then only 0.1 fraction of 
attention is available for monitoring. Thus the control task workload metric and the 
monitoring task workload metric obey the constraint given in Equation (2-6). 


f + f 
c m 


= f. 


TOt 


0.8 


(2-6) 


2.2 THE OPTIMAL CONTROL MODEL OF THE PILOT AT THREE LEVELS 


This section includes a description of the development and use of the 
optimal control model of the pilot which was refined during the course of this study. 

The model can be exercised at three .levels of detail. These levels have been termed 
the "information level," the "display element level, "and the "display format level." 
Each of these levels and the use of the model in these levels will now be briefly 
described. 

2.2.1 INFORMATION LEVEL ■ 

At this level of operation, each of the elements of the state vector is 
assumed to be observable by the pilot. Each of these "observations" is assigned its own 
fraction of attention for control, and the optimization of the performance index. 
Equation (2-5), is carried out for various levels of control task workload f^. This 
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procedure leads to a very simple way of determining the information requirements for 
the task, since those state variables of greater importance will receive the- greater 
amount of attention , 

2.2.2 DISPLAY ELEMENT LEVEL 

The display element level is perhaps the most common level of operation for 
the optimal control model; each of the displayed quantities is assumed to be indicated 
by a display "element." The inputs to the model in this configuration are the position 
of the display element and its rate of change. Each of the display elements is assigned 
a given level of attention, so that the observation noise for both the display element 
position and display element rate are assigned the same fraction of attention for 
control, fg . Attention allocation is found by minimizing Equation (2-5). 

2.2.3 DISPLAY FORMAT LEVEL 

This level of use of the pilot model is very similar to the' display element 
level described above. However, in this case specific characteristics of the display 
format such as indifference thresholds, lack of a zero reference, saturation limits, 
smoothing filters for display augmentation, etc.,' are incorporated into the model. 

The details of the use of the model at these three levels are explained in 
the next section On display/control system design methodology. 

2.3 VTOL DISPLAY/CONTROL DESIGN PROCEDURE - 

The VTOL display/control design procedure developed and utilized in the 
study is a ten-step process. These ten steps can be broken down into four main 
categories; information requirements, control/monitor performance, pilot/automatic 
task allocation, and display format design. The ten steps, outlined In Table 2-1, are 
discussed in detail below. 
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. .. Table 2-J . VTOL Displqy/Control Design Procedure. .. 


CATEGORY- 

— STEP 

Information 
Requirements ' 

1 . Determine max x. and max J from’ 

1 c 

mission requirements. 

2. Select candidate control systems. 

3. Calculate J vs. f at the information 

c c 

level for each control system (Include 
steering commands). 

4. Choose display elements. 

Control/Monitor 

Performance 

5: Calculated vs. f' at the display 

c c 

element level for each display/control 
system . 

6. Determine f = frr\T “ f and J from 

m 1 1 c m 

monitoring model . 

Pilot/Automatic 
Task Allocation 

. 7. Select dlsplby/control system 

configuration . 

. - . Display Format . , 

Design 

8. Select display format candidates. 

9. Determine Jg, vs. f^.. for each 

display format. 

10. Select display format. 


Step 1. Determine max x; and max Jc from mission requirements . 

The maximum deviotions of. state vector components are selected from the 
mission requirements, pilpt acceptance criteria, and passenger acceptance criteria. 
These maximum values of x. are used to determine the coefficients in the quadratic 
cost functional. The maximum of the overall figure of merit of the system control 
performance, max J^, is a function of the upper limit of the maximum acceptable 
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error variances which are def-ermined by the mission requirements; i.e., for a given 

trajectory and attitude tolerance, there is a value of J which must not be exceeded 

for any system to be considered viable. For example^ if the lx. | values are 
' ' I ^ max ‘ 

exceeded with probability of 0.05, then |x. | may be interpreted as twice the 

maximum standard deviation under the Gaussian assumption. If each of the variances 
is at its maximum (a worst case design), then J is n/4, where n is the number of 

•C f maX’ . 

terms in Equation (2-1). 

Step 2; Select candidate control systems. 

The candidate control systems which utilize different levels of augmentation 
are selected and designed. The output of this step of the procedure is a set of 
control led-element dynamics which interact with the pilot. These levels of automa- 
tion may cover the complete spectrum from the raw unaugmented vehicle through 
complete position feedback, i.e., a fully automatic system. The details of the 
control system design procedures used for this study are described in Section 3. 

Step 3; Calculate Jg vs. f^. at the information level for each control system. 

This is the key step in determining information requirements for the candidate 
control systems. Each element of the state vector of the system is considered as a 
measurement, with its own fraction of attention for control, i.e., the display element 
and its rate are treated independently. For a given level of control attention f , the 
performance index is minimized subject to the constraint of Equation (2-4), 



f , f ^ 0 . 
c c. 


Included in the displayed information should be director/steering commands designed 
for each particular level of automation. (See Section 3 for the description of the 
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flight director design used in this study.) The specified level of control fraction of 
attention, f , does not carry the same meaning as it does when operating at the display 
element level of the model. However, when i is varied from 0 to the number of 
states, f. provides the relative importance of each of the states. These in turn are 
used as input to the next step of the design procedure. 

Step 4; Choose display elements . 

At this stage of the design procedure, the information requirements derived 
in the previous step are used to choose display elements. The most difficult part of 
this step is determining whether or not a separate display of rate is required for use by 
the pilot. For the pitch attitude display element, pitch and pitch rate are used as 
inputs to the optimal control model of the pilot; the same is true for altitude and 
altitude rate. However, it is empirically known that separate displays of pitch attitude 
and rate are not required, while a separate display of altitude rate is required. These 
differences are due to the accuracy with which the rate may be inferred from a position 
display; for pitch attitude, pitch rate may be discerned to sufficient accuracy from the 
indicator, whereas this is not triie for the altitude indicator. This is due to the inherent 
del, ay , in the barometric .altitude indication and the dynamic range of the displayed 
signal. New types of altimeter systems and displays mediate this somewhat, but ex-, 
pi icit vertical speed information |s still used in the majority of conventional instrument 
panels. , ; , , , , . 

The decision 'of whether or not to include a separate 'rate display can be 
inferred from the information level results of the previous step by examining the 
required accuracy of display. 


2-8 


AEROSPACE SYSTEMSp INC. • ONE vine brook park • Burlington. Massachusetts oiao3 • (ei7) 37b-7si7 




Step 5; Calculate ys f,. at the display element level for each disploy/ 
control system . 

This step is similar to Step 3, where the control performance metric, J , is 

c 

computed for different values of workload metric f^. This leads to curves of control 
performance vs. workload having the general shape shown in Figure 2-1. This curve is 
plotted for two systems, but in general there will be as many systems as there are levels 
of automation times the number of flight director options, i.e., if there were five 
levels of automation and two flight director options (with and without flight director), 
then there would be curves for ten competing display/control systems shown in Figure 
2-1 . It should be noted that Figure 2-?l is the plot of the minimum value of for a 
given f , i.e., the attention has been allocated to minimize the performance index. 
Thus additionaroutpyts of the model at this stage are the fractions of attention being 
spent on each of the display elements. To prevent overly optimistic predictions of per- 
formance at this stage> indifference thresholds should be incorporated in the model cal- 
culations (see Appendix A), since this is a well-known facet of pilot behavior. (In the 
present study, X. ^^^4 was used as a. representative value.) 
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Step 6: Determine n ^TOT “ and from. monitoring model . 

This step is also shown in Figure 2-1, where the value of total attention for 
the control and monitoring task,(fyQj) is specified. Since: the curves will be monotoni- 
colly decreasing (the more attention for control, the better the job of control one can 
do), a line has been drawn at the maximum value of control; performance, J . . 
consistent with the mission requirements. The intersect ion. of this line of constant con- 
trol performance and the J vs. f curve for each; of the candidate systems determines 

the minimum amount, of control attention required, f The difference between 

1 c, req 

this amount of attention, .and the total available for. the entire task is the residual 

fraction of attention available: for monitoring f 

< t f fTi/ avail 

. Thus of the four important metrics required for the display/control design 

process, three have been determined at this stage:' f , J , and f . The fraction of at- 
* c c m 

tention available for monitoring, is then applied to the monitoring model (see Section 4 
for a discussion of the monitoring models). - Since thO' monitoring performance will be 
monotonic with monitoring attention, one simple monitoring strategy is equal alloca- 
tion of attention to the monitoring instruments. 

Results of applying the monitoring model may be plotted as shown in Figure 
2-2, This shows the contours of constant monitoring performance in the control com- 
plexity-display complexity plane often used in discussions of the VTOL display/control 
design problem. Although curves of pilot workload are normally plotted in this plane, 
monitoring performance is an equivalent measure for our discussion. It should be noted 
that campeting systems, which appear as paints in this plane, are being compared an 
the basis of the same total workload and the same level of system performance. 
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Control 

Complexity 


Constant: 



Display 

Complexity 


Figure 2-2. Monitoring Performance in the Control -Display Plane. 


Perhaps a more meaningful presentation of the monitoring performance data 


is shown in Figure 2-3. This form of the data is easier to interpret when many tradeoffs 
are necessary during the system selection as described in the next step. 



Display System FDq FD^ FD 2 FDq FD-j FD 2 FDq FD^ FD 2 

Control System A B C 

Figure 2-3. Monitoring Performance for Candidate Systems Display/Control. 
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Step 7; Select display/control system configuration . 

Criteria for allocating the tasks between the pilot and the automatic system 
should be based on the following considerations ranked in order of importance. 

• Required workload level less than or equal to approximately 0.8 ' 

• Sensitivity of workload level to changes in display control configuration 
system failure modes 

• Sensitivity of performance to changes in control and monitoring workload 

• Cost 

An approximate level of 80% workload, on which to base the next phase of 
design provides an adequate although not excessive, margin for the more detailed design 
and analysis of the display format. The sensitivity of workload level to changes In dis- 
play and control configurations Is important because it minimizes the impact of model- 
ing errors ori' the ‘final decision and minor display or control mode changes will not 
place excessive workload bn the' pilot. A brief examination of systems failure modes 
is important during this part of the analysis because it provides an indication of 
whether or not the system can be controlled when either the display or the control 
system fails. However, Figures 2-2 and 2-3 are calculated for accuracy requirements 
that may be too stringent during the failed condition. 

The first of ibe allocation criteria is satisfied by ibe conditions under which 
the calculations were made, i.e.. Figures 2-2 and 2-3 hove been obtained for a total 
workload level of 0.8. The sensitivity to changes in display/control configurations 
can be determined by examining Figure .2-3. For example, should System A be chosen 
with the flight director 2 display option (FD 2 ), then under flight director failure a 
change of monitoring performance (while maintaining the control performance) is given 
by the difference in the values of for FD 2 and FDq conditions. Control system 
failures (while maintaining the same level of control performance) can be assessed by 
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examining the change from one configuration to another. For example, if the nominal 
system operates with the display of FD 2 and control system C, then a failure in the 
control system will revert the control task back to that of control system A with 
no flight director. (It is assumed that the flight director mode must be turned off if 
there is a change in control systems because of the difference in flight director gains 
for different control systems. That is, a flight director designed for control system C 
may not be adequate, and may even lead to deleterious performance, when used with 
control system A.) 

Step 8; Select display format candidates . 

The display formats, which may be one of three types (separated displays, 
perspective displays, or combined displays) should be selected according to the follow- 
ing guidelines discussed in Section 7: 

• Operator centered and oriented 

• Geometric "real world" compatibility 

• Naturalness (for high stress situations) 

• "Status at a glance" for situation displays 

• Predictive capability for situation displays 

'• Compactness 

• Lack of clutter 

In addition to these interpretability considerations, operational guidelines should be 
used which include failure mode consideration (the ability to make missed approaches 
with a minimum of control and display augmentation), display options (such as change 
in scale), and flexibility and versatility in trajectory selection. Details of this 
procedure are described in Section 7. 
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Step 9; Determine ,Jc> vs. fg for each display formaf-. 


The basic approach In evaluating the display format candidates is to perform 
the computations of Step 3 with a more detailed model and description to account for 
specific display formats and mechanizations. In general, one should expect the perfor- 
mance to degrade somewhat due to Hie practical aspects of implementing the information 
display. For a given display candidate, one should include the effect of scaling, 
thresholds, and a zero reference (or lack of one) by d detailed nbise-vafiance model. 
Note that if each of the format candidates is chosen to be consistent with the indiffer- 
ence thresholds used In the calculations of Step 5, there should be only minor changes 
in performance due to this effect. 

The specification of the observation noise variance model will take the form 


V (t) = -II 

y] f 

C. L 
I 




(2-7) 


.th . 


where is the observation noise covariance for the i input to the Kalman filter por- 
tion of the pilot model. is the nominal sole-task noise variance,. and f is the 

^ C| 

control-task fraction of attention allocated to this display variable. The quantity in 
the brackets is a detailed noise model that is a function of the particular display format 
being evaluated. Typically there are two effects as shown in Equation (2-7). 


2 

• K represents the describing functjon which predicts and accounts for 
noise varidndes due to thresholds/ saturations, and other nonlinearities 
of the display; 

2 

• is a variance representing the effect of (or lack of) a zero reference 

for the display. In other words if there is no zero reference for the dis- 
play and the pilot' has to judge the value of the displayed variable 
relative to the position in its total range of travel, then the observation 
noise of this variable will be greater thdn if he is provided a zero, or 
commanded reference. 
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A further change in the model should be the incorporation of realistic flight 
director computers; since these signals must be generated from noisy information, which 
often involves high-pass or differentiation of such data, requiring that the navigation 
data be filtered before presentation to the pilot to ovoid his rejection on the basis of 
noisy signals. Thus, one mujt specify the dynamics of the actual flight director to be 
used in the computations at this stage. 

Using the more detailed model of the display format, one may again deter- 
mine system control performance vs control task workload f^ for the different dis- 
play formats as shown in Figure 2-4. These curves are generated for several different 
formats of the display/control system selected in Step 7 and form the basis for making 
the final display format selection in the next step. In a manner similar to that of cal- 
culating the control performance J , one should recalculate the monitoring performance 

c 

metric J for every value of f , which is the residual of the fraction of attention after 
m ' m 

attending to the control task. 

Step TO: Select display format . 

The specific selection of the individual candidate displays should be based 
on the following criteria: 

• Workload level 

• System performance 

• Monitoring performance 

• Sensitivity considerations 

• Adherence to display design principles 

• Operational considerations 
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Primary consideration should be given to system performance and control task work- 
load, since these two quantities are the ones for which the most empirical data exist. 
Furthermore, the pilot will attend to his control tasks in a hierarchal ordering of 
flight control (flight director), tasks and he will then attend to the monitoring tasks 
if there is available time. 

The specification of the system performance level will have an impact on 
whether or not an automatic performance assessment and failure monitor needs to be in- 
cluded in the system for a given set of trajectory and wind conditions. If the system 
performance is fixed at a low enough level, then the amount of attention required to 

perform the control task f , will be close to the pilot's. ultimate capacity, hence little, 

c 

if any, margin is allowed for unexpected distractions. or for monitoring. If these con- 
ditions prevail then an automatic failure monitor and performance assessment system 
must be considered to accomplish these approaches with the desired degree of pilot 
acceptance. 
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SECTION 3 


DESIGN OF CANDIDATE CONTROL AND DISPLAY SYSTEMS 

Realistic helicopter control system models are required in conjunction with 
the analytic representations of the human pilot and the aircraft displays. In accordance 
with the study guidelines the control system models which were developed represented 
levels of automation that varied over a wide range extending from 1) a fully automatic 
system with the pilot in a passive mode with respect to control activity to 2) a system 
with full manual control. This section includes a description of the different levels of 
control automation investigated and the models that were formulated, and a discussion 
of the flight director design process that was followed. 

3.1 LEVELS OF CONTROL AUTOMATION 

The potential control system configurations can be classified in terms of 
1) the aircraft's three dynamic axes and 2) the outermost feedback loop closure. This 
can be seen in Figure 3-1 which shows the level of control automation for the pitch 
or roll channel. With the position feedback loop closed In a particular axis, the system 
is fully automatic and the pilot becomes simply a monitor for that axis. Completely 
automatic control would involve position feedback in all axes. For normal operation, 
the aircraft would have a stability augmentation system (SAS) with attitude rate feed- 
back. However in the event of an SAS failure, the unaugmented configuration would 
exist. 

The output of the control system in each case consists of the four actuator 
displacements that are used to control the unaugmented vehicle. These consist of 

= pitch control (elevator or differential collective) 

6^ = roll control (aileron or roll cyclic) 
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Pifch or Roll Channel 



Aircraff 

Control 


Unaug- 
men ted 
VTOL 


Figure 3-1 . Levels of Control Automation for Pitch or Roll Channel 







6^ = yaw control (rudder or yaw cyclic) 

= vertical control (collective) 

The input commands to the control system ore to be selected. Since it may be desirable 
to design different levels of automation for each control channel, they were considered 
separately. 

For the pitch and roil channels there are four possible input levels for each 
as indicated In Figure 3-1. These are termed horizontal position commands, horizontal 
velocity commands, attitude commands and attitude rate commands, respectively. 

For the vertical channel there are only two input levels to be considered, 
since vertical translation acceleration Is related directly to collective. For the yaw 
channel there are also only two input levels to be considered because yaw motion is 
not directly coupled to translation as are pitch and roll. Yaw control at the attitude 
rate command level is usually implemented such that yaw rate is proportional to yaw 
control displacement. Yaw control at the attitude level is usually implemented in 
either a "Heading Hold" or a "Turn Following" mode. The former is used during low 
speed and hover. The latter is used at higher speed to give coordinated turns. The 
"Heading Hold" mode is implemented by commanding present heading so long as the 
rudders are neutral. When the rudders are out of neutral, the commanded heading is 
changed at a rate proportional to yaw control displacement. The "Turn Following" 
mode is usually implemented by nulling either the sideslip angle, P, or y-axis specific 
force as measured by a body-mounted accelerometer. 

Considering the unaugmented actuator as an additional mode in each channel, 
there are three or five levels of automation possible In each of the control channels of 
the the helicopter (see Table 3-1). Thus the number of possible combinations is 
5x5x3x3 = 225. However, many of these combinations are not practical systems 
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Table 3-1 . Levels of Control Automation. 
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for normal operations. A series of eight systems were selected to represent the full 
range of automation for the CH-47 helicopter ranging from purely manual with direct, 
actuator commands to full position control. These are shown in Table 3-2. 


Table 3-2. CH-47 Automation Levels. 


System 

Control Channel Command 

Pitch or 
Forward 

Vertical 

Roll or 
Lateral 

Yaw or 
Directional 
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3.2 


UNCOUPLED CONTROL SYSTEM MODELS 


This subsection discusses the approximate closed-loop dynamic response of 
the control system for each control channel. The intent is to provide an understanding 
of the design objectives for each level of automation and their overall hierarchy. 
Subsection 3.3 will present the unified design approach for the coupled systems using 
the vehicle dynamics models. 

3.2.1 PITCH AND ROLL CHANNEL MODELS 

The pitch and roll channels can be made to respond similarly in the low speed 
or hover mode. In the high speed mode air speed is substituted for forward ground 
velocity at the velocity command level and heading rate is substituted for lateral ground 
velocity. The switch from hover mode to high speed mode presents several complex 
design choices. For simplicity it will be assumed that those switches are accomplished 
manually by the pilot at an airspeed between 40 and 80 knots. 

The pitch or roll attitude rate model is shown in Figure 3-2. The vehicle 
response to a pitch or roll rate command (q^ or 0^) is approximately that of a first-order 
system: 

-3_ ^ — 1 (3-1) 

q T s + 1 

^c q 

where the time constant is typically about 0.5 second. The maximum pitch or roll 
rate for the CH-47 has been specified by NASA at 25 deg/sec; treating this as a step 
input to Equation (3-1), the maximum pitch or roll acceleration is 

2 2 
q =50 deg/sec = 0.87 rad/sec . 

^max 
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Figure 3-2. Pitch or Roll Attitude Rote Model. 


Figure 3-3 shows the approximate model for pitch or roll attitude control. 
This system merely involves the addition of a pitch or roll attitude feedback loop around 
the attitude rate model. The closed-loop model response is second-order: 


_e_ 

e. 


2 

‘“n 


2 

s 


+ 2C(u 


s + cu 
n 


2 

n 


(3-2) 


with a riatural frequency co^ = K./t and damping ratio ^ = 1/2 V"kTt For a 

n 0 C| 0 C| 


-1 


typical system gain K- ~ 0.5 -1.0 sec , the response has nearly critical damping 

0 

(q = 0.7 - 1 .0) with a natural frequency u)^ = 1 .0 - 1 .4 rad/sec . 


e 

c 



Figure 3-3. Pitch or Roll Attitude Model. 
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The next level in low speed pitch and roll automation is the horizontal 

velocity control model shown In Figure 3-4. Again, this system essentially involves 

an outer feedback loop around the pitch or roll attitude control system of Figure 3-3. 

The pilot input V is compared with the measured velocity V to generate a commanded 

pitch or roll attitude; V provides a trim capability. The closed-loop transfer func- 

*^ref 

tion, using Equation (3-2) for the inner attitude loop, is third order: 


_V 

V. 


gK\/"n 


s^ + 2ctu^s2+u)2s+gK^u,2 


(3-3) 


This can be simplified at frequencies below tu by approximating 0/0 ~ 1, which 

n c 

results in the closed-loop response 


V 


1 


V 


+ 1 


(3-4) 


where = l/g Ky. If we select a gain of Ky - 0.9 deg/ft/sec, the time constant 
becomes Ty ~ 2 sec. 


'ref 



Figure 3-4. Horizontal Velocity Model - Hover Mode. 
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The highest level of automation in the horizontal channels at low speed is 
the hover position model, shown in Figure 3-5, The position feedback loop is closed 
around the velocity control model of Figure 3-4; changes to the reference Hover position 
X are fed forward through a differentiating network to generate the reference velocity 


signal. The closed-loop response of this system is 


X 

X 


s/K^ + 1 


(3-5) 


c TsVKy + s/K^+1 


Selecting the position gain at ~ 0.25 ft/sec/ft, gives a critically damped system 
(q *-0.7) with a natural frequency of ~ 0.35 rad/sec. 

In the higher speed flight regime, the only change in the pitch channel 
models is that forward ground speed V is replaced by airspeed as shown in Figure 3-6. 
The roll channel in the high speed mode is used to control heading rate, as discussed in 
the next subsection. 

3.2.2 YAW CHANNEL MODELS 


At low speed, the yaw channel is in a heading hold mode, as shown in Figure 
3-7. The closed- loop response is given by 

. K 1 

jl = i 

'^c s^ + s/t . + K . T 


(3-6) 


Typical values are K . ~ 2 sec and t, ~ 0.5 sec, which result in a critically damped 

system (£ = 0.7) with a natural frequency u)^ ~ 1 .4 rad/sec. The maximum yaw rate 

command is about ~ ±25 deg/sec. 

max 

At high speed, the yaw and roll channels operate in a turn following mode 
with heading rate control, as shown in Figures 3-8 and 3-9. Normally, there is no 
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Figure 3-8. Turn Following Model - High Speed Mode. 



Figure 3-9. Heading Rate Model - High Speed Mode. 


command input to Figure 3-8, and the system maintains zero sideslip. If the pilot wants 
to slip the vehicle, the input is proportional to sideslip velocity, PV, with a maximum 
value of approximately 60 ft/sec. The closed-loop sideslip response is 

(3-7) 

(PV^c s^ + s/xp+KpV 

^'Typical parameter values are Kp ~ 0.004 rad/sec/ft, and Xp~ 1.4 sec. The resulting 
dynamics depend on airspeed, as shown by the following: 
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V 


c 


U) 


60 kt 
120 kt 


0.6 rad/ sec 
0.9 rad/sec 


0.6 

0.4 


The heading rate loop in Figure 3-9 is nearly identical to the roll attitude 
model in Figure 3-3, since yaw rate is proportional to roll attitude at a given airspeed. 

The closed-loop dynamics are the same as Equation (3-2), except the natural Frequency 

2 -1 
0)^ = gK 0/VTp and damping of C=l/(2TpU)^). For Xp ~ 0.5 sec and ~ 5 sec , the 

resulting dynamics are also functions of airspeed; 


V 

60 kt 
120 kt 


1.8 rad/sec 
1 .3 rad/sec 


0.6 

0.8 


3.2.3 


VERTICAL CHANNEL MODELS 


The vertical speed model is a simple first-order response, as shown in 
Figure 3-10, with a time constant of xj^ ~ 2.5 sec. The maximum input is of the order 
of h ~ 16 ft/sec, which limits the vertical acceleration response to about 0.1 g 


max 


(3 ft/sec). 



Figure 3-10. Altitude Rate Model. 
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The alMtude hold model merely odds on outer loop around the altitude rate 
model, as shown in Figure 3-11. The closed-loop response is second-order; 




‘/■"h 


-1 . 


(3-8) 


Selecting the gain "• 0.2 sec gives a damping of C ~ 0.7 with a natural frequency 


of u) - 0.25 rad/sec . 
n 


3.2.4 BANDWIDTH SUMMARY 

The bandwidths of the uncoupled models in the various modes discussed 
previously are summarized In Table 3-3. 



Figure3-11. Altitude Hold Model. 


Table 3-3. Approximate Bandwidth of Uncoupled Control Model Responses. 


Channel 

Mode 

Rad/Sec 

Horizontal 

Attitude Rate 

2.0 


Attitude 

1.0 - 1.4 


Velocity 

0.5 


Position 

0.35 

Yaw 

Turn Following 

0.6 @ 60 kt 
0.9 @ 120 kt 


Heading Rate 

1.8@60kt 
1.3 @ 120 kt 


Heading Hold 

1.4 

Vertical 

Velocity 

0.4 


Position 

0.25 
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3.3 


CONTROL SYSTEM DESIGN USING QUADRATIC SYNTHESIS 


This section discusses the design of control systems for the CH-47 which 
minimizes cross coupling between controlled axes. From a theoretical viewpoint it can 
be shown that there are insufficient degrees of freedom available to completely elimi- 
nateall cross coupling. The problem then is to use the available degrees of freedom in 
some optimum manner while constraining the direct control to respond with dynamics 
previously determined to be desirable. 

In this analysis the control was found using quadratic synthesis. The cost 
function was established by using those weightings on the state variables that would 
produce the model response in each channel if it were uncoupled. The control deter-^ 
mined will automatically suppress the cross coupling in the optimum manner for the cost 
function selected. Using quadratic synthesis the cost function was determined by the 
maximum acceptable magnitudes of the control and state variables. The maximum 
values of available control are known for the vehicle, while the maximum values for the 
state disturbances are determined indirectly. The ratio of the maximum value of control 
and state in a direct channel is directly related to the bandwidth that results for that 
channel. Since the desired bandwidth has been determined for all the direct channels, 
these numbers can be used to generate the maximum values for each state variable. The 
philosophy used was that every control system is optimum for some cost function. The 
cost function for which the optimum control is the same control desired in each direct 
channel assuming no cross coupling is first obtained; then the same cost function is used 
to optimize the coupled system. 

3.3.1 QUADRATIC SYNTHESIS SUMMARY 

The steady-state quadratic synthesis design procedure provides a convenient 
method for determining the feedback law to minimize the control performance index 
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/-CO 

J n 


(x'Ax + u'Bu) dt- 


(3-9) 


where x is the plant state vector and u is the input vector. The diagohal weighting 
matrices A and B are determined from the maximum desirable variations in the states 
and controls, i.e.. 


A.. = l/(x. ) 

II ' I ' 

max 


(3-10) 


Bji = l/(uj ) 

max 


2 


(3-11) 


The linear system dynamics are given by the standard form 
X = Fx + G(u + u ) 

- - •■ . P . • . . . 

where u^ is the additional input from the pilot. The feedback control law is 
'■ u = -Kx ■" 


where the feedback gain matrix is given by 


K = B"^G'S 


The symmetric matrix S is the steady-state solution to the Ricatti equation 
S = -SF-F'S + SGB"^G'S-A = 0 


(3-12) 


(3-13) 


(3H4)- 


(3-15) 


The closed-loop system dynamics are found from substituting Equation (3-13) 
into Equation (3-12): 

X = (F - GK)x + G u (3-16) 

P 

= F*x + G u (3-1 6a) 


where F* = (F - GK) is the closed-loop system matrix. 
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Occasionally, it will be useful to consider one or more of the pilot inputs as 
a commanded value for an element of the state vector. For example, if it is desired 


that pilot input u correspond to direct control of state variable x. , we can use 
Pj I 

Equation (3-13) to obtain the transformation 


u = -K.. X. 
P; M I 


(3-17) 


3.3.2 UNCOUPLED EXAMPLES 

Several simple uncoupled examples are presented briefly to illustrate the 
design procedure and the relationship of the desired system bandwidth. 


3.3.2. 1 PITCH RATE CONTROLLER 


Consider the quadratic synthesis design of the pitch rate controller shown in 


Figure 3-12, 



Figure 3-12. Pitch Rate Controller. 


The open-loop system is defined by 

X = q, u = q, F = 0, G = 1 
The performance index weightings are scalars: 


A = B = — 

2 -2 


The feedback control gain is 
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K = B‘^G^S = q^S 
*m 


where S, the Ricatti solution is easily found 

1 


Hence, 


S = 




K = !il 


'm 


or 


T = 


’m 


’m 


In this first-order system the system bandwidth is directly related to the ratio 
^m^'^m* ^m from physical constraints, the desired bandwidth can be 

achieved by properly selecting in the design process. 

3. 3. 2. 2 PITCH ATTITUDE CONTROLLER 


The same technique is applied to the pitchattitude controller shown in Figure 3-13. 



Figure 3-13. Pitch Controller. 
The natural frequency and damping of this system are 
'"n " ^2 

C = 1/2VK/K2 


3-16 


AEROSPACE BY8TEM81 IIMC* • one vine brook park • Burlington, in/iasbachusetts oieos • (ei7) s7b-7bi7 





The state vector and the open loop control are 


X = 


q 

e 


, u = q 


and the system matrices are 



o 

o 


T 

F = 


G = 



.1 0_ 


0. 


The performance index weightings are 

"0 0 


A = 


,0 1/eU 


m 


B = 1/q 


m 


where 6, q , 0 are the maximum permissible values of pitch acceleration, rate, 
and attitude, respectively. 

The quadratic synthesis feedback controller is 


q 


u = -K 


e 


where K = [K K„] = [K, K.K^]. 

q 9 I I / 

The Ricatti equation for this system can be solved explicitly to find the 
feedback gains 



Thus the quadratic synthesis design leads to the following dynamic response 
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c = -Jl/l = 0.7 
% = 

Therefore, the only parameter needed to design the system using quadratic 
synthesis is the ratio 9^9 ~ ‘“n* been specified on q as well as 9, 

only the damping in the system would have changed; the bandwidth would still be 
determined by 

3. 3. 2. 3 HORIZONTAL VELOCITY CONTROLLER 

As a final example, consider the horizontal velocity controller designed 
using quadratic synthesis (Figure 3-14). 



Figure 3-14. Horizontal Velocity Controller. 

The open-loop system is defined by 

‘0 0 gi ro 

,u = q,F = 000 ,G = 1 

.0 1 oj Lo 

The performance index weights are 

"1/V^ 0 O' 

m ^ 

A = 0 0 0 B = 1/q^ 

m 

_0 0 0 . 
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The feedback controller is 


u 


-K 


"V" 

q 

_e . 


where K = Kg] = [K^K 2 K 3 , K^, K^Kj). 

Solving the Ricatti equation to obtain Ky, gives 

K., = K,K,K, = q /V 
V 12 3 m 

which could have been deduced intuitively. It is not necessary to solve for the other 
terms since we have sufficient information to determine the ratio. 

There is now sufficient information to use the quadratic synthesis technique 
to determine the feedback control for the coupled equations. 

3.3.3 COUPLED LONGITUDINAL CONTROL 


For the coupled longitudinal helicopter control system design, the state and 
control vectors are defined by 



forward position 
vertical position 
pitch angle 
forward velocity 

vertical velocity 

pitch rate 



differential collective input 
gang collective input 


(3-18) 


(3-19) 


The system matrices F and G describing the kinematic and aerodynamic response of the 
rotorcraft are developed in Appendix B. Numerical values of the CH-47 stability 
derivatives for various flight conditions are presented in Appendix E. 


3-19 


AEROSPACE BVSTEIVie, INC. • one vine bbook park • buhlinston. Massachusetts meo3 . (si7) 



The control weighting matrix in the performance index is 

v-2 n 


B = 


(6e ) 


max 


LO 




(3-20) 


The control input limits 6a and 6^ can be determined approximately for each 
flight condition from the constraints on vehicle angular and vertical accelerations 


’max 


e 


max 


/I 

^ yy 


(3-21) 


max 


w 

max 
Z, /m 


(3-22) 


where M, /l and Z /m are the stability derivatives explaining pitching acceleration 
®e ®c 

due to 6g and vertical acceleration due to 6^/ respectively. From the models in 
Subsection 3.2, typical limits on the accelerations were found to be 


q ~ 0.87 rad/sec^ 
^max 


w ~ 3 ft/sec^ 


max 


The state weighting matrix A takes on different forms depending on the level 
of automation, as shown by the following examples. 
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Pitch Rate and Vet tied I Velocity Comnrtond 



• Pitch Attitude and Vertical Velocity Command 





Forward Velocity and Vertical Position Coinmand 



The maximum values of the states to be used in the "A" weighting matrix can 


be determined from the uncoupled model bandwidths in Table 3-3, and the acceleration 
limits described above: For the horizontal channel of the CH-47: 
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X /u ~ 3.0 sec 
max max 


^max 


max 


max 


max 


0.435 rad/sec 
0.435 rad 

28.0 ft/sec 

84.0 ft 


For the vertical channel of the CH-47; 

w /w ~ 2.5 sec 

max max 

z /w ~ 4.0 sec 

max max 


V « w - 7.5 ft/sec 
z max ' 

max 


z = 30.0 ft 
max 


These design parameters and the quadratic synthesis technique v/ere used to 
generate the closed-loop system dynamics for the seven different levels of longitudinal 
control automation presented in Section 6. 

3.3.4 COUPLED LATERAL CONTROL 

The state and control vectors for the coupled lateral helicopter control system 
design are defined below: 


X = 


y 

0 

' 1 ' 

V, 

p 

r 


u - 


L r 


lateral position 
roll angle 
yaw angle 
lateral velocity 

roll rate 
yaw rate 

roll cyclic input 
yaw cyclic input 


(3-23) 


(3-24) 


The system matrices F and G for the lateral modes are also developed in Appendix B, 
and the stability derivatives for the CH-47 are given in Appendix E. 
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The control weighting matrix is 




('r ) 

max 


(3-25) 


In this cose the maximum control inputs are again related to the angular acceleration 
limits: 


max L. /I 

6/ XX 
a 


(3-26) 


"^max N. /l 

6 / 22 


(3-27) 


From the uncoupled models in Subsection 3.2, the values of the lateral accelerations are 


p =0.87 rad/sec'' 
■^max ' 


r =0.87 rad/sec^ 
max 


As before, the A weighting matrix takes on different forrns depending on the 


level of control automation. Examples are 


Roll Rate and Yaw Rate Command 
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Roll Angle and Yaw Rate Command 


• Lateral Velocity and Yaw Angle Command 

~0 0 0 0 0 o“ 

0 0 0 0 0 0 

0 0 0 0 0 

A = max 

0 0 0 v”^ 0 0 

y 

'max 

0 0 0 0 0 0 

_0 0 0 0 0 o_ 

Maximum values for the state weightings again can be determined from the 
uncoupled models and the acceleration limits above. For the CH-47, these are 


r /r ~ 0.5 sec 

max max 


r = 

max 

0.435 rad/sec 

^ yr 

^ max max 

~ 0.7 sec 


^ = 

^max 

0.30 rad 

^max^Pmax 

~ 0.5 sec 

- 

P “ 

'^max 

0.435 rad/sec 

^max^Pmax 

~ 1 .0 sec 

- 

0 

^max 

0.435 rad 

'^max^^^max 

... 2.0 sec 

- 

V 

^max 

V = 28.0 

max 


~ 3.0 sec 


y = 

'max 

84.0 ft 


3 
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3.4 FLIGHT DIRECTOR DESIGN USING QUADRATIC SYNTHESIS 

This section describes a procedure for using the quadratic synthesis technique 
to design flight director signals for the optimal control pilot model. 


3.4.1 DESIGN OBJECTIVES 

The basic concept associated with the use of a flight director is to provide 
the pilot information that is useful for control, thus rendering the piloting task easier 
in some sense. The general form of the flight director signal is q linear combination of 
vehicle states 


FD = h‘x(t) (3-28) 

t 

with possibly some filtering to remove high frequency components. The flight director 
gains h are chosen so that if FD(t) is kept "small" the resulting aircraft motion is desirable. 
Since the pilot is still in the loop, there are two issues that relate to the harmony between 
FD(t) and pilot response. The first concerns the nature of the control task as viewed by 
the pilot. Thus, the task of keeping FD(t) small should not conflict with the overall 
control task requirements. In the optimal control pilot model the latter is manifested in 
the choice of a control input, u, that minimizes the quadratic cost functional 



(3-29) 


The second issue relates to the required form of the pilot compensation. It 
is well known from experimental results in manual control that one of the easiest control 
tasks Is associated with K/s dynamics, wherein the pilot acts (approximately) like 

H (s) = ^ (3-30) 

T^s + 1 y(s) 

^Note that pilot control is not added directly into the flight director signal. 
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Thus, from a workload point of view, one should design a flight director signal FD(t) 
where the resulting pilot response is 

u(s) = H^(s) . FD(s) = H^(s) . h' x(s) (3-31) 

This is tantamount to picking h so that the open loop transfer function 

V(s) = h' (si - A)'^B (3-32) 

is approximated by K/s. 

3.4.2 SIGNAL GENERATION 

The above concepts and objectives have been recognized in the literature 
(References 8-10). Within the framework of the optimal control model. Reference 9 
has suggested constructing FD(s) as 

FD(s)= 2 hiWy.W (3-33) 

1 = 1 . 

where y.(s) are the displayed outputs and h.(s) are the internal transfer functions gener- 
ated by the optimal control model (Reference 11). This approach has FD(s) compatible 
with pilot (subjective) control requirements, and has inherent filtering built into the 
generated signals via the h.(s). However, only the system outputs are included in FD, 
and the computation of the h. require arbitrary choices for pilot observation noises, 
thresholds, etc. 

The approach that we follow in selecting a flight director signal to be asso- 
ciated with a given control is to recognize that the "optimal" control generated by the 
pilot model is (see Appendix A), 

u. + u. = -j^! X (t) = u^(t) . (3-34) 
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Thus, associating the lag ( s + 1) ^ with the pilot, it is seen that a flight director 
signal 

FD. = X (t) = u Jt) (3-35) 

is precisely the control that the pilot model would apply under ideal conditions.^ Sum- 
marizing, the design of a flight director associated with control i is as follows; 

• Select the nominal cost functional weighting associated 
with the subjective task requirement?, 

• Compute the feedback gains 

• Retain only the important elements in i,. to simplify 
implementation. 

The above approach is simple, is related to the pilot's interpretation of the 
task, and assures that the pilot's transfer function between FD. and control u. is ap- 
proximately as given above. 

3.4.3 VALIDATION OF DESIGN MODEL 


To test the above technique, we perform a design for the pitch axis and the 
power (collective) axis, for the CH-46 in hover mode. The results are compared to the 
flight director described in Reference 12. The cost functional weightings for this task are 
selected as (see Section 4) 



(3-36) 


where the corresponding maxima are 


^This approach is valid if the pilot time^elay tIs negligible with respect to system 
time-constants. 
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AX = 25 


> 1 bo ff/sec 
t 2 = 5 = 2 ft/sec 

0. =2° q = l°/sec 

The control rote weightings ore selected to give t ^ Q.l withq. = 0.25 (in/sec)~^. 

rs: ^u. . . ' ' ' 

• I, min 

The important feedback gains in each of the flight director channels are 
FD =C x(t) + C*V(t) 

FD = C z(t) + C- V (t) 
z z ' ' z z' ' 

Since the overall scale factors of these signals are arbitrary (depending on display gain) 
we normalize the signals by the positional gain. Thisprovides an intuitive feel for the 
ensuing time response as well, since if the flight director signal is zeroed, then the 
system response is exponential with time constant or The normalized 

values of gain ratios found via the optimal control model are 

C. C. 

= 10.3, =2.51 

C C 

X z 

For the flight director design given in Reference 12, which was obtained through extensive 
flight testing, the gain ratio C^/C^ = 10.0 which is in excellent agreement with the 
above result. 

The results for the power flight director command are somewhat different 
though and require further analysis. The form of the flight director given in Reference 12 

FD^ = oz(t) + pz(t) + Yz (f) 

where <y, p, and y are weighting coefficients, and the velocity/displacement weightings 
are given by 
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-^ = 10 sec 
a 

If the values of the accelerometer signal are expressed in terms of the state variables 
through the state variable equation, we obtain 
z = t-z(t) + v6g(t) 

which becomes 

FD =Q'z(t) + (p + Yt)z(t)+Yv6 (t) 

^ c 

Using the numerical values for the stability derivatives at hover, we find that the ratio 
of actual velocity gain to position gain is given by 

=2.64 sec 

Of 

which is in very close agreement with the value obtained above. This gives further con- 
fidence to our flight director signal design technique. 

3.4.4 MODELING THE FLIGHT DIRECTOR EFFECTS 

In order to include the flight directors within the framework of the optimal 
control model it is necessary to select threshold values and cost functional weightings. 

If instruments have been specified (i.e. format level) then thresholds can be chosen 

V 

based on eye physiological considerations. At the element or informational levels, 
thresholds and/or maximum excursions must be chosen via alternate means. 

Before we select cost functional weightings it is appropriate to consider 
whether indeed the signals FD. should even be included in the cost functional . This 
issue has been raised in Reference 9, in a somewhat broader context. Basically, one 
must decide on the nature of the pilot's control strategy. Excluding the FD. from the 
cost functional implies that the pilot's control objectives are basically the same as 
before introducing these signals, i.e., in terms of the situation variables Ax, Ay, etc.. 
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and the flight directors provide only for enhanced state information. Including tf)e FD. 
within the cost functional, in addition to the other terms, implies that one of the pilot's 
direct control objectives is to keep the FD. small. 

In our work we will interpret the piloting task to be the latter. Ouf method 
for choosing the associated cost functional weightings on the FD. is to pick 

2 


\ I, max/ 


(3-37) 


where we determine FD. by 

I, max ' 


FD 


I, max 


Jij'x(t) 


X = X 


(3-38) 


max 


for displacement variables i.e. we substitute into the expression for FDj the maximum 

subjective values for the positional variables, Ax, Az, 9 . The rate variable terms are set 

to zero when evaluating the FD. . The threshold values on the flight director dis-- 

\ fTt\QX 

placements are chosen as 


1 


apr* = -(FD.) 
FD. A I 


max 


(3-39) 


and the corresponding rate thresholds are picked as 


°f6 


1 


i' 2 °FDj . 


(3-40) 
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SECTION 4 


CONTROL THEORETIC MODELS FOR PREDICTING PILOT 
MONITORING PERFORMANCE 

A very important aspect of evaluating different display/automatic systems is 
the ability to predict human monitoring behavior, and to give a metric for assessing 
monitoring performance. With the optimal control model, we have already obtained 
such predictive capabilities for human control response. This section examines different 
approaches for monitoring prediction and choices of monitoring workload performance 
metrics. 


4.1 


CHARACTERISTICS OF MONITORING MODELS 


The approach for determining display requirements and for evaluating systems 
of differing automation levels is discussed in Section 2. Basically, the total workload 
level fjQi- 'S selected; then, for a given automation/display system configuration, the 
fraction of control attention f^ ^ ^TOT required to achieve a desired performance 

level is determined. The excess capacity, fjQj " ~ ^m available to the pilot 

for monitoring the displays. This approach, in which control performance is established 
first, was followed since it a priori limits consideration to those systems that have 
realistic requirements at the initial stages of investigation. Thus, for systems in which 
f >0, the objective is to determine how this monitoring workload is allocated among 


m 


the n displays; i.e., to determine the f ,i-1,...,n. 
y m. y 

The monitoring models that are considered have all been evaluated in 

the light of certain desirable characteristics. These are determined by the overall 

goals of the monitoring process as used by pilots: 

1 . Assess aircraft situation with respect to mission requirements, by 

minimizing the (relative) estimation errors associated with variables y. . 
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2. Assess operation of aircraft displays and controls by cross-checking 
displays for redundancy and self-consistency to detect system failures. 

The first goal is a statement of the well-known fact that pilots desire situation 
information in addition to the control information provided by the flight director and 
other tracking aids. The second goal indicates that the pilots will continue to scan every 
available instrument which is important to the approach, even though straightforward 
application of optimal estimation/control models of the pilot would indicate that some 
instruments may not be scanned because of the correlated information among displayed 
variables. The goal of this scanning process may be thought of as the detection of 
instrument and control system failures. 

Thus, the important characteristics of the monitoring models are assumed to be: 

1 . The ability to assess the performance of the aircraft from situation 
displays. Implicitly derived variable rates are not monitored, but 
their information is used to obtain better estimates of the explicitly 
presented display variables. 

2. No monitoring of flight directors, or other combined state information 

that is geared to aircraft control. Thus, for these instruments, f =0 
sothatf^Q^^ = ^. "^i 

3. Required monitorinq of all primary status instruments. Thus, 

f ^€>0. 

m. 

I 

4.2 MONITORING MODELS 

The two major components of a monitoring model are the monitoring per- 
formance metric (the method by which monitoring performance is evaluated), and the 
attentional allocation scheme (the manner in which the model allocates the attention 
among the various displays). Ideally, the two components can be combined by requiring 

that the attention allocations f be chosen to minimize the given metric subject to the 

m. ^ ' 

I 

constraint T'f = f and that each f be greater than some specified value. However, 
^ m. m m. ^ ^ ' 

i I . ' . . 

it is possible to choose the f^ according to criteria other than optimizing a performance 
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metric. Several methods for choosing the f ore discussed below. In order to keep the 

nij 

presentation simple, these assume that the pilot is monitoring an automatically controlled 

system so that each f =0 and thus frr\T ~ f • Situations in which the aircraft is under 
c. I w I j m* 

manual control, i.e,, f >0 are considered in Subsection 4.3, wherein the f i$ added 

c. m. 



4.2. 1 NON-METRIC BASED MODELS 

4. 2. 1.1 EQUAL ATTENTION 

Choosing the attentional allocation according to 


f =f /n 


m 


m y 


(4-1) 


while simple, does not take into account the relative importance of instruments, nor 
their correlations. However, it does assure that all instruments will be scanned for 
failure detection . 


4.2. 1 .2 PEAK EXCURSION MONITORING 


It is reasonable to expect that a pilot will monitor a signal when its value 

exceeds some multiple, P, of its standard deviation. Thus, usual signal levels are not 

of immediate concern, but signal values greater than usual are monitored. The choice 

of f becomes 
m. 

I 


f =f Pr|ly.|>Pa } 
m. o- ( I ' y./ 


(4-2) 


This technique assumes that the pilot will look at variable y- with probability 

f whenever |y. exceeds . P is a parameter chosen so that/^f = f . However, 
o. y- ^ m. m 


I 


m. 

I 'I II 

the f remain to be chosen in some plausible manner. This scheme introduces the con- 
o. ^ ^ 

I 

cept that monitoring is dependent on relative signal values (i.e., on the ratio of 
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ly. [/^ where is the standard deviation of /.)• Conceptually, the f should be 
I y. y. I o. 

related to the "importance", or information content of the display. However, a pilot 
may not necessarily look at a variable just because ly. 1 if 

1. He can acquire this Information from a second display highly correlated 
with y. . 

2. His estimation error in y. is sufficiently small so that looking at y, serves 

no purpose except instrument verification. ' 

3. Both y. and cr. are less than the display visual or indifference thresholds 
(i.e., the display may be poor). 

Nevertheless, this model with P = 1 .5-3 means that the human will monitor 

the unusual occurrences. In fact, this model with f = 1 and P = 3 requires a minimal 

o • 

\ 

or residual attention be placed on each display of f ~ .01; if P = 2, each f ^ .05. 
By interpreting this to mean a monitoring of unusual or failed situations, this model 
handily provides a lower bound for each f 

^i 

4. 2. 1.3 NYQUIST CRITERIA MODELS AND THEIR EXTENSIONS 


As first postulated in Reference 13, a human samples an instrument 
periodically in an attempt to reconstruct the associated time signal. Thus, information- 
theoretic ideas, particularly Shannon's sampling theorem, were used to obtain the ex- 
pression 

A. 

f = 2c, w. Iog_ — — + 2<A). c- (4-3) 

m. I I 2 p I 2 

1 c • 

I 

where ^ | = signal bandwidth, A. = signal RMS amplitude and E, = permissible rms 

error, c^ and C 2 are constants, the latter used to account for minimum fixation time. 

For multiple instruments where the ratio of signal power to magnitude of significant 

deviations E. are roughly constant, the f would be proportional to signal bandwidths. 

i 

^One might conjecture that f > .05 will enable the detection of system failures in 
short time. "^i 
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(4-4) 


Thus# since 

I 


is proportional to 



witfi 


f 

m. 


= ko, /o 
V{ Xi 



(4-5) 


This simple periodic sampling model would not adequately predict behavior 
in more complex situations with correlated signals and aperiodic sampling behavior. 
Noting that pilots are often concerned only with detection of extreme readings rather 
than with signal reconstruction. Reference 14 proposed a conditional sampling scheme 
that would result in aperiodic behavior. In this approach, the human is considered as 
a channel for the transmission of discrete messages in lieu of a complete time function. 
In this context it is possible to postulate several (somewhat related) sampling strategies. 
Thus, References 14 and 15 hypothesize a strategy in which a sample Is taken when the 
probability that the signal exceeds a prescribed limit is greater than some subjective 

. t 

probability threshold. Reference 15 assumes that a sample is taken when the probabil- 
ity of exceeding the limit is a maximum. On the other hand. Reference 15 suggest^ a 
sampling strategy based on a "Variable Nyquist Interval" . Unfortunately, none of these 
conditional sampling models have been tested against experimental data nor have they 
achieved a high level of acceptance in the manual control field. 


4.2.1 .4 MONITORING FOR FAILURE ANTICIPATION 

The assumption that a pilot monitors an automatic system in such a manner 
as to anticipate a manual takeover provides a basis for an alternate monitoring model . 
For a given failure, one can solve the control problem associated with the pilot control - 


^'The form of this model is somewhat similar to that of Subsection 4.2.1 .2. 


4-5 


AERaaPAce avsTEMa, inc. • ONE VINE BF^OOK PARK • BURLINGTON, MASSACHUBETTe 01S03 • (617) S7S-7917 



ling fhe failed dynamics subiecf to an attentional level of Optimizing the 

control cost for the f will provide the monitoring fractions f by equivalence. 
c» nrij 

Although this method has some intuitive appeal, it suffers drawbacks in 

that 

1 . A specific anticipated failure mode and dynamics must be assumed. 

2, If there were a failure, the f would probably not equal f so that 

c nn 

the f as found may have little relation to the f . 
c. ' m. 

I I 

4.2.2 MONITORING METRICS 

The monitoring models discussed above are not geared. to a performance metric 

with which to evaluate sampling behavior. However, once a monitoring strategy f 

m» 

is determined, any number of metrics may be applied after-the-fact. This suggests an 

alternate, more appealing approach to the monitoring problem: to first specify a 

meaningful performance metric that embodies the. goals of the monitoring function, 

and then to choose the f to minimize the selected metric. Consequently, the dual 

m. 

I 

goals of monitoring as status determination and failure detection suggest the following 
two monitoring cost functionals. . . 

4.2.2. 1 ESTIMATION ERROR COST FUNCTIONAL 

For status determination the pilot's monitoring strategy Is to choose the f^ 
subject to ' 


n 



>0 


to minimize the monitoring cost 


n 




o2 

^i 


(4-6) 


(4-7) 
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where is the rms estimation error in monitoring signal y. . The Yj are scale factors 
that are either 0 or 1 to indicate whether an instrument is of monitoring concern. Ob- 
serve that this cost functional has the properties 

1 . J since it is the relative estimation error that Is weighted, 
m ' 


2. Only if f =0 and no information concerning y. is obtained from other 
nn. on o o I 


I 2 2 2 2 

insTrunrients can otherwise o . 

e. y. e. y. 


I ' I 


Another Interpretation of is obtained by defining 


e. 


k. = — 1- - error fraction for variable y. 

I * I 


(4-8) 


so 


that 


1/2 


J = rms monitoring "error fraction" 
m 


The optimal choice of the f provides a prediction of the monitoring fraction for each 

' 'i 


displayed variable . Note that 


1. All instrument correlations are considered in this formulation, so that 
key instruments have larger f 

mi 

2 . . 2 

2. If u << threshold value, the model will not try in vain to lower o , 

y • o c I 

but accept o' ~ cn by keeping f ~ 0. 

®i ^i "^i 

The error fractions k. are useful in relation to the probabilities associated 


with estimation error criteria. By defining 


E(P) = Pr{le(t)l>p0 } 


(4-9) 


as the probabilify fhaf the estimation error at any time exceeds a fraction P of the sig- 
nal rmsy then (assuming Gaussian statistics) 
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Making the change of variable w = ^~2 yields, noting k = 



(4-10) 


( 4-10 


A reasonable performance level to expect in monitoring is 0 = 1/2, i.e., the estimation 
error should not exceed o/2 of the monitored variable. The percent of time (or 
probability) that this criteria is exceeded is E(1/2), and is monotonically related to the 
error fraction k. This gives another interpretation to the cost functional (Equation 
( 4 - 7 )). 


4.2.2. 2 FAILURE DETECTION COST FUNCTIONAL 

A recent decision model based on Wald's sequential analysis (Reference 16) 
provides an excellent description of the human's ability to detect failures in random 
processes. A single-channel version of the model was successfully validated with 
experimental data. A more elaborate test of the model was made by applying it to the 
task of monitoring a fully automatic ILS approach, with equally encouraging results 
in describing the experimental data. 

The model for monitoring a single instrument uses Wald's Sequential 

Probability Ratio Test to derive a decision function based on the log likelihood ratio 

of each observed residual from the Kalman filter. The decision function ^ at time 

m 

t^ is given by 
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m 


X = 

m 


X > 0 
m 


X < 0 
m 


(4-12) 


w 


here 


X - \ , + 

m m- 1 


r -e/2 

m ^ 


(4-13) 


and r Is the Kalman filter residual at t .o’ is its standard deviation, and 6 is a bias 
m m r 

parameter beyond which the process Is considered to be "foiled". Note that th^ de- 
cision function Is reset to zero should It become negative. 

The decision of "failure" or "no failure" is made from the current value 
of X according to the rule 


m 


If X > — In A, decide "failure" 
m g 


If 0 ^ X < — In A, decide "no failure" 
m 


(4-M) 


The consfanf A is related to the probabilities of the two types of errors: the probability 
of a missed alarm (deciding *'no failure*' when a failure is present); and the 

probability of a false alarm Pp^ (deciding "failure" when no failure is present)^ 

(Reference 4). 

There are two metrics proposed here: the mean time to detect a failure and 
the mean number of looks to detect a failure. These can be determined from the ensemble 
performance of the model of human monitoring and decision making as follows. 

Assume that a constant (bias) failure of magnitude 9 has occurred, I.e., the 

2 

residual r has increased by fl *, and thatall residuals have the same variance of which will 
m ' r 
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be a function of oM the fractions of attention for monitoring ^ f^^ Then the decision 
function is given by 


V +e -e/2' 

\ « X , 

m m-1 


(4-15) 


" -=m a 


(4-16) 


■f r 


where r. is the zero-nrtean component of the residual and is assumed zerp. 

The index m in Equation (4-16) represents the "look nunrtber" sincp a decision cannot 
be made without glancing at tbe instrument, and (m - m^) is the number of looks that 
hove elapsed since the failure occurred. The above approximation is valid when 
is much less than the level of the detection threshold, a reasonable assumption when 
the false alarm rote is not too high. This approximation is also useful in determining 


the mean number of looks to detect the failure since at the moment pf detection 

0 - e /2 


mf-^mD 


A. m pv 
m D 




(4-17) 


a 

i-m^ r 


0 


where iHp Is fhe number of looks to detect the failure. The mean value of m^ is thus 
(approximately) 


In A 


mr 


0 /e* - d/2] 


y2 a2 


(4-18) 


2 . 

where Is the a priori variance of the displayed signal. For the case of multiple 
instruments the mean number of looks to detect failures is therefore 
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(4-19) 


1 / 

J =— i_y Y. 
m n I 


^ ' ^2 
y 1=1 '"v. 


and is a reasonable mefric to minimize by proper choice of the f 

i 

Another important performance metric, is the mean time to detect a failure. 

This can be expressed in terms of the relevant parameters of the optimal estimation model 

th 

as we now show by using a measure of "looks per second". If the i instrument has re- 
ceived n. looks in T seconds, and each look requires seconds = 0,3 - 0,4 seconds 

according to experimental eye movement data), then the fraction of time spent looking 
th 

at the i instrument (dwell fraction) is 


n.A 

dwell fraction = — ! — . 

T 


(4-20) 


Under the assumption that this is close to the fraction of attention f , the scan period 

i 


(seconds/look) can be determined 


J_= __L 

n. f 

1 m. 


(4-21) 


Thus fhe mean time to detect a failure on the i instrument is 


I t 


In A 

e. /e*. -e./2\ f 

mi _L ' 


a \ a 

^i \ yi 


1 

\ 


m. y. 

I ' I 


(4-22) 


so that minimizing the functional 


J 


o2 

Y. r. 

I I 


n .;ri f ct2 

y ' ' m. y. 

' I ' I 


(4-23) 
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minimizes the average mean time to detect failures, as compared to minimizing Equation 
(4-18) which is the mean number of looks to detect failures. The disadvantage of using 
Equation (4-23) is that it does not reflect the differences in bandwidths among the dis- 
played signals. If one postulates, as Senders does, that the look rate is proportional to 
the Nyquist sampling rate, then the number of "looks" required to detect a failure is pro- 
portional to the number of elapsed "cycles" of the displayed signal required to detect the 
failure. This has an obvious intuitive appeal. On the other hand, we suspect that pilots 
may sometimes monitor the lower frequency signals qt a somewhat higher rate because of 
time considerations. One example of this is a relatively smooth, coupled ILS approach in 
which the ILS needles have very low frequency content (on the order of 2 to 10 cycles for 
the entire approach). Detection time would take precedence over detection looks in thjs 
case . 


4.2.3 METRIC BASED MONITORING MODELS 


Having postulated metrics for monitoring performance, one must next consider 

the mathematical problems of minimization with respect to the f . These are con- 

m; 

I 

sidered below. 


4.2.3. 1 RESIDUAL MONITORING FOR FAILURE DETECTION 


In considering the cost functional (Equation (4-19)) within the context of an 
optimal control/estimation model, the Kalman filter residuals are (neglecting the humans 
time delay t), 

= y;(0 - cj x.(t) (4-24) 

•th 

where c| is the i row of the C matrix; it is well known that for optimal filtering the 
residuals are white with covariance equal to the observation noise covariance. Thus, 


^2 _ 2 
= a 


= P.a^ /f 
I y ' m. 


I 


4-12 


(4-25) 


ABROBPACe SYSTEMS, INC. • ONE VINE BROOK PARK • BURLINGTON, MASSACHUSETTS 01803 » (S17) S7S-7S17 



where 


3 =a /N.(cf ) 

y-i y{ y\ 


(4-26) 


and N is the describing function gain for the display variable threshold. 


Substituting into Equation (4-19) gives 


1 y v.P, 


1 


j = — V 

m ^ 

n , r~7 


w • 1' N. f 
y 1=1 I m. 


(4-27) 


This expression is minimized (subject to constraints) at 


*. Vi 


"y Y.p.\-i 


m. 


~ • Z — 

N| Vh ^i 


. f 


m 


(4-28) 


Thus, for well-designed displays (N. = 1, P. = .01 tt)^ each subject to the same 


probability of failure (Yj - Yj), 


f = f /n 


m 


m y 


This provides a new interpretation of the simple, equal attention allocation scheme, 
noted in Subsection 4. 2. 1.1. 

The second failure -related cost functional given by Equation (4-23) can also be 
optimized in a straightforward manner. For the case N. = 1, p. = P., Y. = Y|, an equal 
division of attention is again obtained. This is to be expected, since the cost functionals 
in question have been motivated in terms of (uncorrelated) bias errors on instruments. 

4.2.3. 2 RELATIVE ERROR MINIMIZATION (ASSUMING UNCORRELATED 
INSTRUMENTS) 

As an approach towards minimizing the metric (Equation (4-7)), inter-instrument 
correlations can be neglected, and each instrument can be considered to represent an un- 
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coupled second order system consisting of y. and y.t . The human thus monitors these 
two variables to obtain an estimate of y. with associated estimation error e. . The 
problem is then; "given observations y^(t) = x(t) + Vi(t), y 2 (t) = x(t) + V 2 (t), determine 
the best estimate of x(t)." 

This problem can be solved by designing a Kalman filter for the system 

x(t) = y2(0 - V2(t) (4-29) 

where 

y 2 (t) x(t) + V 2 (t) = velocity measurement 
y 1 (t) = x(t) + v^ (t) = position measurement 
V. = cov (v.(t)] i= 1, 2 
The optimal estimate is 

x(t) = g[y](t) - x(t)] + y 2 (t) (4-30) 

where 

Defining 

= f^2 

9 


Since relations between y. and y. are neglected, information concerning y. can only 
be obtained by viewing y. . 
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yields 


Tx(t) + X(t) = y^(t) + Ty^W 


(4-32) 


as the first-order filter that generates the position estimate x(t). The error covariance Is 


cov[x(t) - x(t)]= = y V|V 


(4-33) 


Note that If Vy -*0 or V2 “* cx>, t -* 0 and x(t) Is obtained directly from y 1 (t) . If V2 0 
or Vi “* 00, g 0 and x(t) Is obtained by Integrating y2(0 . 


For pilot monitoring, 

V. = P. <7?/f.N? 1 = 1,2 

I I r I I ' 


(4-34) 


Thus, for yi and y2 obtained from a single display Indicator 
^2 


_1. = / ^ ' 
2 


X 


1 


N N. 

X X 


(4-35) 


m. 

I 


For the case of multiple Instruments, where It Is assumed that only the position 

(and not the rate) variables are monitored, the monitoring cost functional becomes 

n 2 

1 y ° 

' -4 e: 


m n ;-i ^ 


1 y P. 


= — E 


n . 1 N N. 

y i-l y. y. 


a2 


(4-36) 


y l-l y. y. \ y. m. 


where the sum Is taken over the position variables only. For well -designed displays, 

N = 1, N,. =1 and p. = .01 n and the optimum f are 
y: . y; I f r ^ 


o. 

* y; 


m. 

I 


n 2-1 

y vJ* 

. ( E I -f 


i=i “ 


m 


(4-37) 
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This result/ wherein f is proportional to signal bandwidth/ is precisely the Senders' 

i 

monitoring model arrived at via an entirely different process. 

Despite its simplicity, there are drawbacks with the model as proposed, 

1. Instrument correlations are not treated. 

2. If >> the model would place a high percent of attention on 

instrument i.* In this case x(t) y^(t) +r.|(t)and the white estimation 
error has infinite MSE. This is impossible since we must have 
and suggests that the first order model for x(t). Equation (4-29) may 
need modification. 

3. All of the attendant deficiencies of the original Senders' model resurface. 


4. 2.3.3 RELATIVE ERROR MINIMIZATION (ASSUMING CORRELATED INSTRUMENTS) 

Given the deficiencies of the simple, uncorrelated model we consider 
the entire dynamics of the system being monitored.. Thus, we assume that the pilot is 
monitoring an automatic (stable) system that is driven by white noise 

x(t) = A X (t) + Ew (t) 


y(0“C^(0 = displayed variables (4-38) 

The pilot observes 


yp;(0 = y;(t ~^) + v.(t - t) 


i = l,2. 


. ., n^ (4-39) 


The estimation error covariance for this monitoring task is 


c _ At A't 
S = e Z e 


f 


+ I e^^ EWE'e^'^d? 


(4-40) 


0 


where W = cov(w (t)] and Z satisfies the Ricatti equation 


0 = Z A' + AZ + EWE' - ZC V C2 

y - 


(4-41) 
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and the observation noise covariances are as given in Equations (4r2 5) and (4“26). 

Substituting into Equation (4-7), and neglecting those terms that do not de- 
pend on V (or f ) gives 
y j 


j =-L 

m 

n 


tr [C sc ■] 
e e 


(4-42) 


where 


~ diag 


Y. 


' 'Ce^" 


(4-43) 


The minimization of with respect to the f is indeed a formidable task. Fortunately, 

i 

the results of the optimal control allocation problem (see Appendix B) are applicable. 

The similarity to the control cost functional used to determine the control fractional 

attention f (sampling cost) is readily evident where 
c* 


J = tr [L 2 L '] 
o e e 


(4-44) 


wi 


ith 


Lg =diag (Q^) Le'^'^ 


(4-45) 


Thus, the same techniques used to find 9J ^9 f can be used (with some minor changes) 
to find 9 /b f and to do the subsequent optimizing of subject to the various con- 
straints on f . The computations are greatly simplified in the monitoring case since 
is not affected by f . Thus, using the techniques of Appendix B, 

dj V . 

p.. 


df 


m- 


f 


II 


(4-46) 


m. 

I 
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where 


00 

^ A, 

P = G-y* e^' ''Cg' da . G (4-47) 

0 

G=SC'V ^ = filter gains 

y “ 

* -1 

A -A - SC V C = filter matrix 

y 

The same gradient projection scheme used to minimize Jq is thus directly applicable to 

* 

the moniformg sifuaHon. In order fo avoid a situation where f = 0, a constraint 

* . * 

> e > 0 is imposed to assure each instrument is monitored for failure detection 

purposes, 

4. 2. 3.4 SUMMARY OF MONITORING MODEL 

The monitoring model that was chosen in this study is that of the previous sub- 
section. The monitoring cost functional has various interpretations via the error 
ratios k., and in the case of uncorrelated instruments the model reduces to the Senders* 

I 

model. There is intuitive appeal to the model, plus the ability to solve the optimization 
problem using techniques already developed for solving the optimal control allocation 
problem. 


The one drawback to the model is that the constraint T ^emust be im- 

m. 

I 

posed artificially. As f =0 could not happen in a failure detection oriented metric 

i 

it might seem logical to consider a monitoring cost functional that was a sum of 
estimation and failure detection metrics. Thus, a reasonable modified cost metric 
might be (assuming equal weights to estimation and failure monitoring) 


j =_L 

m 

n 


n CT 2 + p. a 2 /f 

ye. I y- m. 

t 

i=i % 


(4-48) 
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The noise ratio p. is on the order of 0.01; thus for moderate attention levels, minimiz- 
ing Equation (4-48) is equivalent to minimizing Equation (4-7), At low levels of 

attention, however, (say, f ~ 0.01) the second term in Equation (4-48) becomes 

m. 

I 

important and prevents f^ from approaching zero as may happen when minimizing 

2 ' 

. Thus we see that the two following optimization problems will yield nearly 
G* 

I 

identical results. 


min 


n i=i 

y 


a 2 
e. 

I 

a 2 
^i 


+ 


P: 


m. 

1 


(4-49) 



f >e. Ef 

m. I ^ m. m 

I I 


(4-50) 


In conclusion, after an examination of various schemes for predicting human 

monitoring performance, the scheme involving Kalman filter optimization on the full 

aircraft model seems to hold greatest potential. It is not simple mathematically, but 

fortunately the optimization of the f^ can be performed using existing human operator 

i 

computer programs with only slight modification. A default option should be included 

to treat (singular) situations where the optimized filter produces f ~ 0 for some i. 

m* 

I 
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4.3 


SIMULTANEOUS MONITORING AND CONTROL WITHIN A SINGLE AXIS 


The above discussion is relevant to the pilot monitoring assessment for a com- 


pletely automatic system. However, if a pilot is actively controlling a vehicle (or 


axis) there are two reasons for looking at a display variable: 


1 . For control purposes, fraction f . 

^i 

2. For monitoring purposes, fraction f . 

i 

Obviously there is an overlap that must be resolved. The following assumptions are 


made in this situation: 


1 • The pilot allocates his attention first to the control task requirements. 

Spare capacity is then available for monitoring. Thus, f and f are 

c c. 

chosen (subject to ^^TOT ** such that a desired level of ' 
performance is attained. The available monitoring fraction is then 



2. The pilot next allocates f^ among oM displayed variables — 

those being controlled manually and those being controlled automatically — 
to minimize the total monitoring cost. 


J 


m 




(4-51) 


where 


J = (p- / for instruments associated with the manual control 

ml e y 

task . 

J _ = for instruments associated with the automatically 

m2 e' y 

controlled loops. 

3. The optimal f is not dependent on the optimal f (the reverse is not 

^ Cj m. 

generally true). Thus, any additional monitoring of a display does not 
change control strategy . 


For the manually controlled task, it is possible to optimize for f^ and compute 

2 2 ^ 

the relative estimation error fraction o /o^ associated with each displayed variable. 

6 « y . 

I ' I 

Thus, it remains to distribute the f^ by allocating an f^ to each display to minimize 
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Jfn^ + *^012* quite likely that the optimal will be non-zero even for displays 

for which f is a reasonable fraction. This is because the information required of q 

c • 

I 

display for monitoring is generally different than that required for control (i.e., C 


vs L in J vs J ). 
e m o 


Therefore the problem of determining f^ for only a manually controlled loop 


IS solved by first finding f to minimize J . Next find f whileT^f - f to minimize 
' c. o m. m. m 

I III 

J bearing in mind that the total attention to a display for purposes of observation is 

>» 

f,.^T = f + f — i.e., monitoring proceeds over and above control allocation and is 
I vJ I J m. c. 

geared to J . 

^ m 

For the case of an aircraft with both a manual and an automatic control 

loop, the allocation of f is accomplished by considering two "independent" tasks. 

m# 

I 

The following one-dimensional search technique for optimizing and J^ 2 's 

suggested, once the optimal f are found, (Figure 4-1). 

c* 


In some situations it may not be necessary to allocate any additional monitor- 
ing attention to a display used for control purposes already. However It is difficult 
to ascertain this a priori since there may be displays needed only marginally for 
control purposes (e.g., pitch indicator for attitude stability system) but which will most 

likely be monitored to a higher percentage than the control -allocated f . Thus, the 

c* 

I 

above technique makes no a priori assumptions as to whether instruments are used 
primarily for monitoring or for control purposes. 
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Figure 4-1. Flow Diagram for Dual Axis Monitoring Allocations. 
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SECTION 5 
MODEL VALIDATION 


In order to validate the pilot vehicle mpdel, including the means by which 
workload metrics are established, the analysis procedure is applied to a CH-46C VTOL. 
As the CH-46C has been the subject of several studies at NASA Langley Research 

Center (References 12, 17 and 18), there exists a (flight test) data base against which 

1 « ■ * ' 

model predictions may be compared. This analysis will consider the hover task only, as 
this represents one of the most difficult VTOL. control requirements, and is amenable to 
steady-state analysis techniques. The analysis is further restri cted to the display format 
level by considering only the display panel/instruments used in the actual CH-46C 
tests. The Information and/or display element level analysis would be appropriate if 
new display systems were being proposed. The automation levels that we consider are 
the attitude command systems designed dt NASA (Reference 17); the flight directors are 
also the NASA designs (References 12-17). • 

5.1 CH-46C DYNAMIC EQUATIONS OF MOTION 

This section presents the basic equations of motion of the CH-46C, in- 
cluding a representation for the external wind-gusts. It also derives the equations for 
the control stick inputs, taking into account the attitude command system used on the 

CH-46C. 

5.1.1 BASIC VEHICLE DYNAMICS 

The analysis uses the perturbation equations derived in Appendix C, written 
about a hover equilibrium (0 q = 6.8°). The state variables are chosen as the Euler 
angles and body rates to be consistent with the control augmentation and guidance 
schemes employed. Thus, for the longitudinal equations, the state vector is 
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(5-1) 


X = [ax, aZ/ AVj^/ AV^/'q]' 
where ' 

Ax = longitudinal (inertial) hover error 
A z= vertical (inertial) hover error 
A9= pitch deviations (from 9 q) 

= longitudinal (inertial) velocity 
AV^ = vertical (inertial) velocity 
q = body-axis pitch rate 

The control inputs for the basic vehicle are the differential and gang collective 

I 

deviations from trim, so that u = [A 6 g, A 6 g] . The longitudinal equations, in the form 
X = Ax + Bu (5-2) 

are given in Figure 5-1. 

For the. lateral direction dynamics we choose as state vector 

x= (Ay/ 0 , AY, AVy, p, r] , (5-3) 

where 

Ay = lateral (inertial) hovering error 
0 = inertial roll angle r ' 

A y = heading angle error ■ 

AVy= lateral (inertial) velocity 
p = body-axis roll rate 
r = body-axis yaw rate 

These states have been chosen to be compatible with subsequent model following aug- 
mentation schemes. The two lateral controls are u = [A6^, AA^,] where A 6^ = roll 
cyclic input and A 6^ = rudder or yaw cyclic input. The numerical values for the A and 
B matrices are included in Figure 5-1 . 
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Figure 5-1. System Matrices for Basic CH~46 Helicopter at Hover, 



5.1.2 AUGMENTATION SCHEMES 


The basic CH-46 is unstable and virtually Impossible for a pilot to control 

without some form of stability augmentation system or command following system. The 
automation schemes that will be investigated are the NASA designed pitch, roll and 
heading command systems using the differential collective, roll cyclic and rudder, 
respectively. The pitch command model is 



with 0.75, (U= 2.0 rad/sec. Kg - 0.15 rad/in. The roll model 


(5-4) 


P„ = 2 c:uj.p_ + 6 - . 

m ^ m 0 a 


(5.5) 


is identical to the pitch model with K = K • The heading hold model in yaw, 

(J 0 




(5-6) 


has Cq “ 0*7, w = 2.0 rad/sec and K^” 0.35 rad/in. 


Having formulated the desired model characteristics, the next step is to 
develop the compensation that will force the aircraft response to follow the model. In 
the NASA-LoRC design, the compensation used included lead, rate error and integrated 
rate error. The resulting frequency response characteristics showed that the aircraft 
response closely followed attitude commands over a wide frequency range, beyond that 
normally associated with pilot control bandwidth. If one wishes to precisely model 
pilot-vehicle response for the augmented systems it is necessary to include the model- 
following compensation in the system description. This adds additional states in the 
vehicle equations, needlessly complicating the analysis process since the pilot effec- 
tively "sees" a system that responds as the model . The simplified approach we follow 
in studying pilot response is to assume that the system exactly follows the model, and 
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solve for the feedbacks that are required to assure this condition. Additional vehifcle 
states are not needed, but a difficulty arises because of cross-coupling between the con- 
trol actuator dedicated to provide the model response and the effect of that actuator in 
other state variables. For example, a differential collective may be used to satisfy the 
pitch response of the helicopter, but by so doing, it wi|l have an influence in the 
translational equations. 


To show the means by which the augmented equations are developed, assume 
that the basic unaugmented vehicle equations can be written (Appendix D) 



(5-7) 


We hove divided the state vector into two parts with X 2 that portion of the state vector 
to be controlled according to some desired model response. The inputs to this system are 
the actuator inputs & which are dedicated to satisfying the model response. Disturbance 
inputs and other actuators affecting the equations of motion, but not being used to 


satisfy the mode! response can be added to Equation (5-7) but are riot shown here. 


Assume that it is desired to have the partition of the state vector X 2 


follow 


a model equation given by 


2m ^21^1 


+ A*„ 


22 ''2m 


+ B*u 


(5-8) 


Note that the derivative of the model response is given by linear combinations of the 
model itself ^ 2 ^, feedback of the other partition of the state vector x^ (e.g. position), 
and the control stick inputs u. 

The actual response of X 2 is given in Equation (5-:7) as 

X2 ^21 ^22 ^2 ^ 2 ^ (^*^) 
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In order that the system response given by Equation (5-9) follow the model response 
of Equation (5-8), i.e,, X 2 ~ subtract the two equations and solve for the 

actuator activity which is dedicated to satisfy this responses 

& = 82^ t(A|^ - A2P + (A^2 “ A22) ^2 (5-10) 


Substituting this value of the actuator activity into Equation (5-7) provides the following 
form for the state equations when the system follows the model. 


where 



A * i 2 “ A ^2 ^^22 " ^ 22 ^ 


(5-11) 


. V 

(5-12)' 


Note that the partition of the state vector X2 has the desired response of the model, but 
also that the dynamics of- the other partition of the state vector x^ have been altered be- 
cause of the effect of the actuator 6 on these equations. 


The above mode I -foil owing process is easily applied to include the command 
systems of Equations (5-4) through (5-6). For example, in the pitch command model 
(Equation (5-4)), the differential collective is the dedicated actuator, and the par- 
titions A2^, A22/ ^2 are shown in Figure 5-1; the model partitions are 

A*, = [0 0 -0)^ 0 01 

A^2 = 

B. = 
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The resulting A*, B* matrices for the augmented system are shown in Figure 5-2. Note 
that the second column of the B matrix, associated with the power collective, is un- 
modified since this control is not used in the model matching scheme. 

In the lateral case, both controls are used in the model matching process. 
The partitions ^ 2 ^, ^ 22 / B 2 are shown in Figure 5-1; the model partitions follow from 
Equations (5-5) and (5-6) as 



and provide the lateral directional equations of motion for the aircraft with the roll and 
yaw augmentation systems engaged. 

5.1.3 WIND GUST DISTURBANCES 

The principal external disturbances acting on the helicopter are due to the 
wind turbulence or gusts. The representation of these effects in the linear system model 
involves two considerations: a model of the wind spectrum, and the way in which the 
wind affects the vehicle aerodynamic forces. The gust model that is used is the Dryden 
model with power spectral density 






(5-13) 


To obtain a random wind model with this PSD, a white noise T] (t) with variance 
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Longitudinal Direction; 



cn 

I , ‘ 

°° Lateral Direction: 



Figure 5-2. System Matrices for Augmented CH-46 Helicopter at Hover. 



o V 

E [Tl(OTi(t-T)l = -2. 

9 L 


(5-14) 


is passed through the linear system 


V 

5(t) = -^ ?(t) + Tl(t) 

L 


(5-15) 


In our modeling work we will assume moderate turbulence with representative parameters 
a ■= rms gust velocity = 3 ft/sec 

■g ' 

V^L = gust bandwidth =?= 0.1 rad/sec 

We assume that independent gusts pet along all three of the aircraft (body) axes. A 
normal component induces primarily normal accelerations, or perturbations in angle- 
of-attack. An, axial component induces primarily forward speed perturbations, and a 
lateral component 5^ perturbs lateral velocity. For simplicity we assume ~ 

= 3 ft/sec. 

The components, 5 ^(t), Cy(t) and ^(t) are subtracted from the dynamical 
states w(t), u(t) and v(t), respectively, in the equations of motion. The result is. the • 
modifications 


u = - 


w 


m 


w 


-2- I + . 

m u 


v = I + 

m V 


w = 


p = - 


Z 

— ? 
m w 


— 5 , + 

XX 


Z 

- 2 - ? + 
m u 




q = 


M 

T 


M 

. " ?w - ^u 

yy yy 

r ---f— % + ... 
zz 


(5-16) 
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When these terms are carried through the. Euler transformations to the inertial frame state 
variables as defined in Equations (5-l)and (5-3), the end result is the addition of the terms 

''w W (5-I7a) 

to the longitudinal equations, and the term 

' (5-17b) 

to the lateral equationst. For the CH-4^C at hover, these terms are 


" 0 “ 


“ 0 - 


“ 0 “ 

0 ■ 


0 ' i 


0 

0 

0.01535 

F = 

. . , .'W 

0 

0.00595 

F = 

; , . V 

0 • 
0.02663 

-0.0634 


0.3733 


0.00789 

-0.00656 


0.00285 


-0.00097 


Note that these terms are appropriate to both the basic and augmented vehicles as the 
wind gust'states are not included in model following feedbacks. Thus, we arrive at the 
final model for the system/gust dynamics compatible with the pilot model formulation Of 
Appendix A, 

i(t) = Ax(t) + Bu(t) + EC(t) (5-18) 

where, in the longitudinal case 

x = [lu/ V ^x' (^-’9a) 

and, in the lateral case 

X = '^'7 Vy/ P, r] (5-19b) 

and where |(t) is a white noise vector. 

tAdditional states must be added to the system dynamics to include the noise shaping 
dynamics Equation (5-15) for each §. . 
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5.2 CH-46 HOVER DISPLAY/CONTROL PARAMETERS 

In applying the optimal control model of the pilot, it is necessary to describe 
the display informatipn in a form compatible with vehicle states and controls, 

y(t) = Cx(t) + Du(t) (5-20) 

In addition, values must be chosen for the visual and/or indifference thresholds, i.e. 
those limits on displayed items within which the pilot is less likely to take corrective 
action. The values are estimates derived from the display gains> display resolution, 
display markings and pilot opinions. 

5.2.1 PRIMARY DISPLAY INPUTS - HOVER MODE 

The four primary instruments gppropriate to the hover mode- are the map (or 
HSI), radar altimeter, IVSI and ADI. Assuming that the pil,ot perceives both the posi- 
tion and rate of a display indicator, we have the model inputs listed in Table 5-V. 
Alsog^iven in Table 5-1 are the values chosen for the display thresholds, as well as sub- 
jective values for the maximum deviations that a pilot will tolerate in the observed 
quantities. These latter values are used tq generate cost functipnal weightings in the 
optimal control model according to 

. 1 


1 \y\, max J 


(5-21) 


y; \ / ; max 

The large thresholds associated with perceiving horizontal position and veloci* 

» * s , 

ty errors at hover is a reflection of the small map scale (100 ft/in at hover) and the 
relatively large size of the aircraft symbol^. The associated maximum deviations for 
^x. Ay reflect the 50 ft diameter of the landing pad. The small threshold on altitude 
deviations reflects the expanded scale of the radar altimeter below 100 ft, with an 


tAt a viewing distance of 28 in, a 20 ft error is less than 0.5° visual arc. 
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Table 5-1 . CH-46C Primary Display Inputs (Hover Mode). 


INSTRUMENT 

VARIABLE 

THRESHOLD 

MAX. DEVIATION 

MAP 

AX 

20 ft 

25 ft 


Vx 

10 ft/sec 


ALTIM 

tx 

4 ft 

5 ft 


■ V, , . . . 

2 ft/sec ■ 

2 ft/sec 

ADI 

9 

1° 

* ; 

2° 


q 

0.5°/sec 

l°/sec 

IVSI 

Vz 

1 ft/sec 

2 ft/sec 


Vz 

0.5 ft/sec 

— 

MAP 

Ay 

20 ft 

25 ft 

• 

Vy 

10 ft/sec 

— 

ADI 


lo 

5° 


, . _ i . . ; , ; 

0.5°/sec 

2°/sec , 

MAP/COMPASS 


lo 

20 


■ y- ^ 

0.5°/sec 

l°/sec 


instrument scale marking at the nominal 50ft hovering altitude. Finally, we have chosen 
thresholds on the rates of a display indicator equal to 1/2 the corresponding position 
threshold in accordance with previous work (Reference 21), 

5.2.2 FLIGHT DIRECTOR INPUTS - HOVER MODE 

The flight director has undergone many revisions througljout the course of the 
flight test of the CH-46C. Herein, the director form and values for the flight director 
gains at hover are taken from Reference 12; Table 1 1 .There are three flight director 
signals corresponding to pitch, power and roll commands respectively. The form of 
these signals is. (neglecting the effects 6f any additional filtering). 
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(5-22) 


, FD^= .004 AX + .04 V + .56 

FD^ = .0057 Az + .057 + . J 15 V^ 

FD = .004 Ay + .04 V + .56 
y y a 

In order fo express the signal FD^in the requisite form of Equation (5-20), we substitute 
for V^ using the equations of motion. The result is (approximately) 

FD^- .04 + .0057 Az + .04 9 + .0515 V^ - .068q - .85 6^ (5-23) 

or, if we keep only the dominant terms, 

FD^ .0057 AZ + .015 V^ t .85 6^ (5-24) 

The rates of the flight director signals were not included as "observed" 
variables in our application of the pilot model. There were several reasons for this 
decision: 

• The director rates would contain control rate (e) and white 
noise terms, giving=rise to a complex modeling problem 
that could only be solved by adding "pseudo" filters in the 

. equations. 

• The flight director signals are basically high frequency in 
nature so that the value of their rate information in a low 
bandwidth hover task is small. 

• The flight directors have been designed so that the pilot 
acts essentially as a gain (especially at low frequencies) 
on the observed signal . This minimizes the pilot reliance 
on the rate information. 

In order to include the flight directors within the modeling framework we require a 
selection of threshold values and cost functional weightings. The thresholds were 
chosen as 0.1 in for all three directors, corresponding to approximately 1/10 full scale. 
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As discussed in Subsection 3.2, including the FD. within the cost functional, in addition 
to the other terms, implies that one of the pilot's direct control objectives is to keep 
the FDj small. Based on an analysis of the flight director instruments, and piloting 
task, we choose as maximum deviations ’ 


5.3 


K)m AX 

MODEL PREDICTIONS FOR CH-46 HOVER 


Having obtained state-space equations for the vehicle arid display dynamics, 
the ability of the pilot model to analyze the relationship between CH-46 hovering per- 
formance and pilot workload is investigated. We consider pilot control performance for 
the augmented (attitude command) system, both with and without use of the flight di- 
rector signals. Various levels of total attention f from 0.3 to 0.8 are assumed so that we 
may obtain the performance vs. workload curves. The cost funbtional output weight- 
ings, q , are as given in the previous section. The control rate' weightings q* are ad- 
Y\ ^i 

justed to give resulting "neuroi-motor" lags of 0.1 sec, with a minimum value of 

q, = 0.25 allowable. This corresponds to an assumed maximum pilot stick-rate deflec- 
^i 

tion of 2 in/sec. The remaining basic parameters for the pilot model are selected at 

their nominal d priori values: 

> =0.2 sec 




= 0.01 (i.e., -20 dB) for all displayed variables 


~ 0.01 (-20 dB) for all control variables 
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5.3.1 MODEL PREDICTIONS - NO FLIGHT DIRECTOR CASE 

The pilot model is used to predict hovering performance with only primary 
instruments. In the flight tests, the pilots were unable to maintain the CH-46 in a 
hover under these conditions. Table 5-2 gives model results corresponding ,to an allo- 
cation of attention of 0.3 to the longitudinal axis. The results clearly show that accu- 
rate hover is virtually impossible Jongitudinally. One-sigma errors of 26 ft imply that 
about 40 percent of the time the helicopter is not over the pad. Instrument ottentional allo- 
cations are equally divided between map and IVSI, with some attention to the altimetert 
Thus, height information is obtained primarily through integrating the IVSI with cross- 
checks from the altimeter. 

Table 5-2. , Model Results. Longitudinal Axis, f. = 0.3. 


Instrument 

^i 

'c. 

1 

f /f 
c. c 
1 

MAP, X 

25.5 ft 

0.11 

0 .37 

ALTIM, z 

4.6 ft 

0.05 

0 .17 

ADI, 9 

1° 

< 0.01 

<0.03 

IVSI, V^ 

1 .3 ft/sec 

0.12 

0 .40 


The following information relative to the longitudinal task has also been 
obtained from the model analysis: 

• The total cost J(u) = 3.01 whereas the scanning cost (the 
part of J due to human observational processing noise) l(u) = 

2.86. Thus, over 95 percent of the hovering errors are 

due to the human's own injected random errors> and only 
5 percent are directly attributable to the wind gusts. 

• The sensitivity of performance to higher ottentional 
(workload) levels is fairly low. Thus, f^ = 0.5 re- 
sults in only 12. percent improvement. in hovering errors. 
Consequently, even very high workload levels will 

■not result in satisfactory hover. 

■fThe given f^ are accurate to +10 percent as in the vicinity of the optimal. The cost 
functional J(uj is relatively flat with respect to f . 
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• Rate inforrtiation obtained from the displays (i.e., the 
derivatives of displayed quantities) is not very useful 
for control. Running the model with only the four posi- 
tion variables ( 9 , Vj.) gave attentional and 

performance results within 5 jjercent of Table 5-2. 

The results for the lateral axis hover task, corresponding to an attention 
level of 0.3, are given in Table 5-3. As can be seen, hover performance is unaccept- 
able. However, the attitude hold systems are p'erfdrming ds expected . ' As d result, 
little pilot attention is placed on the roll and'cdurse indicators. 


Table 5-3. Model Results. Lateral Axis, f =0.3. 

c 


Instrument 

a 

1 

f 

c. 

1 

f A 

c/ c 
1 

MAP, y 

19.4 

0.27 

0.90 

ADI, 0 


0 - 02 ,. 

0.06 

DG/t 

0 

• 

0 

< 0.01 

0.04 

'a 

0.17 in. 

- 


6 

0.05 in. 

_ 


r 





The results correspond with the pilot's inability to hover on the primary 
status instruments. Although the model predicts a hover rms of 25 feet, this result is 


under the assumption of a well -trained pilot. Thus, an experiment where only several 
runs are tried is likely to result in much larger errors. Finally, the model predicted 



observed experimentally. 


5.3.2 SINGLE AXIS MODEL RESULTS - WITH FLIGHT DIRECTOR 

The results using the, .flight director in the, longitudinal task are given in 
Table 5-4 for a 0.3 level of attention. There is a substantiql improvement in rms 
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hovering error from 25.5 to 17 feet, accompanied by only slight changes in the perform- 
ance of other variables. What is interesting, of course, is the model's prediction of a 
dramatic shift of attention away from the status instruments onto the flight directors. 


Table 5-4. Model Predictions Using Flight Director, f = 0.3. 

c 


Instrument 

a. 

1 . 

fc 

°i 

fcjAc 

MAP, X 

16.7 ft 

1 

<0.01 

<0.03 

ALTIM, z 

. 4.7 ft 

0.025 

0.1 

ADI, 9 

0.7° 

<0.01 

<0.03 

IV5I, 

1 . 2 ft/sec 

0.025 

0.1 

FD 

X 

0.13 

0.10 

0.3 

FD 

z 

0.14 

0.13 

0.45 


A sensitivity aiia lysis was performed for the longitudinal axis by running the model at 
assumed workload levels of 0.3 to 0.7. The optimal attentional allocations to the 
flight directors are evident in all cases, as seen in Table 5-5. 


Table 5-5. Attentional Allocations for Longitudinal Task. 


Total 

Map 

Al tin 

ADI 

IVSI 

FD 

X 

FD 

z 

. 0.3 

<0.01 

0.025 

<0.01 

-0.02 

0.10 

0.13 

0.4 

<0.01 

0.027 

<0.01 

-0.02 

0.12 

0.21 

0.5 

<0.01 

0.037 

<0.01 

.-0.02 

0.17 

0.25 

0.6 

<0.01 

0.043 

<0.01 

-0.02 

0.22 

0.30 

0.7 

<0.01 

0.045 

<0.01 

-0.02 

0.27 

0.36 
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Figure 5-3 shows the sensitivity of the perfonrnance measures (both scanning cost and x 
hover error) to workload level on the longitudinal task. The results show that a level of 
~ 0,6 is required to keep the aircraft oyer the pad in the x-direction 95 percent: of the time . 

The results for the lateral axis are given in Table 5-6. RMS hovering error 
improves to 13 feet at 0.3 attention level with a total shift of attention to the flight 
director. A sensitivity analysis of lateral axis performance is shown in Figure 5-4. For 
attention levels greater than 0,15, scanning cost is reasonable insensitive. For each at- 
tention level, the allocations among instruments remain the same - full attention to the 
flight director! 

Table 5-6. Model Results. Lateral Axis with Flight Director, f^ =0.3. 


Instrument 

■ a. 
1 

f. , 
1 

f /f 
^i ^ 

MAP, y 

13.0 

<0.01 

< 0.03 

ADI, 0 

1.1°^ : 

<0.01 

<0.03 

6 

a 

0.4^ 

<0.01 

< 0i03 ' • 

FD 

’ 0 . 1 in 

0.27-1- 

0.9+ 

y 





5.3.3 DUAL AXIS MODEL RESULTS - WITH FLIGHT DIRECTOR 

The curves of lateral and longitudinal cost, l(u), vs attentional allocation 
can be used to determine performance for the combined task given a level of total con- 
trol workload f^ = ^ Thus for a given f^ the optimal and f^^ 

are determined to minimize 

'(“'total = + '(“’long (5-25) 

This minimization is accomplished here using a simple one dimensional search procedure. 
It could also be done via computer by constructing a state space representation for the 
entire two axis task. In this example, the first approach is easier. 
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Figure 5-3, Scoresvs Attention, CH-^46 Hover, Longitudinal Axis. 
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RMS Hovering Error, 



Cost l(u) 



0 0.1 0.2 0.3 


Attention to Lateral Task 

Figure 5-4. Scores vs Attention, CH-46 Hover, Lateral Axis. > 
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The results are shown in Figure 5-5 which represents the overall predictions 
of the model for the hover task. The following comments are pertinent: 

• The optimal attention to the lateral task fc*Ipt Q-15 -' 0.2 
or in the vicinity of the "knee" in Figure 5^4 as might Be 
expected. 

• At the optimal point (f..*u* / f *i ) the RMS hover errors 

a ~ a towithin5perclnt. ; 

^ y 

• The sensfitivity of tQtol l(u), in the vicinity of the optimum 

is extremely low. Thus lateral attention could be increased by 
0.05, while longitudinal attention is decreased by 0.05, with 
little change in l(u), . ' ' 


The "overage" hovering error 


2 4 . „ 2 
a + C7 

X y 


V 1/2 


ovg 


(5-26) 


shows the same relative insensitivity as does l(u). How- 

ever> there is a fairly sharp trade-bff between ^ and (j . 

X y 

• Thus, it is quite likely that different pilots will adopt dif- 
ferent operating points with respect to lateral vs. longitudinal , 
trade-offs. Total performance J(u) is insensitive here but dif- 
ferent distribution of errors is to be expected. 

It is possible to investigate the probability of hovering within a circle of 
radius R for different workload levels. Since we assume Gaussian statistics, the proba- 
bility densities of x and y hover errors are, respectively. 


2/„ 2 

1 /^^x 

P(x) = b 


(5-27) 


2/0 2 

1 -y /2 c 7 
p(y) = = e 


y 


a \/2n 

y 


(5-28) 


9 99 •• 22 2 . 

Assuming ~ cr'^ = ct^, it is easy to show that the probability x + y < R is 

X y 


Pr + y^ < R 3 ~ 1 - e 




(5-29) 
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Figure 5-5. Predicted Performance, CH -46 Hover, Flight Director Used. 
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Thus, the probability of being within a 25 ft radius at any time is easily found for any 
given "overage" hovering error. Several probabilities are noted in Figure 5-5. 

We see that workload levels of 0.7 - 0,8 are required to maintain the aircraft 
over the pad 80 to 85 percent of the time.. Thus at this criterion level on performance, we - 
find that the workload level Ts barely acceptable,^ The hover task can be accorhplished,but 

no time is left for monitoring the status instruments. (Recall most attention goes to the 
flight-directors). If we require a tighter hovering performance, the model predicts that 
the pilot will be overworked. At lower workload levels (e.g., 0,6) the aircraft is over 
the pad only 75 percent of the time. 

The conclusions of our modeling effort thus indicate that 

• Hovering cannot be satisfactorily accomplished us3ng the 
existing status display. 

• The flight directors make hovering possible with high 
workload 

• Virtually full control attention is given to the flight 
director needles 

• There is little or no remaining pilot capacity to monitor the 
status instruments. 

These conclusions are very much the same as those from the flight experiments. 


t 


We consider fj. 


0.8 to be full capacity. 
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SECTION 6 


, . CH-47 RESULTS 

The synthesis and analysis techniques described in the preceding section were 
applied to the Langley CH-47 tandem rotor helicopter, which will be used as the flight 
research aircraft for the VALT prpgram. The numerical results presented in this section 
are limited to the longitudinal axes due to time constraints and uncertainty in the 
lateral directionaf stability derivative data. Two flight conditions were analyzed to 
indicate thq differences in performance over the expected flight environment:' 

• hover (sea level) 

' • approach (l000,ft/min descent at 60 knots) 


6. 1 ^ CONTROL AUTOMATI ON LEVELS 


As discussed previously in Section 4, each of the control channels of the 
helicopter has three or five possible levels of automation, shown in Table 6-1 . 


Table 6-1. Levels of Control Channel Automation. 


Forward 

Vertical 

Lateral 

Directional 

6 


e 

6 

e 

c 

a 
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The total number of automation levels from Table 6-1 is 5x3x5x3 = 225. Obviously it 
would be impractical to analyze all combinations of automation levels, and indeed 
many of the combinations in Table 6-1 would not be practical systems. Consequently, 
several levels of control augmentation were selected to perform the control/display 
analyses. Table 6-2 shows the control channel commands for each of. the eight selected 
systems. These range from the completely manual, basic vehicle (System A) to a full 
position feedback system for automatic hover (System H). Note that for the longitudinal 
axes, Systerns D and E are identical. 


Table 6-2. Selected CH-47 Automation Levels.. 


System 

Control Channel Comrriand 

Forward 

Vertical 

Lateral 

Directional 

- A 

■ • 6 

6 V 


-6 - 


e 

c 

a 

r 

B 

' q 

' ■■'6c 

p ' 

r ' 

C 


■'■‘V 

0- . 

r 

■ D 

9 

V : 

. 0 ■ 

r 



z 


i 

E 

e 

V 

z 

0 

t ' 

F 

6 

h : 

0 

' ^ 

G 

V 

h 

V 

■ I' 


X 


y 


H 

X 

h 

y 

i '1' ■ 


The quadratic synthesis technique described in Section 4 was used to obtain 
the closed loop longitudinal system matrices for each of the systems in Table 6-2. As 
discussed in Section 4, the weighting functions used in the quadratic synthesis can all 
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be related to the specifications on the maximum allowable pitch acceleration and 
vertical acceleration. For the CH-47, the maximum state and control variable excur- 
sions which define the weighting functions are shown in Table 6-3. 


Table 6-3.. Weighting Functions for Longitudinal 
CH-47 Automation Levels. 


State Variables - Hover and Approach 

System 

X 

(fo 

z 

(fo 

e 

(rad) 

. Vx , 

(ft/sec) 

V, 

(ft/sec) 

(rad/sec) 

A 

B 

C 

P>E 

G 

- > . 

84,0 

' 30.0. 

30.0 

30.0 

0.435 

0.435' 

0v435 

28.0 

7.5 

0.435 



Control Variables: 

Systems A-H 



FI ight Condition 

6 3 (in) 

6 

(in) 
c ' ' 


Hover 




2.19 


0.336 


Approach 



2.45 


0.404 
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Recall that the dynamics of the system to be controlled, by the pilot are 


given by 

X = Ax + Bu + Ew 


( 6 - 1 ) 


where the state vector x consists of the vehicle states shown in Table 6-3 augmented by 
shaping states for the disturbances. From the CH-46 results it was concluded that a 
first-order wind gust in both the forward and vertical directions provided a sufficient 
disturbance input to the system. Consequently, the system state x in Equation (6-1) is 
augmented by two wind gust velocities given by 


w = -w /t + w 
gx gx' X X 


w = -w /t + w 
gz gz' 2 z 


(6-2) 


Again from the CH-46 results the time constants and were each selected qs ten 
seconds, and the power of the white driving noises (w^, w^) were chosen to gjve qn 

rms gust velocity of 3 ft/sec. 


The quadratic synthesis technique was applied for the above weighting func- 
tions at hover and approach to generate the closed loop system matrices for Systems B 
through H (the wind disturbances were not included in this process). The resulting ' 
closed loop dynamics for each system are shown in Figures 6-1 through 6-7. 


6.2 DISPLAY LEVELS 

The optimum control pilot model was used to determine the attention alloca- 
tions and system performance at the information level and at the display element level 
for each of the systems described previously. Without flight director signals, the 
control performance metric weights quadratic terms for each of the basic aircraft states 
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System Dynamics Matrices, System A 
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Figure 6-2. System Dynamics Matrices, System B. 
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Figure 6-3. Sysfem Dynamics Matrices, System C. 
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Figure 6-4. System Dynamics Malrices, System D, 
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Figure 6-5. System Dynamics Matrices, System F 
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Figure 6-7. System Dynamics 'Matrices, System H 



(excluding gust states) and the pilot control input rates, as discussed in Section 2. 
Recall that the displayed information vector presented to the pilot Is 


y = Cx + Du (6”3) 

For the information level and the display element level v^ithout flight director, the 
matrices C and D are shov/n in Figure 6-8. Desirable maximum allowable values (2 a) 
are presented in Table 6-4 for each of the basic vehicle states and’ for the control rates. 
Table 6-4 also shows the indifference thresholds for each of the displayed states. As 
discussed previously, these are taken to be 1/4 of their respective maximum allowable 
values in the cost functional. 

Flight director algorithms were obtained for each automation level/flight 
condition as described in Section 4. The flight director commands were simplified by 
neglecting cross coupling terms to give the following general forms: 




FD 

z 


K^Z + Ky 
z 


V 

z 


Table 6-4. Control Cost Weighting Functions and Indifference 
Thresholds for CH-47. 


(6-4) 


Variable 

- 

Units 

Maximum Values for 
Cost Function 

Indifference Thresholds 

Hover 

60 kt Descent 

Hover 

60 kt Descent 

X 

ft 

5. 

25. 

1 .25 

6.25 

z 

ft 

5. 

25. 

1 .25 

6.25 

e 

deg 

I. 

I . 

.25 

.25 

V 

X 

ft/sec 

I. 

2.5 

.25 

.625 

V 

z 

ft/sec 

?. 

2.5 

.25 

.625 

q 

deg/sec 

0.5 

0.5 

.125 

.125 
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e 

in/sec 

2. 

2. 
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In/sec 

2. 

2. 
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Figure 6-8. Display Vector for Basic Information. 
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The flight director signal rates were obtained by simply differentiating Equation (6-4) 
and substituting from Equation (6-1). The effect of the flight director signal presence 
on the pilot strategy is taken into accounf by including the flight director signals in 
the control cost functional . The maximum allowable values used to weight these sig- 
nals in the cost functional are obtained from the maximum allowable excursions of the 
position states, e.g., 

= l'<xMmaxlHKg(e)„„J 

(6-5) 

(FD ) = K (z) 

' z^max ' z' 'max ' 

Finally, fhe indifference thresholds for the flight director signals are obtained in the 
standard 1/4 of the associated maximum value, i.e. 

(FD ) . = (FD ) /4 

(6-6) 

(FD ) . = (FD ) /4 

' z'min ' z'^max'^ 

The resulting flight director algorithms, indifference thresholds, and control cost 
v/eightings are shown in Figures 6-9 through 6-15 for the seven automation levels. 

6.3 PROGRAM PIREP 

Program PjREP is an interactive program that predicts human operator per- 
formance and response characteristics in complex control tasks involving varying levels 
of control automation and pilot monitoring. The human operator model is based on 
optimal control and estimation theory coupled with a mathematical description of the 
human*s limitations. The basic underlying assumption is that the well -motivated, well- 
trained human operator behaves in a near optimal manner, subject to his inherent 
limitations and constraints and his control task. The program subroutines are written 
in Fortran IV and have been operated on PDP-10, IBM 360-65 and CDC 6600 computers. 
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Figure 6-9. Flight Director Algorithms, System A. 
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Figure 6-10. Flight Director Algorithms, System B. 
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Figure 6-1 1 . Flight Director Algorithms, System C. 
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Figure 6-12. Flight Director Algorithms, System D. 
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Figure 6-13. Flight Director Algorithms/ System F. 
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Figure 6-14. Flight Director Algorithms, System G. 
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Figure 6-15. Flight Direcfor Algoririims, System H. 








The inputs to the program include the vehicle description (system dynamics, 
input disturbances and displayed variables), task description (cost functional weight- 
ings), human description (subjective weightings, neuromotor time constants, time 
delay, motor noise and observation noise), and the available capacity for control an d/ 
or monitoring tasks. The outputs include the state, output and control variances; the 
optimal allocation of pilot attention to fhe prescribed displays; control and monitoring 
performance metrics; pilot describing function and remnant spectra; and the state, 
output and control power spectra. 


Figure 6-16 summarizes the important parameters of program PI REP for this 
analysis. The seven control systems, four displays, two flight conditions and varying 
control attention levels required over 200 runs of program PIREP. The following two 
subsections present selected results obtained with program PIREP for the control and 
monitoring models. 


System Dynamics 
(Automation Level, 
Disturbances, etc.) 

Display Model 
(Elements, 
Thresholds, etc.) 

Total Attention 
(Control/ 
Monitoring) 


PIREP 


Attention Allocation 
(Control/ 
Monitoring) 

Performance Metrics 

(J , J , Gradients) 
'em' ' 


RMS Predictions 
(State and Control) 


Figure 6-16. Principal Inputs and Outputs of Program PIREP. 


6.4 CONTROL MODEL RESULTS 


This subsection contains comprehensive results obtained using program 
PIREP to analyze the control performance of the aforementioned control/display 
configurations selected for the CH-47. For the reader's convenience, it is suggested 
that Table 6-2 on page 6-2 be marked for convenient reference during the ensuing 
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discussion. The pilot's total allowable attention to the longitudinal task was assumed 
to be 0.6, leaving at least 40 percent of his time for control/monitoring of the lateral 
axes and for miscellaneous duties (communications, etc.). 

6.4.1 ATTENTION ALLOCATION VS WORKLOAD 

Figure 6-17 depicts the relative attention (in percent of total control 

attention f ) paid by the pilot to the display elements for the unaugmented helicopter 
c 

(System A) with only situation information. At both flight conditions, most of his 
attention is on the pitch attitude indicator. From hover to approach, some of his 
attention is shifted from attitude to the position displays. For this simple control/ 
display configuration, his relative allocation of attention remains essentially constant. 
For comparison. Figure 6-18 shows the relative attention allocation for System A with 
a full flight director. Now the forward flight director is the most Important instrument, 
followed by the pitch indicator; the vertical flight director is significant only during 
hover. For this configuration, the pilot spends a lower percentage of his time on the 
X flight director, and more on the attitude indicator, as his total available control 
increases. Figure 6-19 shows the predicted rms system performance for the unaugmented 
system without a flight director. In general, the performance improves as the pilot 
spends more time on the longitudinal control task. The hover errors are consistently 
lower than those during approach, reflecting the cost weightings of Table 6-4. 

Figures 6-20 through 6-31 compare the fractional control attention alloca- 
tions for the other six automation levels with ho flight director and with full flight 
director. In general, as the control system automation increases, the pilot is able to 
direct more attention to the vertical information and thereby improve the total perform- 
ance cost (e.g. Figures 6-28 through 6-31). Also, as his total control attention 
increases, he tends to spend slightly less time on the flight director signals and more on 
the raw data. 
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Figure 6-17. Attention Allocation vs Total Control Attention 
(System A - No Flight Director). , 
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Figure 6-18. Attention Allocation vs Total Control Attention 
(System A - Full Flight Director). 
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Figure 6-20. Attention Allocation vs Total Control Attention 
(System B - No Flight Director). 
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Figure 6-21 . Attention Allocation vs Total Control Attention 
(System B - Full Flight Director). 
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Figure 6-22. Attention Allocation vs Total Control Attention 
(System C - No Flight Director). 
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Figure 6-23. Attention Allocation vs Total Control Attention 
(System C - Full Flight Director). 
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Figure 6-24. Attention Allocation vs Total Control Attention 
(Systems D, E - No Flight Director). 
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Figure 6-25. Attention Allocation vs Total Control Attention 
(Systems D, E - Full Flight Director). 
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Figure 6-26. Attention Allocation vs Total Control Attention 
(System F - No Flight Director). 
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Figure 6 - 27 -. Attention Allocation, vs Total Control Attention 
(System F - Full Flight Director). 
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Figure 6-28. Attention Allocation vs Total Control Attention 
(System G - No Flight Director). 
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Figure 6-29. Attention Allocation vs Total Control Attention 
(System G - Full Flight Director). 
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Figure 6-30. Attention Allocation vs Total Control Attention 
(System H - No Flight Director). 
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Figure 6-31. Attention Allocation vs Total Control Attention 
(System H - Full Flight Director). . 
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6.4.2 


SYSTEM PERFORMANCE VS WORKLOAD 


Figures 6-32 through 6-35 present the overall performance cost variation 
with workload (e.g. control attention) at hover for each control/display configuration. 
Increased automation improves performance for a given workload, regardless of the 
display provided. Comparison of these figures shows that the vertical flight director 
alone helps only system G's performance, while the forward flight director improves 
all systems' performance. However, the full flight director gives only a marginal 
improvement over the x flight director by itself. 

Similar results are obtained during the approach, as shown in Figures 6-36 
through 6-39, In nearly all configurations the overall performance is better during the 

approach than at hover. However, the performance of System H (full position feedback) 
is approximately the same for both flight conditions. Moreover System H shows very 
little sensitivity either to the flight director level or to the total control attention at 
both flight conditions. This result is as expected, since System H is essentially a pure 
automatic system. In general, the more automated systems become considerably less 
sensitive to workload variations.^ 

6.4.3 PREDICTED RMS PERFORMANCE 

Figures 6-40 and 6-41 show the rms position and attitude error contributions 

to the overall system performance, at a "comfortable" workload of 0.4. The position 

errors are considerably lower at hover than during approach, while the attitude errors 

are nearly the same; this reflects the higher penalties on x,z in the performance metric 

at hover (Table 6-2). In general (exceptSystem D), the errors tend to decrease with 

system automation and the flight director also reduces the errors. 

^An exception of these generalizations is System D, whose predicted performance is 
extremely good during approach and very poor at hover. The erratic results are not 
yet fully understood, but it appears that numerical difficulties with this system are 
producing two local minima in the optimization algorithm. 
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Figure 6-33. Control Performonce Cost vs . Total Control Attention 
(Hover - z Flight Director), 
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Figure 6-35, Control Performance Cost vs Total Control Attention 
(Hover - Full Flight Director). 
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Figure 6-36. Control Performance Cost vs Total Control Attention 
(Approach - No Flight Director). 
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Figure 6-37. Control Performance Cost vs Total Control Attention 
(Approach - z Flight Director). 
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Figure 6-38. Control Performance Cost vs Total Control Attention 
(Approach - x Flight Director), 
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Figure 6-39. Control Performance Cost vs Total Control Attention 
(Approach’ - Full Flight Director), 
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Figure 6-41 . RMS Performance vs System Automation 
(Full Flight Director). 
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6.4.4 ATTENTION ALLOCATION VS CONTROL AUTOMATION 


Figures 6-42 and 6-43 compare the fractional control attention for the various 
levels of control automation, \with no flight directors and with full flight directors 
respectively, and at a "comfortable" pilot workload (f^ = 0.4). These charts indicate 
the comparative importance of the instruments as the flight conditions change and/or 
as the system automation changes. For example, in Figure 6-43 the raw position 
instruments are used infrequently in all s^tems, particularly during the approach. As 
another example, adding pitch rate feedback (System B) and pitch attitude feedback 
(System C) allows the pilot to shift some attention from the attitude indicator to the 
vertical instruments (z and FD^), thereby improving his overall performance. 

6.4.5 ATTENTION ALLOCATION VS DISPLAY SOPHISTICATION 

Figures 6-44 through 6-50 show the change in fractional control attention with 
display sophistication for each system at f^. = 0.4. They illustrate the importance of the 
available instruments at a fixed workload for each level of control automation. The 
information level results demonstrate the importance of the rate information from each 
instrument. It is interesting that although the pilot spends a large fraction of his 
available attention to the z flight director (when provided), this does not achieve a 
significant improvement in performance. 

6.4.6 CONTROl/DISPLAY CONFIGURATION EVALUATION 

To evaluate the various control/display configurations, the workload for 
each can be normalized to a specified minimum acceptable performance. The 

performance cost for longitudinal control of any system is defined as 

6 - 

Jp = E ( a, /(x,)^,I^ (6-7) 
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Figure 6-43. Control Attention Allocation vs System Automation 
(Full Flight Director). 
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Figure 6-44. Control Attention Allocation vs Display Sophistication 
(System A). 


6-53 


AEROaPACe? SY8TEIVI8| ' INC* one vine brook park • Burlington. Massachusetts oisos • (sit) ste.tbit 












H = Hover 


A - Approach 


f =0.4 
c 



Display SophisHcation 


Figure 6-45. Confrol Attenfion AllocaHon vs Display Sophistication 
(System B). 
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Figure 6-46. Control Attention Allocation vs Display Sophistication 
(System C). 
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Figure 6-47. Control Attention Allocation vs Display Sophistication 
(System D). 
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Figure 6-48. Control Attention Allocation vs Display Sophistication 
(System F). 
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Figure 6-50. Control Attention Allocation vs Display Sophistication 
(System H). 
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where a are the predicted rms state errors (x, y, 0, etc.) and (x.) are the corre- 
X* I mox 

spending cost functional weightings given in Table 6-2, which correspond to 2a values. 
If we specify the minimum acceptable performance as 1 a errors in all states, i.e. 



(x.) 

' I 'max 


/2 


the maximum allowable value in Equation (6-7) is: 

(J ) = 60/2)2 = 1 .5 

^ max 


(6-8) 


(6-9) 


For the longitudinal control task, the maximum capability for both control and monitor- 
ing is about 0.6, and a "comfortable" value for control is f^ «. 0.4. By defining an. 
acceptable workload to be ± 0.1 from the comfortable level, the control/display con- 
figurations can all be classified into three categories (unacceptable, acceptable, 
excel lent) as shown in Figure 6-51. 

Applying this normalization technique to the candidate control/display 
systems for the CH-47, we obtain the performance evaluation shown in Figure 6-52 
for hover and approach. Besides quantifying the well-known conceptual plots of 
control automation vs display sophistication^. Figure 6-53 provides a great deal of 
informatiori on the various configurations. For example, System C is acceptable in 
both flight conditions if at least ah x flight director is provided. However, if the 
flight director should fail, the performance would be unacceptable in the hover 
condition . 


•f 

If the System D results are neglected, it is possible to sketch contours of constant 
workload as functions of control automation and display sophistication in Figure 6-52. 
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Figure 6-51 . Normalization of Performance and Workload Metrics. 
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Figure 6-52. 


Normalized Control Display Configuration Evaluation for CH-47. 
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Figure 6-53, Average Gain Flight Director Algorithms. 










6.4.7 


FIXED-GAIN FLIGHT DIRECTOR 


The flight director algorithms presented in Subsection 6.2 were designed sper 
cificolly for thegiven flight conditions. This implies the flight director computer would 
hove to vary all the gains as the flight conditions change.. As an alternative, a preliminary 
analysis of a fixed-gain flight director was conducted. The technique used was simply 
to average the "optimum” flight director algorithms for hover arid approach. The control 
cost weightings and indifference thresholds are based on the maximum permissible position 
and attitude errors, which are the approach values. The resulting flight director signals 
are shown in Figure 6-53 for Systems C and G. 

Figures 6-54 and 6-55 compare the respective system's performance at hover and 
approach. For comparison, the corresponding performance with no flight director and 
with the optimum (variable-gain) flight director is shown. Obviously, the optimum flight 
director is significantly better than the fixed-gain flight director for both systems. For 
System G, the fixed-gain flight director performance is about midway between the optimum 
flight director . and that with no flight director. However, in System C, the fixed-gain 
flight director is no better than having no flight director at all. 

■ Obviously, it is impossible to draw general conclusions from this brief exam- 
ination, but it does appear that some form of flight director gain adjustment or sensitivity 
switching may be required to meet the desired system performance throughout the flight 
regime. 
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Figure 6-55. Fixed-Gain Flight Director Performance (System G). 
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6.5 


MONITORING MODEL RESULTS 


Having established the necessary control workload and corresponding per- 
formance, we can turn to the pilot's monitoring capability. Referring to Figure 6-51 

the minimum control workload is determined by J =1.5. If this is less than f = 0.6, 

P ^ 

the difference between 0.6 and f is available for monitoring the status instruments, as 

c 

described in Section 2. Applying the monitoring model described in Section 4 to each 
acceptable system, we obtain the monitoring performance shown in Figure 6-56. It is 
interesting to note that monitoring performance does not always improve with increased 
automation. However, looking at the combined control and monitoring performance. 
Figure 6-57, we see that increased automation generally does improve total perfor- 
mance. 
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Figure 6-56. Monitoring Performance Cost Summary. 
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SECTION 7 


DISPLAY CONCEPT AND FORMAT 

Although no unique transformation from analytically determined information 
requirements to display layout is available, there are nonetheless a number of important 
design principles which must be taken into account. This is of particular importance 
when dealing with integrated displays and with the problem of VTOL control. 

Some of these principles, outlined in' Reference 20. and reiterated below, ore 
generally useful in the design of flight control displays.. They are not in any sense 
fundamental rules, and Individual display situations may call for variations in their 
application. These principles are supplementary to the generally accepted conventional 
display criteria relating to instrument design, including location, size, contrast, 
quantization, and display-control compatibility. 

• Operator Centered and Oriented Display - This is an extension of the 
inside-out display principle for integrated displays and favors a pre- 
sentation with the aircraft position and orientation fixed in the dis- 
play and the other pictorial information (horizon, glide -si ope, hover 
point, velocity impact point, altitude reference, etc.) moving with 
respect to this reference. 

• Geometric Real World Compatibility for Pictorial Displays - Although 
the Integrated display is not in general a contact analog (and 
typically includes command and/or situation information not present 
in the VFR view), such pictorial information as is present should be 
compatible with a view of the real world situation. The integrated 
disp ay is at its best when the information contained in it is perceived 
by the pilot as a single complex picture giving all the attitude and 
position information required, rather than as a densely packed code 
through which he can successively determine the aircraft flight path. 

• "Status at a Glance" for Situation Displays - In keeping with geometric 
real world compatibility, the essential elements of the display must 

be clearly delineated by size, shape, or color and coordinated with 
respect to one another so that the status of the aircraft, especially 
in unusual attitudes, is immediately obvious and does not require 
element decoding. 


7 - 1 

AEROSPACE SY8TEMS| INC* • one vine eRooK park • Burlington, Massachusetts oisoo • (017) S7a-7si7 


• Predictive Capacity - In addition to indicating the current state of the 
aircraft, the integrated display must readily show the dynamic situation 
so that the futpre state can be easily surmised. This kind of information 
is necessary for lead generation in fast loops such as attitude control and 
for planning maneuvers in guidance or collision avoidance. Display 
quickening, explicit rate symbols, display prediction, and historical 
trail markers may all be used to this end, and should follow the practice 
of derivative information "leading" the variable on the display. 

• Geometric Sensitivity and Scaling - The symbols and elements in an 
integrated display must move far enough and fast enough so that the 
pilot will be able to notice the motion and estimate Its magnitude in 
performing the appropriate control response. Maximum range and 
desired pilot gain in each loop must be considered In scaling the integrated 
display elements for various phases of flight. 

• Use of Digital Information Where Required, - An exception to the pictorial 
compatible principle is in the display of information which is slowly 
varying and which must be read accurately over a large range . In this 
case, the judicious use of some digital presentations on the integrated display 
is appropriate. As a simple example, altitude may be displayed digitally, 
although altitude deviation from a glideslope is best handled by analog 
motion of lines In the display. The amount of digital Information should 

be kept as small as possible, displayed only when necessary (perhaps on 
pilot demand), be legible, and contain as few digits as absolutely required. 


These design principles Were adhered to in the development of straw-man display format 
concepts for the implementation of Systems C and G designed and analyzed in the 
previous sections. 

7. 1 STRAW-MAN DISPLAY CONCEPT - PANEL LAYOUT 

The general layout of the display panel is illustrated in Figure 7-1. Details 
of the ADI and HSI differ between SystemsC and G. However, since both systems assume 
dependence upon o flight director, they are treated as integrated command displays 
rather than situation indicators. Although a full complement of flight instruments is 
anticipated, inthis discussion we will concentrate only on the instruments essential 
for the flight control functions considered in the study and will therefore limit the consid- 
eration to the quantities displayed on the ADI, HSI and altimeter. The assumption is 
made that the pilot derives vertical rate information from the altimeter rather than 
from an IVSI in both Systems C and G. 
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7.2 


ATTITUDE DIRECTOR INDICATOR (ADI) 


In Sysf-emC fhe pllof's control task consists of primarily following the flight 
director commands for aircraft attitude and collective pitch. The ADI format shown in 
Figure 7-2 is based on the well accepted moving aircraft symbol flight director and 
artificial horizon ball with modifications for VTOL applications. Indicated and 
commanded roll angle are scaled one to one. Pitch angle is scaled to permit display of +5® 
of pitch on the attitude ball permitting the pilot to resolve changes in pitch angle of 
less than 0.25 degrees. The turn rate iridicator at the bottom of the ADI has the 
capability of including a commanded turn rate independent of the commanded roll angle 
if so desired. It also indicates lateral acceleration in the conventional manner. The 
additional indicators on the right side of the ADI are for use in the commanded collect- 
ive for independent control of altitude. Consistent with the other elements of the 
flight director system is a pursyit display with commanded and actual collective stick 
indicated on either side of a fixed scale. 

The principal difference in the ADI format for System G is the attitude informa- 
tion displayed on the right side of the instrument. The ADI as used in System G becomes 
strictly a secondary monitoring instrument for attitude and turn rate. The flight director 
symbols are normally not displayed, and as envisioned in -System G, not utilized. The 
collective indicator on the right side of the ADI in System C is not required for System 
G . 

7. 3 HORIZONTAL SITUATION INDICATOR (HSJ) 

The HSI is discussed first for System G where it is used as one of the two ■ 
principal command instruments for the pilot's direction. In System C the HSI is used 
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Roll Angle 



Figure 7-2. Straw-Man Display Concept - ADI Layout for System C, 



as a secondary situation indicator only while the pilot flies the vehicle using attitude 
commands from the ADI, Consequently the HSI for System C will be simpler. 

The HSI format shown in Figure 7-3 is a heading-up mbvlng-map display with 
command information on required x and y velocity. The aircraft symbol in the center 
of the display is surrounded by range rings of a variable scale which may be set manually 
to a pre-selected sensitivity indicated by the number of feet for the first ring, or in an 
automatic mode. The automatic mode, selects the scale sensitivity which is the highest 
possible one that does not place the landing pad off the screen. The landing pad indicator 
which is used principally in the final phases of descent does not move entirely off' 
scale when a range sensitivity too great to permit it to be displayed is selected. Instead 
it is moved to the edge of the display and indicated with a semi-circle rather than a full 
circle, to indicate the relative direction to the pad. The compass rose fixed to the 
nroving map indicates heading in the conventional manner for System G. A command 
heading bug is also available as a supplementary flight director output. 

The pilot uses the HSI in System G principally as a horizontal velocity 
command display. The command velocity is shown as an open arrow with origin at the 
aircraft symbol, indicating both magnitude and direction of the commanded horizontal 
velocity. Full scale on the command velocity vector is taken as 100 kt. The actual 
velocity vector indicated by the solid arrow with origin at the aircraft symbol is the 
quantity directly controlled by the pilot's control stick in System G. To follow the flight 
director commands the pilot must line up the tip of the actual velocity vector with the 
command velocity vector. In this pursuit display the accuracy with which the align- 
ment can be resolved is of the order of 2 percent of full scale or approximately 3 ft/sec 
which is greater than the assumed indifference thresholds for velocity. Small errors in 
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velocity matching, however, are detected by the rate of change of the track position 
error crosspointers. 

When the aircraft is within 50 ft of the commanded position for hover or 
descent the crosspointers indicate cross-track and along-track position errors. Each 
ticmark is a nominal error of 25 ft. The crosspointers are not on the same scale 
as the moving map range rings and are used for a fine error adjustment on velocity 
control. They are consistent with the other flight director information in being 
a ”fly-to" display in which the aircraft symbol is flown towards the intersection of the 
crosspointers in order to null the error. Drift of the crosspointers is the most sensitive 
indication of a lack of perfect alignment of the command velocity with the actual 
velocity. 


7.4 ALTIMETER 

In System G the altimeter is one of the principal flight director instruments 
enabling the pilot to use the vertical command to achieve the altitude hold or altitude 
track required by the flight director. At the same time, however, the altimeter must 
serve as an accurate situation indicator to indicate actual altitude for monitoring pur- 
poses in both Systems C and G. The solution to these conflicting demands of range 
and precision is a type of combined digital and moving tape display as illustrated in 
Figure 7-4. 

The first three numerals are a digital display giving altitude In hundreds of 
feet. To the right is a moving tape giving altitude to the nearest foot. The portion 
of this tape centered within the viewing window is the actual altitude at the mid-point 
with a visible range of + 10 ft. This accuracy of altitude information is available from 
the radar altimeter in the hover case. To eliminate difficulty of reading the moving 
tape, the information will be filtered at a corner frequency of about 1 Hz prior to being 
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Figure 7-4. Straw-Man Display Concept - Altimeter Layout. 



displayed. Altitude rate is not assumed to be displayed explicitly on an IVSI in these 
systems/ but rather is assumed to be determined by the pilot implicitly from the altimeter. 
This is particularly easy for low altitude rates using the moving tape which moves upward 
or downward at a speed proportional to altitude rate. (Transfer of a numeral through the 
moving tape window in one second indicates a 20 ft/sec ascent or descent rate. Drift 
of less than 0.25 ft/sec will be observable). 

The commanded altitude is indicated In System G in two ways : by a 
numerical six-digit command which Is set either manually for an altitude hold or by the 
flight director system for control descent, and by a command altitude bug which moves 
up and down to the right of the altitude moving tape. The command bug is at the same 
scale as the moving tape altimeter when the commanded altitude is within +10 ft of 
the actual altitude; however, the. command altitude bug moves over a vastly reduced 
scale (approximately 100 ft per division) when the commanded altitude is more than 
10 ft above or below the actual altitude. The predictive altitude symbol shown by 
a solid circle is controlled directly by the vertical control stick in System G. The 
pilot's task is to place the predictive altitude indicator at the same point as the 
command altitude bug. In the situation indicated in Figure 7-4, for example, the 
pilot would press down on the vertical control stick to move the predictive altitude 
indicator down to the level of the command altitude bug which Is indicating a commanded 
altitude of 00225 ft. The predictive altitude and command altitude indicators would 
both move up until they are aligned with 00225 inthe center of the tape window. 

The altimeter layout is the same for System C with the exception that the 
predictive altitude Indicator, command altitude bug and command altitude numeric 
display are not present . 
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SECTION 8 


SUMMARY AND CONCLUSIONS 

This section presents a summary of the principal accomplishments achieved 
during this effort and a review of the most significant results. The next section will 
contain some suggestions for additional research and experimental investigations based 
on the study results. 

The primary accomplishment of this investigation has been the development 
of a systematic design methodology for pilot displays in advanced VTOL aircraft for 
commercial operations. This design approach accounts for various levels of control auto- 
mation and display sophistication. It is based on the optimal control model for the human 
operator, but includes several significant extensions in the state-of-the-art of pilot model- 
ing. An explicit attention allocation procedure has been established which determines 
the optimal division of the pilot's total attention between monitoring and control ta^ks, 
and among the various displays available to him for each task. 

The design methodology separates the model into three levels of detail. At 
the "information level," all of the state variables are assumed to be perfectly displayed 
to the pilot. Thus, the pilot has perfect knowledge of each state variable and his al- 
location of attention among these indicates their relevant importance in the ideal 
situation. At the display element level the effects of pilot indifference thresholds are 
introduced, and the pilot's ability to detect both position and rate from a given dis- 
play element is included. At this level the relative importance of each display ele- 
ment is determined, and a more realistic estimate of the overall system performance 
is obtained. Finally, at the "display format level" realistic performance estimates 
due to display thresholds, maximum deflections, instrument noise, scan frequency, etc. 
are determined for an actual display format which has been designed from the display 
element results. 
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The design mefhodology includes a model for simultaneous monitoring and 
control, which is based on the premise that the pilot first attempts to control the air- 
craft to a given level of performance, and then uses any additional capability for 
monitoring status information and/or automatic system performance. The model uses as a 
metric for control performance a quadratic function, of the state errors. The control work- 
load metric is the pilot's total control attention to all the displayed elements that is re- 
quired to achieve a desired level of system performance. The model optimizes the 
control performance metric by allocating his total control attention among the available 
displayed elements. Then his available attention for monitoring is determined as the 
difference between his total capacity and that required for control to the given perfor- 
mance level. The model next determines the optimum allocation of monitoring attention 
among the available status displays and the overall monitoring performance metric which 
is a quadratic index similar to the control performance metric. 

In summary, the VTOL display/control design: methodology is: 

• Determine (xj) and (J^) from mission requirements. 

max max 

• Select candidate control systems. 

• Ojiculate Jj. vs fg at the information level for each control 
system (include director/steering commands). 

• Choose display elements. 

: • Calculated^ vs f^ at the display element. 

• Determine f^ = fjoT ” ^c '^m monitoring model. 

• Select display/control system. 

• Select display format candidates. 

■ • Determine -^m format. 

• Select display format. 

A computer program entitled PIREP has been developed to implement the 
extended optimal control/monitoring model for the pilot. It can be used at the informa- 
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tion level/ at the display element level, or at the display format level, to determine 
the optimal* allocation of pilot's attention for either rnonitoring or control, as w6ll as 
the associated system performance. The principal inputs of PIREP are the system dynamics 
(automation level, external disturbances, etc,), the display model (display elements, 
threshold, etc,), and the total attention (control/monitoring). The primary outputs of 
PIREP are the optimum attention allocation (control/monitoring), the system performance 

metrics (J , J , cost gradients), and the rms predictions (state, display, control). 

c 'in 

The extended optimal control model for the pilot was validated by attempting 
to reproduce flight results obtained by NASA/LaRC with the CH-46 tandem'rotor heli- 
copter. Descriptions of the CH-46 model ^following control system, evaluation display 
panel, and the flight director algorithms were obtained from NASA publications. The 
optimal control model was exercised at the display format level for a hover flight condi|'ion 
and the results were compared with limited flight data. Both the analytical and experi- 
mental results show that the pilot could not. adequately hover without the flight director, 
but that he had very little difficulty with the flight director. 

A flight director design technique using quadratic synthesis was developed 
as a straightforward means of generating flight director algorithms. These algorithms 
were designed to relate to the pilot task objectives, to minimize his workload and/pr 
improve his control performance, and to satisfy the pilot's desired goal of behaving ap- 
proximately as a gain and time delay. The flight director signals are obtained as linear 
functions of the system states as a direct fallout from, the optimal control model . When 
applied to the CH-46 helicopter, the flight director design technique produced nearly 
identical time constants to those of the empirically determined flight director algorithms 
used by NASA/LaRC . 

A similar approach using quadratic synthesis was applied to determine flight 
control automation levels for the helicopter. By appropriately specifying the control 
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and state weights in a, quadratic performance index, various levels of automatic feed- 
back. control systems can be systematically designed. These range from a totally manual 
basic vehicle with no feedback to the fully automatic system with complete position 
feedback . 

The display/control design methodology was then applied to predict the 
performance of the LaRC GH-47 helicopter which will be used as the VALT research 
aircraft. ’ Two flight conditions were investigated, hover at sea level and an -approach- 
ing condition at 60 knots and 1,000 ft/min. descent. Seven levels of control automa- 
tion were considered: . . 

• Basic system without any feedback 

• Pitch -rate feedback 

• Pitch-attitude feedback 

• Pitch-attitude and vertical speed feedback 

• Pitch-attitude and vertical position feedback 
0 Forward-speed and vertical position feedback 
0 Full position feedback. , 

Five display system levels were also considered: 

0 Information level 

• Element level with no flight director 

• Element level with vertical flight director 

• Element level with horizontal flight director 

0 Element level with vertical and horizontal flight director. 

Cost weighting functions and indifference thresholds for the CH-47 were selected based on 
the desired performance requirements for an advanced VTOL commercial helicopter. 

In general, the numerical results indicated that the flight director does improve 
system performance. Although this is a fairly obvious and expected result, the model does 

8-4 

AEROSPACE SYSTEMS^ ll\IC» • ons vine brook park • surlington, Massachusetts oieoa • (ai7) sts-tsit 


provide quantitative indications of the performance improvement with the flight di- 
rector. The results also showed that the vertical flight director provides marginal ; 
performance, improvement; most of the performance gain is produced by the forward 
flight director. The CH-47 results also showed that control automation generally 
improves performance. Again, this.is aniobvious and expected conclusion, but the ■ 
model provides quantitative measures of the performance improvement for various auto- 
mation systems. Moreover, In order to achieve the desired system performance, some . ■ 
leyel of ..automation will be. required for most advanced VTOL missions. The numerical 
results showed that at hover and approach the CH-47 helicopter without control' 
automation cannot be flqwn to an acceptable performance level. They also indicoted 
that the hover condition Is considerably more, difficult than the approach. However, 
increasing system automation tends to reduce the difference in difficulty between.the 
two flight conditions. Also, the more automatic systems tend to be less sensitive to 
pilot workload variations; as automation increases, the slope of the performance curve 
versus workload is lower. This means that other temporary demands on the pilot's at- 
tention will cause less deterioration In system performance as automation increases. 

A.fixed-gqin flight director was used to examine the effects of changing ■ 
flight conditions on system performance. This flight director used the average of the 
hover and approach gains and was examined at both flight conditions for two automa- 
tion levels. At the lower automation level, the fixed-gain flight director was no better 
(perhaps slightly worse) than no flight director at all; whereas in the other case, the 
fixed-gain flight director performance lay about mid-way between the optimum flight 
director and no flight director. From this limited examination, it appears that in 
general the flight director gains will probably have to be adjustable to handle a wide 
range of flight conditions. 
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The relaHve im porta nee of the individual display elements was clearly 
demonstrated -for all display sophistication and control automation levels. The model 
provides a quantitative measure of the relative importance of each display element by 
means of the optimum attention allocation. For example, at a given level of control 
automation and control attention, the information level results show that the pilot is 
far more interested in forward and vertical velocities and pitch rate than he is in the 
position and attitude elements. As the display sophistication increases ~ that is as ver- 
tical, forward, and full-flight director signals are jsfovided — he spends less attention 
on the situation displays and more on the flight director signals. Sirnilarly, as system 
automation increases for a given display configuration, the pilot adjusts his attention 
accordingly. For example, as pitch-rate feedback and pitch-attitude feedback are 
included, he spends less time on the pitch display and more time on the position display. 

The monitoring model confirms the a priori conjecture that more monitoring 
generally improves system performance, since more monitoring time implies less control 
workload. However, it is interesting to note that monitoring performance itself does 
not necessarily improve either with increased system automation or with display 
sophistication. Several questions still remain regarding the interpreting of the monitor- 
ing model and merging the monitoring results with the control results. 

The actual design of a display format is still far more of an art than a science. 
However, there are several design principles that should be used to simplify the 
translation of display element analytical results to the instrument format: 

• Operator Centered and Oriented Displays 

• Geometric Real-World Compatibility for Pictorial Displays 

• "Status at a Glance" for Situation Displays 

• Predictive Capacity 
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• Geometric Sensitivity and Scaling 

• Digital Information Where. Required 

Using these principles, a straw-man display concept was designed in an attempt to 
satisfy the results of the optimal control/monitor model for the CH-47 with two levels 
of control automation. These control systems were pitch-attitude feedback with no 
vertical feedback and velocity command with altitude hold. The three primary instru- 
ments of interest for both cases are: the attitude director indicator (ADI), the horizontal 
situation indicator (HSI), and the altimeter. For the pitch-attitude feedback system, 
the ADI provides all four flight director commands in addition to pitch and roll attitude 
information, while the HSI and the altimeter provide accurate situation information. 

For the velocity command system, however, the ADI is primqrily a monitoring instrument 
and provides no flight director commands, while the HSI and the altimeter provide the 
principal flight director commands as well as situation information. 
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SECTION 9 
RECOMMENDATIONS 

During the course of f-he study, several areas hove been identified in which • 
additional research is needed to define, test and validate the use of the pilot In an 
automated VTOL aircraft operated as a. short-haul, commercial aircraft. These subject 
areas are outlined below. 

9.1 COMBINED LATERAL AND LONGITUDINAL ANALYSIS 

The pilot/vehicle/control system models, computer programs and design 
procedures developed in the study should be utilized to conduct a combined lateral 
and longitudinal analysis of pilot/automatic system task allocation for the VALT Project 
CH-47 helicopter. Such an analysis wiM extend the' results obtained in the present 
study and will investigate.a number of flight conditions including hover, straight-in 
approach, turning approach and cruise phases of flight. Effects of wind gusts and . , • 

wind shears should be Included in the analysis. Different levels of control automation 
ranging from fully automatic to fully manual should be investigated. Consideration 
should be. given, to use of an advanced stability augmentation system such as that 
designed for the VALT Project CH-47 helicopter. The work effort under this task 
would culminate in a recommended display/control system concept for the CH-47. 

9.2 MONITORING MODEL INVESTIGATION 

The investigation of the monitoring model developed in the study for 
varying levels pf automation should be extended. This effort would utilize the model 
to determine the monitoring attention allocation to each of the primary instruments, 
and to the longitudinal and lateral axes. Analysis and interpretation of the monitoring 
model results obtained should then be conducted to determine 1) the effects of different 
monitoring strategies (i.e., current status or future status) on attention allocation and 
2) the extent to which attention used for control (f^) can also be used for monitoring. 
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9.3 


EFFECTS OF SYSTEM FAILURES ON THE PILOTING TASK 


The effects of system failures on the piloting task should be investigated. 
Representative failure modes should be defined (i.e., actuator failure, sensor failure, 
etc.) consistent with the generic augmentation systems. The optimal. control model 
developed in the study should then be utilized to examine the, pilot control workload 
and performance under the assumed failure mode conditions. For example, results . 
could be obtained that indicate whether a flight director designed for a velocity 
command system can be used when that augmentation system fails. Consideration should 
be given to investigation of transient conditions to determine the time required by the 
pilot to recover from different assumed failures. 

9.4 FIXED-BASE SIMULATION EXPERIMENTS 

A series of fixed-base pilot-in-the-loop simulation experiments should be 
planned and conducted to T) validate the extensions in the pilot model accomplished 
during the initial phase of the study and 2) evaluate and verify the display/control 
system concept for the CH-47. The experiments should be conducted initially on a 
relatively simple and inexpensive interactive display facility such as a PDP-ll/lO 
and on the NASA VALT fixed-base display research facility with a good cross section 
of subject pilots. Consideration should be given to experiments for measuring the 
performance differences and subjective differences between integrated displays and 
separate displays. Based on the results of the experiments, methods for representing 
appropriate changes in the analytical model should be made. Experiments to measure 
monitoring strategies in the simulation should be included to determine whether the 
pilot actually uses a normative strategy (i.e., how he should do it) as represented in 
the present model . 
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9.5 


PILOT INTERACTION WITH AUTOMATIC SYSTEMS 


The pilot's interaction with on automatic system has been only briefly 
addressed by control theory models. Studies of this interaction should be conducted 
to include factors such as pilot acceptance of the system; the "harmony" of the system 
(does it respond the way the pilot thinks he would?); does the pilot interact with the 
system by monitoring a closed-loop system or by "controlling" the open-loop system? 
Actuator movement information can be presented to the pilot by dedicated displays, 
or by control stick motion. In one case only visual information is presented (actuator 
monitoring is done through the visual channel); and in the other case the monitoring is 
done through the kinesthetic channel (thus allowing more time to visually monitor 
other displays). The model should be examined to see whether differences in failure 
detection times using these two methods can be accurately represented by the model. 


9-3 


AEROSPACE SYSTEIVISi IIVC* • one vine BROOK park • Burlington, Massachusetts 01B03 • (si7) 27S-7S17 


REFERENCES 


1. Hoffman, W. C.; Hollister, W. M.; and Howell, J. D.: Navigation and Guid- 
ance Requirements for Commercial VTOL Operations. Aerospace Systems, Inc. 
Interim Technical Report, ASI-TR-74-17, NASA CR 132423, January 1974. 

2., Hoffman, W. C.; and Hollister, W. M.: Navigation and Guidance Require- 
ments for Commercial VTOL Operations - Navigation System Sensitivity Analysis. 
Aerospace Systems,. Inc, Interim Technical Report, ASI-TR-75-25, May 1975. 

3. Hoffman, W.C.; and Hollister, W. M.: A Spiral Guidance Approach Concept 
for Commercial VTOL Operations. Aerospace Systems, Inc. Technical Report, 
ASI-TR-75-21, March 1975. 

4 . Curry, R. E.; Kleinman, D, L.; and Hoffman, W. C.: A Model for Simultaneous 
Monitoring arid Control. Paper presented at Eleventh Annual Conference on 
Manual Control, NASA Ames Research Center, May 1975. 

5. Kleinman, D. L.; Baron, S.; and Levison, W. H.: An Optimal Control Model 
of Human Response, Parti: Theory and Validation. Automatica, 1970, Vol . 6, 
pp. 357-369. ' 

6. Kleinman, D. L.; Baron, S.; and Levison, W«. H..: A Control Theoretic Approach 
to Manned-Vehicle Systems Analysis. IEEE Trans. Auto. Control. Vol. AC-16, 
No. 6, December 1971, pp. 824-832. 

7. Levison, W. H.: A Model for Task Interference. Proceedings of the Sixth At 
nual Conference on Manual Control, April 1970, pp. 585-616. 

8. Weir, D.H.; Klein, R. H.; and McRuer, D. T.: Principles for the Design of Ad- 
vanced Flight Director Systems Based on the Theory of Manual Control Displays. 
NASA Cr- 1748, March 1971. 

9. Baron, S.; and Levison, W. H.: A Manual Control Theory Analysis of Vertical 
Situation Displays for STDL Aircraft. Bolt Beranek and Newman, Inc. Report 
No. 2484, April, 1973. 

10. Chen,; Lebacqz,; and Aiken: A Preliminary Look at Flight Director Design Phil- 
osophies for Application to a VTOL Landing Approach Flight Experiment. Pro- 
ceedings of the Tenth Annual Conference on Manual Control, April 9-11, 1974. 

n. Kleinman, D. L.; and Baron, S.: Manned-Vehicle Systems Analysis by Means 
of Modern Control Theory. NASA CR- 1753, June 1971. 

12. Garren, J. F., Jr.f Kelly, J. R.; Sommer, R. W.; and DiCarlo, D. J.: Flight 
Investigation of VTOL Control and Display Concept for Performing Decelerating 
Approaches to an Instrument Hover. NASA TN D-6108, February 1971. 

13. Senders, J. W.: The Human Operator as a Monitor and Controller of Multi- 
Degree of Freedom Systems. IEEE Trans, on Human Factors in Electronics, HFE-5, 
No. 1, September 1964. 


R ^ 1 


AEDOaPACE SYSTEIVISi IIMC* • one vine brook park • burunoton. Massachusetts oiBoa • (017)273-7317 


14. Senders, Ward, J. L.; and Carbonell, J. R.:^. Human Visual Sampling Pro- 

cesses: A Simulation Validation Study. NASA-CR-1258, January 1969. 

15. Senders, J. W.; Ellcind, J. 1.; Grignetti, M C.; and Smallwood, R.: An Invest- 
igation of the Visual Sampling Behavior of Human Observers. NASA Report No. 
CR-434, April 1966. 

16. Gai, E.: Physchophysical Models of Signal Detection with Time Varying Stimuli. 

Ph. D. Thesis, Dept, of Aeronautics and Astronautics, M.I.T., 1974. 

17. Kelly, J . R.; Niessen, F. R.; and Sommer, R. W.: Evaluation of a VTOL Flight- 
Director Concept During ConstantrSpeed Instrument Approaches. NASA TN D-5860, 
June 1970. 

18. Kelly, J. R.; Niessen, F. R.; Thibodeaux, J. J.; Yenni, K. R.; and Garren, 

J. F., Jr.: Flight Investigation of Manual and Automatic VTOL Decelerating In- 
strument Approaches and Landings. NASA TN D-7524, July 1974. 

19. Kleinman, D. L.; and Killingsworth, W.: A Predictive Pilot Model for STOL Air- 
craft Landing. NASA CR-2374, March 1974. 

20. Young, L. R.: Integrated Display Principles and Some Applications to V/STOL Air- 
craft. AGARD Conference Proceedings No. 96 on Guidance and Control Displays, 
Oct. 1971. 

21 , Kleinman, D, L.: Optimal Control of Linear Systems with Time Delay and Observa- 
tion Noise. IEEE Trans. Autom. Control, Vol. AC-14, No. 5, October 1969. 

22. Kleinman, D. L.: Suboptimal Design of Linear Regulator Systems. MIT-ESL-R-297, 
February 1967. 

23. Bartels, R. H.; and Stewart, G. W.: Solution of the Matrix Equation AX + XB = C. 
Comm, of ACM, Vol. 15, No. 9, September 1972. 

24. Press, H.; Meadows, M. J.; and Hodlack, I.: A Reevaluation of Data on Atmo- 
spheric Turbulence and Airplane Gusts Loads for Application in Spectral Calcula- 
tions. NACA Report 1272, 1956. 

25. Schlundt, R. W., et al.: STOL Traffic, Environment and Operational Procedures. 
R-717, C. S. Draper Laboratory, March 1972. 


R -2 


AEROSPACE SYSTEM S| live. • one vine brook park • Burlington. Massachusetts oieo 3 • (017)273-7317 


BIBLIOGRAPHY 


1. Adams, G. D.: V/STOL Approach System Flight Tests: Final Report. FA A Report 
No. RD -66-95, December 1966. 

2. Allen, R. W.; Clement, W. F.; and Jex, H. R.: Research on Display Scanning, 
Sampling, and Reconstruction Using Separate Main and Secondary Tracking Tasks. 
NASA Report CR-1569, Washington, D. C., July 1970. 

3. Anon.: Space Vehicle Displays Design Criteria. NASA Space Vehicle Design 
Criteria (Guidance and Control), NASA SP-8086, March 1972. 

4. Anon.: Theoretische und experimental le Untersuchungen von Verfahren ein 
Kombiniertes Navigations- und Allwetter-Landesystem (KNL) fur V/STOL Flugzeuge. 
Teldix Luftfahrt-Ausrustungs-GmbH , Heidelberg, Germany, June 1967. 

5. Anon.: V/STOL Displays for Approach and Landing. AGARD Report 594, 1972. 

6. Ashkenas, I. L.; and McRuer, D. T.: A Theory of Handling Qualities Derived from 
Pilot-Vehicle System Considerations. Aerospace Eng., 1962, Vol. 21, No. 2, 
pp. 60 , 61, 83-102. 

7. Baron, S.: Analysis of Response to Wind-Shears Using the Optimal Control Model 
of the Human Operator. Proceedings of 9th Annual Conference on Manual Control, 
May 1973, pp. 419-428. 

8. . Baron, S., et al.: Application of Optimal Control Theory to the Prediction of 

Human Performance in a Complex Task. AFFDL TR-69-81 , March 1970. 

9. Baron, S.; and Kleinman, D. L.: The Human as an Optimal Controller and Infor- 
mation Processor. IEEE Trans. Man -Machine Sys., Vol. MMS-10, March 1969, 
pp. 9-17. 

10. Baron, S.; and Levison, W. H.: A Display Evaluation Methodology Applied to 
Vertical Situation Displays. Proceedings of 9th Annual Conference on Manual 
Control, May 1973, pp. 121-132. 

11. Bergstrom, B.: Interpretability Studies of Electronic Flight Instruments. Saab 
Aktiebolag, .Linkoplng, Sweden. Saab TN 61, January 1967. 

12. Bondurant, R. A., Ill; and Kearns, J. H., Ill: V/STOL Vertical Situation Display. 
Paper presented at the SAE National Aeronautics and Space Engineering Meetings, 
Los Angeles, CA, Society of Automotive Engineers. Paper No. 690694, October 
1969. 

13. Clement, W. F.; Hofman, L. G.; and Graham, D.: A Direct Procedure for 
Partitioning Scanning Workload in Control Display System Design. Proceedings of 
9th Annual Conference on Manual Control, May 1973, pp. 389-400. 


b - 1 


AER08PACB 8 YSTEMB 1 INC* • 0N6 vine brook park • Burlington, Massachusetts 01803 • (ei 7 ) a 7 a- 73 iT 



14. 


Clement, W, F.; and Hofinan, L. G.: A Systems Analysis of Manual Control 
Techniques and Display Arrangements for Instrument Landing Approaches in 
Helicopters, Vol. I: Speed and Height Regulation. JANAIR Report 690718, 
Systems Technology, Inc. Tech. Report 183-1, July 1969. 

15. Clement, W. F,; McRuer, D. T.; and Klein, R, H.: Systematic Manual 
Control Display Design. AGARD Conference Proceedings No, 96 on Guidance 
and Control Displays, October 1971. 

16. Dunham, R. E.; and Sommer, R, W.: Evaluation of a Moving-Graph Instrument 
Display for Landing Approaches With a Helicopter. NASA TN D-6025, 

October 1970. 

17. Elkind, J . 1 Characteristics of Simple Manual Control Systems. MIT Lincoln 
Laboratory, Tech. Report III, April 1956. 

18. Etkin, B.: Dynamics of Flight: Chapter 10, Flight in Turbulent Air. Wiley & 
Sons, New York, 1962. 

19. Garren, J. F., Jr.; and Kelly, J. R.: Description of an Analog Computer 
Approach to V/STOL Simulation Employing a Variable-Stability Helicopter. 
NASA TN D-1970, January 1964. 

20. Garren, J. F.: VTOL Flight Investigation to Develop a Decelerating Instrument 
Approach Capability. SAE Paper 690693, October 1969. 

21. Garren, J. F., et al.: Flight Investigation of VTOL Control and Display Con- 
cept for Performing Decelerating Approaches to an Instrument Hover. National 
Aeronautics and Space Administration, NASA Technical Note TN D-6108, 
February 1971. 

22. Greicey; W., et al.: Evaluation of Cross-Pointer-Type Instrument Display in 
Landing Approaches With a Helicopter. NASA TN D-3677, November 1966. 

23. Gracey, W.; Sommer, R, W.; and Tibbs, D. F.: Evaluation of a Moving-Map 
Instrument Display in Landing Approaches With a Helicopter. NASA TN D- 
3986, May 1967. 

24. Gracey, W.; Sommer, R. W.; and Tibbs, D. F.: Evaluation of a Closed- 
Circuit Television Display in Landing Operations With a Helicopter. NASA 
TN D-4313, February 1968. 

25. Gracey, W.: Comparison of Information Display Concepts for Landing of VTOL 
Aircraft. NASA Technical Note TN D-4861, November 1968. 

26. Hoffman, W. C.; Zvara, J.; and Bryson, A. E.: Automatic Guidance Concept 
for VTOL Aircraft. J. Aircraft, Vol. 8, No. 8, August 1971, pp. 637-642. 

27. Hofiman, L. G.; Clement, W. F.; and Blodgett, R. E.: Further Examination of 
Pilot Instrument Scanning Data and Development of a New Link Value Estimator. 
Proceedings of 9th Annual Conference on Manual Control, May 1973, pp. 379- 
388. 


b -2 


aerospace SYSTEIVISt IIMC. • one vine brook park • Burlington. Massachusetts oieoa • (si7) 


28. Hollister, W. M,, et al,; Effect of the City Center Environment on VTOL 
Terminal Operations. American Helicopter Society Preprint No. SW-70-48, 
Presented at the Joint Symposium on Environmental Effects on VTOL Designs, 
Arlington, Texas, November 1970, 

29. Illingworth, J, K. B.; and Shazler, J, E,: Blink Flying of V/STOL Aircraft, 
With Particular Reference to Takeoff and Landing, AGARD Report, June 1960, 

30. Kelly, J, R.; Garren, J. F,, Jr.; and Deal, P. L,: Flight Investigation of 
V/STOL Height-Control Requirements for Hovering and Low-Speed Flight Under 
Visual Conditions. NASA TN D-3977,'May 1967. 

31. Ketchel, J. M,; and Jenny, L, L.; Electronic and Optically Generated Air- 
craft Displays: A Study of Standardization Requirements. JANAIR Report No. 
680505, May 1968. 

32. Kleinman, D. L.: Optimal Control of Linear Systems With Time-Delay and 
Observation Noise. IEEE Trans. Auto. Control, Vol. AC-14, October 1969, 
pp. 524-527. 

33. Kleinman, D, L.; and Baron, S,: Manned Vehicle Systems Analysis by Means 
of Modern Control Theory. Bolt, Beranek and Newman, Inc., BBN Report 
1967, June 1970. 

34. Kornstadt, H. J.; and Pfennigstorf, J.: Evaluation of an Integrated Flight 
Display for the Manual IFR Landing of VTOL Aircraft. AGARD Conference 
Proceedings No. 96 on Guidance and Control Displays, October 1971 . 

35. . Levison, W, H.; Baron, S.; and Kleinman, D. K.: A Model for Human 

Controller Remnant. IEEE Trans. Man-Mach. Syst., 1969, MMS-10, pp. 101- 
108. 

36. Levison, W. H.: The Effects of Display Gain and Signal Bandwidth on Human 
Controller Remnant, Report No. AMRL-TR-70-93, Aerospace Medical Research 
Laboratory, Wright-Patterson Air Force Base, Dayton, Ohio, March 1971. 

37. Levison, W. H.; Elkind, J. I,; and Ward, J. L.: Studies of Multivariable 
Manual Control Systems: A Model for Task Interference. NASA Report CR- 
1746, Washington, DC, May 1971. 

38. , Levison, W. H.; and Tanner, R. B.: A Control -Theory Model for Human 

Decision-Making. NASA CR-1 953, December 1971 . 

39. McRuer, D. T.; and Jex, H. R.: A Review of Quasi -Linear Pilot Models. 

IEEE Trans. Hufian Factors Electronics, 1967, HFE-8, pp. 231-249. 

40. McRuer, D.; and Weir, D. H.: Theory of Manual Vehicular Control . IEEE 
Trans. Man-Mach. Syst., 1969, MMS-10, Vol. 4, pp. 257-291. 


b - 3 


aerospace SYBTEIVISi IIMC* • ONE VINE BROOK PARK 


• BURLINGTON* MASSACHUSETTS 01B03 • (S17) 372-7017 


41 . Murphy, M. R., et al.: Simulator Evaluation of Three Situation and Guidance 
Displays for V/STOL Zero-Zero Landings. Presented at 10th Annual Conference 
on Manual Control (Wright-Patterson Air Force Base), April 1974. 

42. NASA Langley Research Center: VTOL & STOL Technology in Review. Astro- 
nautics & Aeronautics, Vol. 6, No. 9, September 1968, pp. 56-67. 

43. Reeder, J. P.: The Impact of V/STOL Aircraft on Instrument Weather Operations. 
NASA TN D-2702, February 1965. 

44. Reeder, J. P.: The Changing World of the Pilot in the Next Decade - V/STOL 
Aircraft Operations. Paper presented at the Human Factors Society National 
Symposium, San Francisco, CA, October 1970. 

45. Ringland, R. F.; Stapleford, R. L.; and Magdaleno, R. E.: Motion Effects on 
an IFR Hovering Task - Analytical Predictions and Experimental Results. Systems 
Technology, Inc. Tech. Report 188-1, October 1970. 

46. Russel, L.: Characteristics of the Human as a Linear Servo Element. M.D. 

/ Thesis, MIT Dept. Elec. Eng., May 1951. 

47. Seckel, E., et al.: A Note on the Effect of Helicopter Dynamics on Steep 
Instrument Approaches. Report No. 600, Department of Aeronautical Engineer- 
ing, Princeton University, February 1962. 

48. Simpson, R. W.; and Miller, R. H.: V/STOL in the Northeast Corridor. 
Astronautics & Aeronautics, September 1968. 

49. Simpson, R. W.: The Outlook for Future Commercial VTOL Transportation. 

AIAA Paper No. 69-198, AIAA/AHS VTOL Research, Design, and Operations 
Meetings, Georgia Institute of Technology, Atlanta, Georgia, February 17-19, 
1969. 

50. Sommer, R. W.; and Dunham, R. E.: Evaluation of a Contact-Analog Display 
in Landing Approaches With a Helicopter. NASA TN D-5241, June 1969. 

51. Stepniewski, W. Z.; and Schmitz, F. H.: Possibilities and Problems of 
Achieving Community Noise Acceptance of VTOL. Presented at the Eighth 
Congress of the I CAS, Amsterdam, September 1972. 

52. Sweeney, J. S.; and Runner, G.: A Comparison of Quickened and Unquickened 
Displays for the Monitoring of Vehicle Performance Under Full Automatic Control. 
Proceedings of the IEEE International Congress on Human Factors in Electronics, 
Long Beach, Mqy 1962. 

53. Trant, J. P.; and Algranti, J. S.: Investigation of VTOL Approach Methods by 
Use of Ground-Control led-Approach Procedures. NASA TN D-1489, October 
1962. 

54. Van Houtte, N. A. J.: Display Instrumentation for V/STOL Aircraft in Landing. 
Ph.D. Thesis, MIT Dept, of Aeronautics and Astronautics, June 1970. 


b -4 


AEROSPACE BVSTEMSi INC. • one vine brook PARK • Burlington, Massachusetts oisoa • (si7) a7s-7Bi7 


55. Vreuls, D., et a!.: All-Weather Landing Flight Director and Fault Warning 
Display Simulator Studies. Proceedings of Fifth Annual Symposium: Human 
Factors in Aviation, 1968 (Los Angeles), June 1968. 

56. Walters, D. J.; and Boyer, R.: V/STOL Displays for Approach and Landing. 
AGARD Conference Proceedings No. 96 on Guidance and Control Displays, 
October 1971 . 

57. Weir, D. H.; and Johnson, W. A.: Pilot Dynamic Response to Sudden Flight 
Control System Failures and Implications for Design. NASA Report CR-1087, 
Washington, DC, 1968. 

58. Winn, A. L.; and Lewis, R. B., II: Pilot Workload During Instrument Flight. 
American Helicopter Society, Preprint No. 820, May 1974. 

59. Winter, F. J., Jr.: Controls and Displays for Helicopter IFR Operation - Pilot 
Factor Considerations. American Helicopter Society, Preprint No. 825, May 
1974. 

60. Wolf, J. D.; and Barrett, M. F.: IFR Steep -Angle Approach: Effects of System 
Noise and Aircraft Control -Augmentation Variables. JANAIR Report No. 
700810, April 1971 (AD 724 336). 

61 . Wolf, J. D.: Display and Related System Requirements for IFR Steep Approach. 
Honeywell, Inc. (St. Paul), AD-736 247, January 1972. 

62. Wolkovitch, J.; LaMont, C. W.; and Lochtle, D. W.: Fundamental Limitations 
on V/STOL Terminal Guidance Due to Aircraft Characteristics. NASA CR-1901, 
National Aeronautics and Space Administration, Washington, DC, December 
1971. 

63. Youens, F. P.: Approach and Landing Prpblems in Jet VTOL Aircraft. AGARD 
Report 489, October 1964. 

64. Young, L. R., et al.: The Adaptive Dynamic Response Characteristics of the 
Human Operator in Simple Manual Control. NASA TN D-2255, April 1964. 


b -5 


AEROSPACE 8Y8TEMS| INC* • one vine brook park • Burlington, Massachusetts oieoa • (ai7) S7S- 


APPENDIX A 


EQUATIONS FOR THE OPTIMAL CONTROL 
MODEL FOR THE HUMAN OPERATOR 

A.l BACKGROUND 

The human operafor's basic task is to control, in some prescribed way, a 
dynamical system that is subject to external random disturbances. It is assumed 
that the system dynamics, which may include actuator, sensor, and noise dynamics, 
are described by the linear, time invariant equations of motion. 

x(t) = Ax(t) + Bu(t) + w(t) 

where the n-vector x(t) represents the vehicle state, u(t) is the human's control 
r.-input to the system, and where w(t) represents random external disturbances, 
w(t) is assumed to be a vector of independent zero-mean, Gaussian white noises 
with autocovariance 

E {w(t)w'(cr)} = Wfi(t - O’). 

In controlling the vehicle, the human perceives a delayed, noisy replica of the m 
displayed outputs y = Cx(t) + Du(t). Thus, 

yp(0 = Cx(t - r) + Du(t - t) + Vy(t - t). 

It is assumed that the control task is adequately reflected in the human's choice of 
a feedback control u(-) which, in the steady-state, minimizes the general quadratic 
cost functional 


A - 1 


(A-1) 
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J(u) = E 


lim 1 

T-*eo J 


y'(0Qyy(0 + 


u(t)Q^u(Odt| 


(A-4) 


If has been shown that fhe control which minimizes J(u), conditioned on the observations 

y (‘)f is generated by the linear feedback law 

P ’ - 


u(t) + u(t) = - Lx(t) =u^(t) 


(A-5) 


where x(t) is the best estimate of the system state x(t) based on the observed data 
yp(o’)/0’< t. The matrices Tj^ and L are obtained as 

I’m " ^22 ®r 


1 = p“^ P' 

L ^22 t'j2 
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where 


P = 
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-D 

^ 2 ' 

P* ' 

^12 * 
L 1 

P 22 


satisfies the Ricatti equation 
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and C = [C ; D] 
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Nofe that the Control Equation (A-5) can also be written as 


where 

I, = Q-' = [q;'p’,2 ; q;'p J = [rj l | rj 

The human operator injects noise into the generation of the control input. 
In the model, a Gaussian white noise v^(t) is added directly to the "commanded 
control" u^(t) in Equation (A-5) with covariance 

E IV') v;„W| «(»-■’) (A-8) 

Thus the human's control input is assumed to be generated by 

T j^u(t) + u(t) = u^(t) + v^(t) (A-9) 

with 


u^(t) = -L x(t) 


The estimate x(t) is obtained from the cascade combination of a Kalman 
filter and predictor. Define the augmented state vector x(t) = 
new x(t) satisfies 


[u(o] 


x(t) - Ax(t) + B u^(t) + w(t) 


(A-10) 


where 
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The Kalman filter generates x(t - t), the I east-mean -square estimate of the delayed state 
from 


-^x(t - t) = /S^(t - t) + 2 C Vy Vp W " Cx(t - t)1 + Bu^(t - t) 


where V is the covariance of the observation noise v (0# and S satisfies 
y y 




0 = AZ + 2A' +W -2C V"’ C2 


(A-n) 


(A-12) 
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t“?' 


' N m 
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W = 


The predictor generates the best estimate x(t) from the estimator output x(t - t) 
according to the equation 

x(t) = e^’^ x(t - r) +y ~ B u^ (a) d (T (A-13) 

t-T 

A differential equation for x(t) can be derived from Equation (A-13), 


^i(t) = Ax(t) + e'^'^ZCNT^ [yp(0 - Cx (t - t)1 




+ Bu (t - t) 
c' ' 


(A-14) 


A. 2 STATISTICAL CHARACTERISTICS 

A closed form expression for the covariance of x(t) may be derived 
(Reference 21) as 
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(A-I5) 


ow T ^ 

X = E |x(t)x'(t) I + J' 


+ f e^°’e^^C'V"^ C2e^‘^ da 
0 ^ 

— rs-» />-> 

where A = A - BL = A^ - is the closed-loop system matrix. 

There are three terms in X ; they arise from writing 


A 

X (t) = Ef(t)+Ep(t) + X(t) 

A 

where Ej.(t) = filtering error, ~ prediction error and X(t) = “pilots" state 

estimate. The covariance matrices for these terms are all affected by the time delay 
with 

= ; i = 6a (A-16) 

o 

A 

and X = E|x(t) x'(t)| is given by the last integral term in Equation (A-15). Note 
that X satisfies 

^ A rs^ . 

AX + XA' + e^^ZC'V“^C2e^'’’ = 0 (A-17) 
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APPENDIX B 


ATTENTIONAL ALLOCATION USING THE 
OPTIMAL CONTROL MODEL 


The basic approach that we follow in applying the optimal control pilot model 
to predict attention allocation among a set of display indicators is to optimize a quadratic 
cost functional with respect to pilot attentional constraints. In order for the entire scheme 
to be computationally attractive, the following must be accomplished: 

• Relate attentional model parameters (f.) to pilot model parameters. . 

• Obtain an expression for J* = minimum J that shows explicitly how the 
f. affect the various cost functional terms. 

• Obtain an expression for the gradient terms that will be needed in 

subsequent optimization algorithms. ^^i 

• Develop an algorithm to minimize J* with respect to the f., subject to 
total workload constraints on f.. 

I 

In this appendix we discuss these items. The discussion is pertinent to a pilot control 
task, although the concepts can be applied to a pilot monitoring situation as well. 

B.l INCLUDING ATTENTIONAL CONSTRAINTS 


In the optimal control model, a fractional allocation of attention fj to the 
informational variable y. modifies the "observation" noise Vy. (t) associated with that 


variable (Reference 7). Thus, the noise covariance associated with y. is 

P? 2 
V = JL 

yi f. ' 

I 

where CT. = cr./N(a.) 

a. = RMS value of y. 

N(a.) = describing function gain of threshold a. 

a. 

m 

= erfc 


(B-1) 

(B-2) 



(B-3) 
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The noise/signal ratio p? is the "full attention" noise ratio, which typically is O.Oltr or 

“20dB. Note that Equation (B-1) represents an implicit relationship for the actual noise 

2 

variance Vy. since a. is itself a function of Vy, For a given Pj/fj/ ^he requisite Vy is 


solved for via on iterative procedure. The quantity P?/fj ~ P- is called the modified 
noise/signal ratio. 


The method that is used to determine how a pilot allocates attention among the 
various yj is to minimize the optimal control cost with respect to fj subject to constraints. 
This step will require an iterative process to arrive at the optimal f., where successive 
iterates'fj* result in lower values for the cost functional. The cost functional that is used 
to determine the f. is 

J* = min J(u) (B-4) 

u 

where J(u) is the basic cost functional in the optimal control pilot model. 


B.2 EXPRESSION FOR OPTIMUM QUADRATIC COST, J* 

In lieu of attempting to minimize the entire expression for J* numerically with 
respect to f., it is more efficient to isolate those terms in J* that are affected by f. . 
Since changes in f. are reflected as changes in observation noise Vy, we first obtain 
an expression for J* that shows the Vy dependence. 

The cost functional J(u), in the steady-state is given by 

J(u) = ECx'Qx + u'Q|.G} (B-5) 

where Q is defined in Equation (A -7) and x is the augmented vector x = col [x, u]. 
Equation (B-5) may be rewritten as 

J(u) = tr[QX + Q^E[uu'}] (B-6) 
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where X = E{x (t)x'(t)}. Assuming optimal pilot control, X is given by Equation (A-15). 
The control rate, u, may be approximated by the pilot's own estimate of u (as the actual 
u is modeled to contain a white motor noise). Thus 


^ = -Tn 

= -T-JO - T-^' Lx - T;;' 

' -h * - ■'n “e 


(B-7) 


where u^ is the error(u - u). Thus, since x and u^ are Uncorrelated for optimal linear 
estimation 


E{uu'} = L^Xlj + 


TN'](Ef ^ Ep) 


'n 


(B-8) 


A 

where E^, Ep and X are the filtering error, prediction error and (augmented) state estimate 
covariance matrices respectively. They are given by Equations (A-16) and (A-17). 

Substituting Equations (B~8) and (A-16) into Equation (B-6) yields with 
optimal control. 


where 


0 

+ tr|(Q + L' Q^L^)X} 


Aa 


We' 




Q, = Q + 


0 


(ENQr''TN>'' 


(B-9) 


(B-10) 


The observation noise affects only the error covariance matrix S which appears in the 
first and third terms and J 2 in J*. To show this dependence more clearly, we rearrange 

A 

J 2 » Substituting for X and using the cyclic property of the trace gives 
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Jj = (ry” (S + L'l Q, L,)e^‘'dcr . V'’ Cle*"" (B-'O) 

0 

But the integral term is identified as P, the solution of the Riccati Equation (A-7), thus 
using Equation (A-12), 


J3 = t,p'"Pe^^ (lA’+AZ + W)] 


(B-ll) 


which may be combined with and rewritten as 

+ J 3 = tr [e^''^Pe^'^W] + tr[e^''"(Q^ + PA + A'P)e^'" S] (B-12) 

The first term in the above expression is independent of and may be combined with 
J 2 which is also independent of V^. Thus, the only term of interest that remains is 
called the "scanning cost": 

Jq = part of J* dependent on 

= tr[e^''^(Q| + PK + X'P)e^'^i:] (B-13) 

This expression can be simplified by noting that 


+ pX + A'P = 


'’l2®r‘''’l2 I 0 
1 

0 

0 I -P22^r^*’22 

ffNQr''TN>’'. 


butQ”^P '^2 ~ ~ P22 

Jo = 

where L = "Equivalent" gains 




(B-14) 


(B-15) 
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and T, satisfies the variance equation 


0 = AS + 2 A' + W - SC V"^ C2 


(B-16) 


B.3 


GRADIENT EXPRESSIONS 


Minimizing Jq (and hence J*) with respect to f. represents a difficult nonlinear 
optimization problem. The difficulty is two-fold. First, fj affects in an implicit 
manner; second, affects Jq through the Riccati solution 2 . As it is unlikely that a 
closed-form solution for the optimum f. can be found, the optimization process will be 
carried out numerically via some form of gradient algorithm. 

In order that the numerical process be reasonably efficient, it is. desired to 
obtain closed-form expressions for the gradients - - ^ ■ or . Thus, the Jirne- 


^f: 




consuming process of numerically evoluoting these derivatives can be avoided. 
We wish to obtain the gradient vector 


9f 


_ a 


af 0 


0 - 17 ) 


where the i^^ element is ^ ^ . Consider a small change in f. . This will result primarily 

of* ^ * 

_ Po o. 2 


O'. ) but will also cause changes in the other 


in a change in V (since V - 
noise variances V , j- since or. will be modified* Thus, 

y* I 

af. av, * af. av, at. 

I I I 2 I 


= < 


av 


av 


^av \ aj 


af. 


JL> = 


0 


and so the gradient vector becomes 


,af. / av 

‘ I / y 


(B-18) 


3f 


^•^0 = r .^‘^0 


af 


av 


(B-19) 
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where is the gradient vector of Jq with respect to Vy., i = 1,2, . .m, and F 

O Vy I 

is a "transformation" matrix where 


SVy. 

(r)::= ■ 


■I af. 


(B-20) 


The gradient vector can be evaluated using a technique of Kleinman 

dVy 

(Reference 22) for derivatives of trace functionals. From Equation (B-16), the 


change in Z due to a change in to + 6V is given by 
6SA' +A62 +2C' V"\6V)V"’C2 =0 


(B-21) 


where we note (V + 6V) ^ = V ^Tl-6V*V~^ + (6V • V . . .1. The matrix 

y yL y y j 

A = A ~ZC V C is the closed-loop system matrix for the Kalman filter and has eigenvalues 

'^-•1 

with negative real parts. The term G =2C‘ is the Kalman filter optimal gain matrix. 
The solution to Equation (B-21) is given by 


OC 


6 Z - / e 

0 


Act 


GfiVG'e" ‘^dff 


Substituting into 6 Jq = tr (L 6 2L ') gives 


(B-22) 


S Jr 


av 


= diag (G'M G) 


(B-23) 


w 


here 


M 




e e 


(B-24) 


and is easily computed by solving the linear (Lyapunov) equation 


AM + MA + L' L =0 
e e 


(B-25) 
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P* A 2 

The elements of the matrix Fmust be computed. Since'V = — — O' . # 

n ^ ' 


we obtain 


bV V 
)1= -ZL 


^f. f. 


a. 


Sf' < 5 h 
«- I 


V V 

= T..=j[Ls^y\ 
■I f. . 

I 


(B-26) 


where is the used sensitivity coefficient, defined as 

y 

gX _ 5x/x 

y hy/y 


(B-27) 


If the terms happen to be 1 , then F is approximately a diagonal matrix 


V n 

r = r.. = YL. = - ^i 


II 


f. 


f. 


(B-28) 


I I 

I ' 

To determine more precisely the matrix F, rewrite Equation (B-26) as 


Vy. . 2 

Sf ^' = sf i - 6 .. 

'I 


I / I 1 m 


(B-29) 


"'vi 


so fhat if we can obtain Sr , the matrix Fis readily obtained. Since a change in f. 

i . 

induces a change in all of the , we. can write 


* 2 
<s . 


m 


2 

of 


Sf Sy • 

' k=i k . i 


(B-30) 


using the chain rule for sensitivity coefficients. Defining the m x m matrices 


Q = 


V 

/: _ 

Sf. ■ '’m 


(B-31) 
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and 


R = S, 


r.. 

'I 


(B-32) 


we obtain, substituting Equation (B-30) into Equation (B-29), 
Q = RQ - 1 


or 


Q = -0 - R) 


-1 


(B-33) 


The matrix R is next obtained since the terms r.j involve changes in the additive observa- 
tion noises and not in the multiplicative ratios. 

.2 

B.4 COMPUTING THE SENSITIVITY MATRIX, S^' 

yj 

The main computations required above are in obtaining the matrix R. Below, 
the method used to obtain R is described. We first consider the case of zero observational 
thresholds for simplicity. (This case is pertinent to informational level studies or labora- 
tory situations with high resolution displays*) Modifications for the case when visual 
thresholds a. > 0, e.g., when considering attentlonal allocation with realistic displays, 
are then presented. 

The augmented state covariance matrix is given by 


_ At„ A,t 


X ^ e" 'Ee""'' + / 


f 




CO 


+ f (EX’ + Xe + W)e^' da 

0 


(B-34) 


B - 8 


AEROSPACE SYSTEMSi INC. • ON6 vine BROOK PARK • Burlington. Massachusetts oieoa • (sit) 27a-7si7 


and the covariance of the outputs is 


Y= E {y(t)y'(t)[ = C X C 


Consequently, the variance O', is given by 


<J. = cl X c. = tr j X c. c! } 

I I I ( I I ) 


where cl is the i-th row of the £ matrix. Substituting Equation (B-34) into Equation 
(B-35) yields, after some manipulation, 

O'? = tr / fc. cl + A'Z. Z. a 1 e^^2 e^ + terms inde- 
' H ' ' ' J / pendent of Vy 

where 


OO^ 


- f A'o- , Act , 

- / e c. cl e do' 

J II 


and satisfies the linear equation 


A'Z. ■+■ Z.A + c.cl= 0 

I I II 


Since A - A - BL, the bracketed term in Equation (B-36) may be rewritten as 


H. =Z. BL + L'B'Z. 

I I I 


Finally, since 



G6V G' e^'^^do 


we obtain 




W ^ A 

/^, f A'a A'^ Lj Xt Act, 
|G / e e H. e e dcr 

' 0 
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so that 


r..=!ri 

i| 2 

C7. 


2 

daf 

I 


= 5y 


The above equations form the basis for the following algorithm to compute R, 
one row at a time. 











To reduce the computational requirements in computing R, the n’X r matrices 
(r = number of control inputs) 




I 1 ^ ^ 


are precomputed and stored. Note that these are dependent only on the optimal feed- 
back gains, and not on any noise statistics. Further reductions in computation are 
possible by using the linear equation algorithm of Bartels and Stewart (Reference 23) 
that permits efficient multiple solutions of 

XA + A'X = -Q 

with different right hand sides. However, even without this qigorithip, the computa- 
tional requirements for R are less than those needed for one Riccati equation solution, 
and must be termed slight. 

The above developments assumed a. = 0, i.e., no visual and/or indifference 
thresholds, so that the noise covariance associated with output y. is 


w _ Pi 2 
V - — a. 

y\ f. ‘ 


(B-41) 


When a. > 0, the effective noise covariance is determined from 

I 


V - Pi -2 

V - a . 

I 


I 


f. 


(B-42) 


where 


a = 


a/erfc / = a -J f 

[oJ-2 ) J 


and b = a/a f/2 . 


oo 9 

2 f -y- , 

e dx 
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In fhe non-zero threshold case, the matrix R in Equation (B-32) is given by 


‘"i V ho. 

(R):, = Sy ' — 

'I >. .2 3V 

' I rr \i 


(B-43) 


To relate this to the results of Equation (B-40), we write (since a. does not depend on 


^2 '^2 2 
Soj Saj doj a. 

hW 9a? 9V ha. ho. hV 


9a. doj 


or 


' (R)|i = S°' 
M 


'V 


(B-44) 


So that it is only necessary to obtain the m scalar quantities ,5a ./ba. to include the 

2 

threshold effects, since the sensitivities of. a. have been deternriined from the zero 
threshold Equation (B-40). From Equation (B-42) 


A 

5(J: 


5ct. 


1 - 


2b 


-b 


2 1 


|_ erfc (b) 


erfc (b) 


1 -1 


or 


^ F 

S = 1 - 


-b2 
2 b.e ‘^i 


*/jr erfc (b.) 


(B-45) 


where b. - 


a. 


o.j2 


. Thus, multiplying the i-th row of the zero-threshold R by 


Equation (B-45) gives the required sensitivity matrix for the situation a. > 0. 
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The above development assumed multiplicative observation noises. While this 

is true in general for manned vehicle systems, there are situations in which one might 

consider purely additive noises Vw (t) with covariances Vy. . Both additive and multi- 

/j 'i 

plicative cases can fit within the above framework as follows. 

The effective obser^^ation noise in the pilot model is 
o 

(B-46) 


_ ^ I -L 2 




where p. = noise/signal ratio if noises are multiplicative or p. = V? additive noise 
variance if observation noises are purely additive. In addition. 


a. 


a. = Verfc(b. ) ; b. = ^ 

' ' ‘ ojl 


(B-47) 


where x = ct. if noises are multiplicative and x = 1 if they are additive. To obtain the 

I r . 

gradient of J with respect to f, it was shown necessary to compute the matrix R: 

' ' a' 1 


R.. = S 

i| a 


(B-48) 


'I 


a. 


The computation of S. is identical for either additive or multiplicative noises. 
However, the matrix S will be different for these two cases. For the multi- 
plicative case. 


cr. 


1 _ 


au -W 
2b. e ' 




= 1 - 


(B-49) 


Jn erfc (b.) 


while for the additive case (x - 1 in Equation (B-47)) 

2 


CT. 

s ' = - 


OU -b:' 
2 b. e I 


(B-50) 


i erfc (b.) 
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Thus, it becomes relatively simple to compute dJ/bf for either type of noise cose by select- 
ing the proper form of either Equation (B-49)or Equation (B-50). Note that if a. =0, then 


s' = 1 0/ additive case 
^i ( 1, multiplicative case 


B.5 GRADIENT PROJECTIONS AND OPTIMIZATION 

SJq 

The above algorithm is used to obtain g, = — - which is the unconstrained 

t 

gradient vector. However, the attentional allocations f. are not free but are con- 
strained by ‘ 


m 


2 f. - f* - total attention and f . > 0 
i=l ' 


(B-51) 


The constraint. Equation (B-51), describes a portion of a hyperplane, 

< c, f> = f* (B-52) 

with c = col ll, 1,...,11. Thus, in order to determine the feasible direction for 
cost reduction, it is necessary to find the projection of g^ on the hyperplane. Equation 
(B-52). This is given by 


or. 


p. <9f.c> 

9f ■ — ' 

< C, C> • 


9f 9f 


(B-53) 



(B-54) 


In other words, g^ is obtained by subtracting the average of the g^, from 
each element of the vector. The angle between g^^ and is 
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(B-55) 


< 9r/ 9f > 

COS0 = — L__ 

Bg^D • llg^ll 

The gradient pro{ection. Equation (B-54) is appropriate at the informational 
level, i.e., where no distinction is made that two "outputs" come from one display 
indicator. In the display format studies, the fact that a pilot obtains both position and 
velocity observations from a single display indicator provides yet another constraint on 
the f.. Assuming outputs are ordered in position'Tote pairs, 

f| — f2/ fg “ f^/ • • » (B“56) 

(assuming /2 “ etc.). For this case, secondary projections of g^ are necessary on 
the planes “ ^2 ~ 
with (g? + gf )/2, 

display level studies are easily treated insofar as gradient projections are involved. 

In order to develop a reduced (projected) gradient optimization scheme, 
assume we are at iteration n, with attention vector f*^. A small change Af^ such that 
fT + Af'^ still satisfies the constraints of Equation (B-51) will cause a corresponding smalj 
change in Jq. Thus, at iteration n, 

jS+<9^Af"> (B-57) 

If Af*^ is selected as 

2 

Af" = -eJ^gJ/lig^'l! e< 1 (B-58) 

then, 

Jj^'- (1 - . ) (B-59) 
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0, f 3 “ 0/ etc. This is easily done by replacing and gP^ 

etc. Thus, the various constraints imposed by information and/or 


Thus, each successive iteraHon will result in a lower cost (to first order terms only) 
ofl00e%. 


Equation (B-59) serves as the basis for a gradient optimization scheme. We 


e - ' P cos 0 


(B-60) 


and pick P < 1/2 and sufficiently small such that f*^^^ = f*^ + Af*^ satisfies > 0 and 
also < Jq • Convergence occurs when Is arbitrarily close to Jq . If > 
Jq, a smaller step is taken by reducing /3, Note that since Af*^ is in the direction of 
g^,-fhe resulting f”^^ must necessarily continue to satisfy the constraints imposed on f*^, 

. When using the above gradient projection algorithm to optimize attentional 
allocation, it is necessary To modify . the scheme whenever 


g^ > 0 and c > f . > 0 


(B-61) 


In such a situation the mathematics do not want to.place any attention on Instrument I . 
Thus, one cannot move in the direction of the gradient since f. would go negative. The 
modification that is obviously called’ for is to set 

9f. = 0 

I 

whenever the above condition. Equation (B-61), is encountered. Thus, f. remains 
fixed for the iteration and the re-projected gradient is modified to guarantee that 

® 

I 

I 

as required. 
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APPENDIX G 


LINEARIZED ROTORCRAFT EQUATIONS OF MOTION 


In developing a mathematical model for use in the analysis of an automatic 
approach and landing system for a tandem rotor helicopter, linearized equations not 
normally used in aircraft stability analysis are required. 

C.l DERIVATION OF THE NONLINEAR EQUATIONS OF MOTION 

Before developing the linearized equations of motion for analysis purposes, 
the nonlinear equations are developed. 


C.l.T ROTARY MOTION 

The equations are written in the conventional body axis system, with the 
origin at the helicopter's center of mass. The x-qxis (roll) is directed forward, the 
y-axis (pitch) is directed out the starboard side, and the z-axis (yaw) is directed 
"downward" to complete the orthogonal set. Assume that the x-z plane is one of 
symmetry and define the inertia matrix as 



0 

I 

yy 

b 


J 

xz 

0 

I 

zz 


where 



J = - 
xz 


J* xy dm 
J' zy dm 
J " xz dm 



due to x-z symmetry 
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The components of the angular velocity vector In the body axis system are given by 



where P, Q, R are the components of the angular velocity resolved along the body axes 
X, y, and z, respectively. Similarly, let the applied torque be denoted as 



With the assumption that the angular momentum of any rotating machinery onboard, 
such as engines and, in helicopters, the rotor, is negligible, the equations in inertial 
space are written simply 

T = A(Iu,).' ■ (C-1) 

dt . 

On the assumption that the Earth is inertial, then Equation (C-1) is written in the 
body oxis system as 

T = l-iiiL +«>xla> (C-2) 

dt 

where the prime denotes “as seen by an observer fixed in the body axes," and therefore 
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where the dot denotes total time derivative. Substituting the appropriate definitions 
into Equation (C-2) results in the equation set 


:l “ 


ri p + J Ri 

XX xz 


rj PQ + 1 OR - 1 OR 1 

xz zz yy 

M 


1 Q 
yy 

+ 

1 RP + J (R^-P^)-l RP 

XX xz^ zz 

-N. 


J P + 1 R. 

xz zz 


J PQ - 1 PQ - J RQ , 

yy XX xz 


These equations define the rotary motion of the helicopter. 


(C-3) 


C.1.2 TRANS LA TORY MOTION 

Assume that the mass of the helicopter is constant and denoted m. Let the 
force vector be denoted . 


■'A 

F = 

where X, Y, and Z lie along the body axes. Note that in the conventional body axis 
system the total velocity 

U' 

V 
W. 



rx 

Y 

-Z_ 


is resolved along the body axes. 


The equation in inertial space may then be written 


F = 1 (mVJ. 
dt ' 


(C-4) 


On the assumption that the Earth is inertial, then Equation (C-4)is written in the body 
axis system as 
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d’Vy , . 

F = m +u)xmVT . (C-5) 

dt ' 

Making appropriate substitutions into Equation (G- 5) gives the equation set 

X "mU" "mQW - mRV 

Y = mV + mRU - mPW 
Z . _mW_ _mPV - mQU 

These equations define the translator/ motion of the aircraft. 

CJ .3 ^EULER AI^GlV TRANSFORMATIONS 

For an aircraft in the conventional body axis system, the attitude can be 
described with respect to a set of axes fixed in the Earth. To do this we define an Euler 
angle set denoted Y, 6, §. These three angles are the azimuth change, elevation 
change, and roll required to arrive at the aircraft attitude from the inertial axes. 

They must be taken in the given order. If a vector is denoted in the inertial coor- 
dinates as C. and viewed from a coordinate system which has been slewed through Y, 
we get 

c, = T^(f)C, (C-7) 

where 

cos Y sin Y 0 

A 

T (Y) = -sin Y cos Y 0 

0 0 1 


C -4 

aerospace SVSTEIVISi II\IC« ONE VINE BROOK PARK • BURLINGTON. MASSACHUSETTS 01803 • 





(017) 272-7317 


If fhe observer is fhen elevated through 9, he sees 
Cj = Ty(e)C^ = 

where 



‘cos 0 

0 

-sin 0 

A 

T (0) = 

0 

1 

0 

y 

_sin 9 

0 

cos 0 . 


(C-8) 


Finally/ if the observer is rolled through he sees from aircraft body axes 

^b = = y«)Ty(9)T^(Y)C. (C-9) 

where 

0 

sin f 
cos $- 


A 

T (9 = 

X 


1 0 
0 cos $ 
LO -sin 9 


C.1.4 GRAVITY FORCES 

With Equation (C-9) the gravity force may be written in body axes as 



= T^(§)Ty(e)T^(Y) 


" 0 ■ 


'-sin 0 

0 

= mg 

sin f cos 9 

_mg- 


cos f cos 0 _ 


(C-10) 


C.1 .5 EULER ANGLE RATE EQUATIONS 

The body axis rates may be written as functions of the Euler angle rates: 
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■p" 



rii 


"O" 


foi 


Q 

.R. 

= T (§) 
x' ' 



+ \(S) 

• 

e 

A 

+ T^(’f) 

— 1' 
O 




■] 0 -sin 0 


V 


0 cos 0 sin $ cos Q 


• 

e 


_0 -sin 0 cos ^ cos 0_ 


_Y. 


C.].6 INERTIAL VELOCITY EQUATIONS 

If the velocity in inertial coordinates is denoted 





then we may write the body axis velocities as functions of the inertial velocities 


"U" 


fv 1 

X 

V 

= T^(S)Ty(e)T^(Y) 

''y 

_w. 

♦ 

V _ 
z 


V cos 0 cos Y + V cos 0 sin Y - V sin 0 j 
X y z I 

V^(sin $ sin 0 cos Y - cos $ sin Y) I 

V^(cos $ sin 0 cos Y + sin f sin Y) | 

1 0 

V (sin $ sin 0 sin Y + cos$ cos Y) + V sin $ cos 0 
y z 

V (cos $ sin 0 sin Y - sin § cos Y) + cos 9 cos $ 


(C-12) 
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C.2 LINEARIZING OF EQUATIONS USING DELTA PERTURBATIONS 

The purpose of this linear model is to onalyze.an automatic approach and 
landing system. As such, the system is tied to the Approach Navigation Frame (ANF). 
This is an Earth-fixed coordinate system with the origin at the desired touchdown point, 
the X axis along the runway, and the Z axis down along the local vertical. The system 
is considered inertial. The variables to be commanded will be in the ANF such as Vx, 
Vy, and Vz and Y, 6, f . The control system will therefore be feeding back th(ese 
quantities from an inertial platform. The linearized equations are then desired jn terms 
of these variables. The derivation of these equations follows, 

C.2.1 BODY TRANSLATORY EQUATIONS 

Rewrite Equation (C-6) with the forces being divided intp aerodynamic forces 
and the gravity forces of Equation (C-10) 

X. 

U = RV - QW + — - g sin e 
m 

Y^ 

V = PW - RU + + g sin f cos 9 (C-13) 

m 

^A 

W - QU - PV + + g cos § cos 0 

m 

Assume the following perturbations 


u = 

Uq+U 

P = p 

0 = 

CD 

+ 

o 

CD 


V = 

V 

o 

II 

$ = 


I (C-14) 

w = 

Wq + w 

R = r 

Y = 

Y 
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Note that the perturbations are simply added to the Euler angles. If steady 
flight is assumed, the steady flight equations are obtained from Equation (C-13) 


^AO 

0 = g sin 


m 


0 = !^ 


(C-15) 


m 


^AO 

0 = + g cos 0Q 

m 


Next note that for small 0 


sin( 0 Q + e) = sin 0 q + (cos 6 q) 8 


(C-16) 


cos(0Q + e) = cos 0Q - (sin 0 q) 0 

Substituting the definitions of Equation (C-14) into Equation (C-13), using Equation 
(C-16) and subtracting off steady flight. Equation (C-15), yields the following perturba- 
tion equations 


u WQq - (g cos 9 q)0 + 


6X, 


m 


V = WqP - UQr + (g cos 0 q)§ + 


6Y. 


(C-17) 


m 


w = UQq - (g sin 0 q)0 +■ 


6Z, 


m 
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C.2.2 BODY RATE EQUATIONS 


Substituting the perturbation definitions of Equation (C~14) into Equation (C~3) 
and dropping second-order terms gives 


• _ _ ~’xz • ^ AL 
I I 

XX XX 


aM 


(C-18) 


yy 


• _ _ xz • ^ AN 


C.2.3 EULER ANGLE EQUATIONS 

The relation between Euler angle perturbed rates and body perturbed rates is 
desired. This is obtained by inverting Equation (C-lT) and applying the perturbation 
definitions to give 


0 = q 

§ = p + (tan 0 q) r 


(C-19) 




cos 




C.2.4 EQUATION SET WITH AERODYNAMIC PARTIAL DERIVATIVES 

The motion of the aircraft is described by Equations (.C-T7)/(C-18), and (C-19). 

The force and moment perturbations are, in general, functions of the motion and can 
be written to first order as a linear function of the motion variables. Motions in the 
longitudinal plane are assumed to separate from the lateral-directional to give the 
total equation set which follows. 
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X. X. 

u = -WQq “ (g cos 0 q) 6 + — u - w + — Sq + — + — 2c A 6 


m m m m 


m 


Y Y Y Y Y 

= WqP - (jQf + (g cos 6 q)§ + —V + —£p + — r + — ^ A6g + — ^ 

m m m m m 


w 


Z Z Z Z. ^6 
= UQq - (g sin 0 q) 9 + — u + — w + — 3. q + — ^^6^+ — 2 A 6^ 


m m m m 


L L L '-6 *-6 

p= -iE-r+^v + -P-p+-Jlr+_3A6a+~lA6, 

^xx *xx ^xx 'xx 'xx ^xx 


. M M M ^6 ^6 

q = — iiu + — ^ w + — 9- q + 9 A 6 + 9 A 6 


yy 'yy 'yy ’yy 


I 


yy 


N 


J N N N N ’’'6 

r*= ^p+!j^v + -Pp+-^r+^^^^ 

I I 1 I I i 

zz zz zz zz zz zz 


0 = q 

§ = p + (tan 0Q)r 


Y = 


cos 0, 


> (C-20) 


In the above expressions 6^, 6^, 6^, 6^ are the four control (displacements 
corresponding to differential collective, gang collective, roll cyclic and yaw cyclic, 
respectively. 
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C.2.5 EQUATION SET WITH BODY RATES ELIMINATED 

Equation (C-20) may be reduced in variables by eliminating the body rates 
with the last three equations. The result is 


u = 


X X 
u . w 
— u + — w 

m m 




(g cos 9 q)0 + A 6 


m 


— V + ^Wq + $ + (g cos 9 q)« - ^sin 0 q |Wq + + cos 9 q ^Uq - 


+ _A 8 
m 


w 


Z Z / : Z \ . Z^ 

= -JiL u + — ^ w + (Uq + — 3. J 9 - (g sin 9 q)9 + — A 
m m \ m / m 


L ^ L /J 

^_^v + _p_^ -/JE 

I I \l 

XX XX \ XX 


cos 0Q - sin 0 


' \ XX XX ' XX 


.. M M . M 

9 = — 3. u + AW + — 3 A 0 + A6 


I I 

yy yy 


I I 

yy yy 


Y = 


I cos 0r» + J s'fi 6ft 
zz 0 xz 0 


(N V + J $ + N $ - (N sin 0« - N cos 9^)Y + N . A6) 
' V xz D ' p 0 r U' 6 ' 


where ( A8 indicates the summation of partial derivatives and actuator deflections, 

C.2.6 REDUCED EQUATION SET IN VECTOR-MATRIX NOTATION 

Putting Equation (C-21) in vector-matrix notation results in the following 
longitudinal equations: 
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C.2.7 INERTIAL VELOCITY EQUATIONS 

If Equation (C-12) Is Hh^arlzed, the result Is given 

u = cos 0Q - AV^ sin 9 q - Wq 9 \ 

V = AVy + Wq f - (Uq cos 9q + Wq sin 0 q) 'll I (C-24) 

w = sin 9q + AV^ cos 9 q +, Uq0 / 

With these equations, the body axis velocities can be eliminated from Equation (C-21) 

and the result is given 
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( XX \ /x X \ 

_ii + -Ji; tan 0 Q) LV + (tan 9o)^V - I _ii tan Sq - ~)aV 
mm / ' . ni m / 

cos0«\m \m m / m / 


Y Y . 

L V = — AV + — £■ A0 + 1 — - Wa + g cos 0 


y y 

' m ' m 


^ ^ — Wq + g cos 0 

- — cos 0Q J i - -1 |uq cos'Oq + Wq sin 0 q^ i]( + — 
m / m \ / m 


A (l-on ^ { ~ ^ ~ ®oV^x " f “ ®0 " "^) 

\mm / \m m/ 


AV 


cos 0^ \ m \ m 


+ , 


( M M V 

^ cos 9 q + ~ sin 00 ^ "( 

yy yy ' ' yy 


z \ 

Wq - — Uq + g sin 0o)0 + — 6 
m 




yy 

M . /M 

+ -^0 

'yy X yy 


M 

L 

yy 


M 


sin 0Q ^ cos 0 qJaV^ 

yy 




10 +-^ 6 

I 

yy 


i(C- 


(C-25) 


95 = 


-JLaV +-E- 

I y I 


XX 


/M M \ 

■ r""-r“» 

\ yy yy ' 

L . L /J /L L 

0 + Wq 0 +j cos 0Q + sin 0Q J ij( - I sin 0 q cos 0 q J'lr 

^XX ^XX ' ' ^XX ^XX f 


- 7^ (^0 ®0 P 

XX ' f XX 


'1' = - [Ny^Vy + j^^0 + Np0 + NyWQ0 - (Np sin cos 0q) 

- Ny(L)Q cos 0Q + Wq sin 0 q)+ + N g6] 
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C.2.8 


INERTIAL VELOCITY EQUATIONS IN VECTOR-MATRIX NOTATION 


Equation (C-25) may be put into vector-matrix notation giving the 
folloNving longitudinal equations: 



(C-26) 


and the following lateral equations: 



(C-27) 


ti ' 

where V^q ~L)q cos 0q + Wq sin 0 q, 
differential operator. 


t . L ^ 

B = 1^^ cos 0Q + sin 0Q and S = — the Laplace 

dt 


zz 
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C.2.9 STATE-VECTOR EQUATIONS 


The equations above are not in a form suitable for the state-vector formulation 

• • • • 

because the equations containing and ^ , (t> are each simultaneous sets of 

• • • • • - . : ' ^ . 

equations. To solve for and ^V^, we note that these two equations in Equation:-. 
(C-27) can be written 




e, 9, 4 6 

f2(4V^/4V^, 9, 9, 46y 


(C-28) 


where 4 6 are the actuator deviations from trim settings. 
Inverting these equations leads to 



(C-29) 


where the matrix D is given by 



(C-30) 
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id 


/x X 

* u . w • 

— cos + 5**^ Qf] 


L\ 


m 

k 

u 

m 


nri 

m 


- — sin 0Q - — cos 0Q 
, m m 


m m 


ro ) 

— WO^ — 

I m m 


— • VVq - ■— Uq + 9 sin e£ 
m m 


X 

Ji 


Similarly, the equations containing ijl and *0 can be written 
-J. 


xz 


XX 


cos 0Q - sin 0Q 


-J 


xz 




zz 


cos 0Q + sin 0Q 
*zz 


cp 


/ 9, 9, Ijf, if, A6) 

^q2(AVy, 9, 9, i|f, ijt, A6) 


Inverting these equations to solve for 9 and yields 




1 - 


xz 


'xx'zz 


1 + tan 0, 


ZZ 


xz 


XX 


+ tan e. 


xz 


] 


'zz 


cos 0< 


\’'2l 


or 



♦ 

A6 


where the matrix E is given by 
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1 + — tan 9/ 


+ tan 9f 


(C-35) 




— I 


1 1 

1 1 'zz 


XX 

1 


L 

L 

L 

V 

V 

P 


1 0 
XX 

1 

XX 

N 

N 

N 

V 


p 

1 

r ® 

1 

zz 

ZZ 

zz 


cos 0/ 


V 

1 

cos 

®0 

- p 
1 

sin 0Q 

6 

'xx 



XX 


'xx 

N 



N 


N, 

V 

1 

cos 

®0 

- P 
1 

sin 0Q 

6 

1 

zz 



z? 


zz 


C.3 WIND GUST DISTURBANCES 

C.3.1 GUST MODEL 

The wind gust model is derived from the Press-Meadows analytical representa- 
tion for the power spectral density (PSD) of random turbulence (Reference 24): 


2 L ' ^ 
9(w) — <y ■ — — ^ 


,^^o 1+(1/VJ^U)^ 


(C-37) 


where ?((b) = power spectral density, ft /rad/sec 

2 2 2 
a = mean-square gust velocity, ft /sec 
9 

L = scale of turbulence, ft 
= equivalent airspeed, ft/sec 
u) = frequency, rad/sec 

To obtain a random wind model, w , which satisfies the Press-Meadows power 

9 

spectral density^ a wind filter can be used to shape a unity rms white noise input. Since 
the PSD is 


$(m) = Iwg(ju))l^ = Wg(ji«)w* (iu)) 


(C-38) 
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where w is the complex conjugate of w , spectral factorization can be used to 
obtain w (w). Replacing jw by the Laplace transfer function equivalent for the wind 
filter: 


w (s) (s+Qf) 

—M — = a Y 

n 3 (s + P)2 

where 

Y = yjzyyi 

O' = L 

V/L 


(C-39) 


and n is the driving white noise. 

Using partial fraction expansions. Equation (C-39) can be written as 


W (s) 


O' - p 1 

+■ 

(s + p)2 (s + P) 


(C-40) 


A block diagram of the wind gust model resulting from Equation (C-40) is shown in 
Figure C-1 . The corresponding model in the time domain is obtained by defining 
the two wind gust variobles 



Figure C-1 . Block Diagram of Press-Meadows Wind Gust Model. 
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X 


w 


w 

w 


gi 

9^. 


The linear equations describing the wind gust model become: 



X + B u 
w w 


E n 
w 


where . 


“1 B=n ,Jv] 

b-p -pJ LoJ LoJ 


The total gust velocity is 


(C-41) 


(C-42) 


(C-43) 


W = W , + W /K 

g gi g2 

and the wind gust rate is 


(C-44) 


(C-45) 


The previous development is applicable to gusts along any of the three aircraft 
i* 

axes' . However fhe appropriate values must be selected for the model parameters a , 

g 

L and V^. The magnitudes of the turbulent velocities are represented by their rms values 

0,0 . and a (downwind, crosswind, and vertical, respectively). The rms com- 

gu' gv' gw ' 

ponent of a is approximately 

y V 

a = 0.2W (C-46) 

gu 


^The Dryden model of atmospheric turbulence PSD is equivalent to Press-Meadows 
for turbulence normal to the airspeed; for turbulence parallel to the airspeed, the 
Dryden PSD is 2 ct 2(l/Vj/[] + 
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where W is fhe mean wind speed (Reference 25), The ratios of the rms speeds near 
the surface are about 


gu gv gw 


= 1.0: 0.7: 0.5 


(C-47) 


These increase with altitude so that all components are nearly equal at 500 ft (Reference 25). 
The turbulence scale L is approximately equal to 1,000 ft for altitudes above 1,000 ft; 
below 1,000 ft, the scale for vertical turbulence is approximately equal to the altitude, 
while the scales for the horizontal components decrease to about 500 ft qt the surface. 


3.2 AERODYNAMIC FORCES DUE TO WIND GUSTS 

The atmospheric turbulence provides disturbance inputs to the aircraft through 
the aerodynamic forces. Each of the wind gust components may induce important air- 
craft responses, but the vertical component w is primarily responsible for normal 

9 

accelerations. This component is generally considered the most important disturbance, 

and is the only one treated explicitly in the following discussion. However, the analysis 

of u and v follows the same general pattern. 

'99 

For most rigid-body analyses, it is sufficient to omit the small portion of the 
turbulence spectrum for which the gust velocity gradients are not adequately represented 
by equivalent rates of aircraft pitch and roll (Reference 24). The aerodynamic effect 
of the vertical gust is to modify the angle-of-attack and the pitch and roll rates. 

Effective angle of attack = a + 0/ =(w-w )/V 

g g o 


Effective pitch rate = q+ q -q+w VV 

9 9 ^ 

Effective roll rate = p+Pg = p- w 
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These modified values are used in calculating the aerodynamic forces and moments. 


e.g.f instead of (X^^m)w 

in the equations of motion, there will appear the term 

(X^^rn)(w-Wg). Of course 

, none of the inertia terms in the equations of motion 

are affected. 




The result of the wind input appears as additional terms on the right hand 

side of Equation (C-20). 
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(C-48) 


When these terms have been carried through the derivation, the end result is the addition 
of the following terms to the inertial frame equations: 
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where J - (1 ~ J 

xz 


AEROSPACE SYSTEMS, 



/'xx 'zz> 


-1 


The gust stability derivatives have been defined as 
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X. =x/v 

wg q o 

Y. = -Y /V 

wg p o 

Z = -Z 
wg . w 

'■w9 = V''o 

M = -M 

wg w 

M. = M /V* ' 

wg cf o 

N • = -M /V 

wg o 


(C-50) 


The gust rate coefficients in Equation (C-50) are Inversely proportional to the 
airspeed and become infinite as approaches zero. To avoid this difficulty in 
hovering conditions, an equivalent airspeed can be used in the wind calculations. From 
Equations (C-46) and (C-47), the rms vertical turbulence velocity is approximately 0.1 of the 
mean wind; therefore, for a specified turbulence level, the equivalent airspeed can be 
defined as 


V = 1 Oa 
o gw 


(C-51) 


C.3.3 LONGITUDINAL EQUATIONS WITH GUST INPUTS 

The state variables for the longitudinal equations must be augmented by the wind 
states in Equation (C-41): 


X - [w , w , Az, A6, AV , AV , q] 
9] 92 X z 


(C-52) 
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Using Equations (C-44), (C-45), (C-49), the wind inputs alter the vehicle equations 
as follows: 
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APPENDIX D 

CH-46 DYNAMICS WITH MODEL-FOLLOWING AUGMENTATION 

This appendix presents the derivation of the equations of motion for the 
CH-46 to be used with the optimal control model of the pilot. The basic equations of 
motion of the aircraft with the four actuator inputs are contained in Appendix C. 

The first task is to derive the equations of motion for the control stick inputs taking 
into account the attitude command system used on the CH-46C whose flight test data 
will be used to confirm the pilot model. The second task is to specify that the output 
equations of the controlled element are consistent with the displays implemented in 
the CH-46C. 

D.i ROTORCRAFT DYNAMICS WITH MODEL-FOLLOWING AUGMENTATION 

A block diagram of the LaRC Model -Following Control augmentation system 
for the CH-46 is shown in Figure D-1 . In this section we derive the effects of this 
type of augmentation system on the dynamics controlled by the pilot. The difficulty 
arises because of cross-coupling between the control actuator dedicated to provide the 
model response and the effect of that actuator on other state variables. For example, the 
differential collective may be used to satisfy the pitch response of the helicopter, 
but by so doing, it will have an influence in the translational equations. 

The basic unaugmented vehicle equations are written in the form 

x = Ax + B6+Fw (D-1) 

or 
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In Equation (D-2) the state vector has been partitioned such that X2 represents that 
portion of the state to be controlled according to some desired model response. The 
Inputs to this system are the actuator inputs ft which are dedicated to satisfying the model 
response. Disturbance inputs and other actuators affecting the equations of motion, but 
not being used to satisfy the model response, are included in the vector w. 

Assume that it is desired to have the partition of the state vector X2 follow a 
model equation given by 


.m _ .m 
X 2 A2^x^ 


. A m m 
^22 ^2 


□ m 
B2 u 


(R-3) 


Note that the derivative of the model response is given by linedr combinations of the 
model state itself x^ , feedback of the other partition of the state vector x.j (e.g,. position), 
and the control stick inputs u. 


If the control augmentation system is well designed, then it will force X2 ~)<2 
and )<2 ~ :<2 . In this case, we may subtract the second part of the state Equation (D-2) 
from the model Equation (D-3) and solve for the actuator activity which is dedicated to 
satisfy this response: 

62^ ^^^21 " ^21^ ^ ^^22 " ^22^ X2 + B2 u - F2W] (D-4) 


Substituting this value of the actuator activity into Equation (D-1) provides the follow- 
ing form for the state equations when the system follows the model: 

*2 


• *1 * - 

A 'A 

Mil Mi2 

—5 — r 

A ' A* 

M21 , M22 


"1 


U + 


w 


(D-5) 
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where 


A* = A + B B (A^ - A ) * A* = A*^ • A* = 

^11 '^r2 ^^21 2r' 2i ^2i ' ^22 ^22 


^12 ^ 12 '^® 1®2 ^^ 22 '^ 22 ^ 


(D-6) 


B, = B,B2 B^ B 2 =B-; F^ = Fi-B, B2-F2 j 

Note that the partition of the state vector X 2 has the desired response of the model, but 
also that the dynamics of the other partition of the state vector :jC| have been altered 
because of the effect of the dedicated actuator 6 on these equations. 

D.2 LONGITUDINAL EQUATIONS OF MOTION 

The equations of motion for the longitudinal axes are derived in Appendix C. 
We use as state variables in these equations, the Euler angles and body rates to be 
consistent with control augmentation schemes and guidance schemes employed. When 
we partition the longitudinal equations to account for the attitude command system in 
pitch, we have 


aVx / = 


(D-7) 


D - 4 


AEROSPACE SYSTEMS, INC. • ONE VINE BROOK PARK • BURLINGTON, MASSACHUSETTS 01S03 • (S17) S7S>7ai7 


where 
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The coefficients Dj. are given in Appendix C. 

The desired pitch response system may be vsrritten, in terms of the present 
notation, as 


= (0, 0, 0) 


li) = 2.0, 1 .43 

n 


Aoo = -2Co)^ 


C = .75 


(D-12) 


BX = K.w 
2 9 n 


Kg = .15, .10 rad/in 


Substituting these results into Equations (D-5) and (D-6) will provide the equations of 
motion in the longitudinal axes under the attitude command system. 

D.3 LATERAL DIRECTIONAL EQUATIONS OF MOTION 

In a manner similar to that for the longitudinal case, the lateral directional 
equations of motion can be written to account for the attitude command systems in 
roll and yaw. Using the formulation above, these equations become 


A ! A 
j ^12 


^21 j ^22 




C:) 


(D-13) 
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(4x4) 


where 


0 

0 

0 

0 


0 

0 

0 


0 

0 

0 


°42 °43 


1 

0 

0 

Yv 

m 


Yv 


3 - W + g cos 9 

.42 m o o 


U3 


- “(U cos e + W sin e ) 


Yv 


(D-14) 


o m 


‘12 


0 

1 

0 

Y 

m 


tane^ 

_T_ 

cos 9^3 
Y 

r 

m 


(D-15) 


^21 


'22 


15 


-25 


11 

^12 

E ,3 

21 

^22 

CO 

CN 

UJ 


16 


■26 


(D-16) 


b, = Y 
1 f) 

c 

^2 = \ 


(D-17) 
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^2 = 


■16,6 


-26,6 


■16,6 

• ir 


■26, 6 ^ 


(D-18) 


See Appendix C for the coefficients E”. 

For simplicity, we show the model response equations in yaw assuming that 
the heading hold mode is employed. In this case the model response equations become 


- 0 ) 


'21 


0 


-u) ^ 0 

o 


-2CU), 


'22 


where 




= 2 , 1.43 


(D-19) 


(D-20) 


C 

U) 


= .75 

= 2.0 




= .70 


2 

u) 

n 0 


62*6 


where 


2 

u) K 

o 


K = .15, .10 rod/in 

0 

= .35 rod/sec/in 


y 

6 . 


(D-21) 
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D.4 PHYSICAL AND AERODYNAMIC CHARACTERISTICS OF NASA/LaRC 

CH-46C HELICOPTER 

Table D-1 summarizes the pertinent physical characteristics of the vehicle, 
including control displacement limits. Tables D-2 through D-7 give the trim conditions, 
fuel flow and stability derivatives as functions of forward speed for three rates of 
descent. The data in tables D-2 through D-7 were calculated for a gross weight of 
13,400 lbs, normal c.m. position and sea level flight. In Tables D-4 through D-7 
(non-zero rate of descent), the data for forward speeds greater than 80 knots were 
extrapolated from 80 knot values by assuming constant offset from the corresponding 
zero rate of descent (Tables D-2 and D-3) data. 


Table D-1. Physical Characteristics of 
LaRC YHC-IA Helicopter. 


Parameter 

Value 

Units 

English Units 

Operating mass, m 

6078 

kg 

(13,400 lbs) 

Rolling moment of inertia, l^ 

12,474, 

L 2 

kg-m 

(9203 slug-ft^) 

Pitching moment of inertia, I 

yy 

102,898 

. 2 
kg-m 

(75,914 slug-ft^) 

Yawing moment of inertia, I 

zz 

97,238 

2 

kg-m 

(71,738 slug-ft^) 

Cross-pr<xluct of inertia, 

-9638 

kg-m^ 

(7144 slug-ft^) 

Reference area, 

341 

2 

m 

(3670 ft^) 

Reference chord, c 

0.4572 

m 

(1.5ft) 

Rotor radius, R 
' r 

7.Z67 

m 

(24.17 ft) 

Control travel limits 



' 

Collective stick, 6 

c 

0-32.512 

cm 

(0-12.8 in.) 

Longitudinal stick, 6^ 

+13.97 

cm 

(+5.5 in.) 

Lateral stick, 6 
. ' a 

+9.14 

cm 

(+3.6 in.) 

Pedal, 6 

+5.85 

cm 

(+^2.3 in.) 

r 
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Table D-2, YHC-IA Longitudinal Stability Derivatives (Rate of Descent - 0 ft/min). 
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Table D-3. YHC-IA Lateral Stability Derivatives (Rate of Descent = 0 ft/min). 


1''0;<U'AUD 

VET,OCl'i,’Y 

KTS 

0 

20 

40 

60 

80 ! 100 ], 120 

L. ! 

140 

0 

in 

.12983 

.08462 

.09191 

.13466 

.17931 

.29754 I ,4415^1 

.49785 

\ 

In 

-.17764 

-.04701 

-.08508 

-.33705 

-.61293 

-.75397 

• 

-1.04713 

-1.14123 . 


ft/scc^ 
Tl:/ sec 

-.02663 

-.05408 

-.07986 

-.11823 

-.16121 

-.18040 

-.19766 

-.23442 

Yp/m 

J^_/sGcf7 

xaif/ ^5C;jp 

-.76514 

-.95618 

-1.15981 

-1.24824 

-1,23341 

-1.02295 

-.76619 

-.29372 

Y^/»i 

it/sec® 
rnd/ soc 

-.12517 

-.18179 

- . 20082 

-.12813 

-.05819 

. 

.13380 

.12514 

,22125 

. 


f t/sec^ 

.99794 

.99700 

.97673 

.96426 

.96955 

.97123 

L._ . • 

1,01613 

1.10225 


r- 

ft/sec^ 

In 

.14652 

.13634 

.11652 

.10583 

.09849 

.11641 

,13883 

.15262 


rad/ sec^ 
iU/sGc _ 

-.00778 

-.01305 

jjjjHQQQ 

-.01720 

-.02442 

-.03188 

• 



-.04724 

-.05080 


rad/ sec® j 
rad/ sec- ' 

-.50730 

-.57483 


-.65271 

._1 

-.62549 



-.23476 


rad/sec ^ 
racf/^fi'c 

-.02297 

-.04571 

-.05943 

-. 02270 j 

.01172 

.09531 

,11659 

.19955 

^63/^xx 

rad/ sec^ 
in _ 

.46536 

.46558 

.45954 

.45575 

.45805 

.45820 

.47354 

wmm 

IBimi 

rad/ sec*® 
In 

-.12638 

- ;. 13036 

-.13422 

-.13544 

-.13910 

-.13256 

-.13248 

-.14433 

N /I 
V zz 

rad/sec'^ 
it/ sec 

.00013 

.D0037 

-.00123 

-.00290 

-.00415 

. .j 

-.00536 

-.00377 

-.00060 

V'xx 

rad/ sec^ 
xns]lji\\c^ 

-.01831 

-.02076 

-.02598 

-.03947 

-.05344 

-.06864 

-.07747 

-.07471 


rad/se^ 

racr/ircc 

-.05847 

-.05450 

-.04597 

-.05020 

-.05438 

1 

-.08078 

-.08642 

-.12617 

hn„. 

rad/ sec® 
in 

.03001 

.02923 

.02764 

,02663 

.02633 

.02696 

, 02857 j ,02970 


rad/ sec® 
Ld 

.17534 

.17578 

.17261 

,17032 

.17141 

.17146 

. 17933 ' .19536 
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Table D-4. YHC-IA Longitudinal Stability Derivatives (Rate of Descent ^ 1500 ft/min). 


' FORWARD 
VKI.OOTTV 

KT8 

0 

20 

, 90 

60 

80 

100 

3. 20’*' 

140 

“o ■ 

clcy 


8.13867 

6.99471 

5.34.340 

3.21352 

2.32445 

-1.37813 

-5.74252 

In 

.03060 

-1.69234 

-.56815 

i 11550 

.60565 

-.23371 

-.00860 

.17046 

6 

% 

i.n • 

3.89872 

3.25038 

1.25071 

.90347 

1.15881 

1.98741 

3.35947 

5.32989 

f 

Ib/lu- 

791.8A9 

711.302 

500.687 

470.262 

495.961 

568 , 665 

610.625. 

971.518 

Xy/.n 

i^/secf_ 
it /’see 

.0117+3 

.02237 

-. '01281 

-•.03179 

-.04306 

-.05243 

-.06120 

-.06870 

X„/m 

ft/'se^' 
it /sec 

.04360 

.05111 

.10637 

.10919 

.10653 

.12448 

.10628 

.07771 

Xg/.n 

'it7s'ec2"““ 
raclT’sccT • 

1.02536 

.93903 

1.55534 

1.45965 

1.37043 

1.88836 

l;41565 

.89147 

X^^/m 

Tt'/sec'-‘‘ 

Tn 

.16625 

. 12347 

.08223 

.12593 

.14599 

-.06972 

-.05583 

■ -.03761 

X6j./m 

it/scc^ 

1.13928 

.94352 

.82126 

.880,28 

.82702 

1.02567 

.86407 

.62703 


it/sec^ 

Tt7se_c“ 

- . 13435 

- . 23640 

-.12619 

-.02151 

.01777 

. .06458 

.06927 

i 05361 


Tir/’sfTc'^ 

itr 

- . 28748 

-.26321 

-.82254, 

-.90575 

-.98052 

-1.06060 

-1.10549 

-1.16513 

2,|/» 

ft7sec^ 

radTsec 

-.79039 

.30060 

-2.45731 

-2.40596 

-2.42975 

-2.83028 

-2.86944 

-3. OS 9 24 

hjm 

it/sec^ 

xiT 

.00615 

.19996 

,94841 

.77315 

.71536 

.65350 

.60489 

.5-W52 


ft/sec® 
"7. n 

-7.45257 

-7.38250 

-6.35173 

-8.00904 

-9.05358 


-10.65093 

-11.27463 


racl/sec^ 

JJ:7's1-clT. 

.01308 

.01863 

-.00906 

-.00696 

-.00635 

-.00238 


' - .^0142 ’ 

W,/>yy 

racl/scc“ 

.00368 

.03420 

.00228 

"’65^69 

.00683 

.004W 

!6o303 

.00220 

v: 



rncl/ 

)?aci/.s7ic * 

^'758241 

- .81810 

-1.48310 

-1.56553' 

-1.62309 

-1.72093 

-1.69157 



-1.61660 
. 56768 


raci/soc^ 
ij) 

.34363 

.30612 

.43861 


”'”.'50338 

.53/75 

. 55538 

>'6,/i.yy 

rn cf/ soe^' 
"ITf 

-’.'05591 

- . 16472 

- 71536^ 

!'iTg2o 

7i.0884 

'.oo’ool' 

"'6s84.3’ 

.OS/: 00 


D.'ifa c>:l:):n|)olal:(:‘cl from valuoa a(:.80 knof.a. 
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Table D-5. YHC-IA Lateral Stability Derivatives (Rate of Descent = 1500 ft/min). 


I'OllWARD 

VELOCITY 

KTS 

0 

20 

40 

60 

80 

100* 


"a 

o 

in 



.01404 

.05177 

.09882 

.21705 



in 

-.07915 

.32928 

.06918 

-.19696 

-.43185 

r -.57289 

-. 86605 !-. 96015 

f 

wsm 

ft/ sec® 
fU/sec 

-.08635 

-.04068 

-.11126 

^-.15170 

-.18262 

-.20181 



Yp/m 


-1.10440 

-1.30270 

-1.78540 


-1.73270 

-1.52224 

-1.26548 

-.79301 

Yp/m 

ft/sec® 
rad/ sec 

-.18435 

-.10798 

- .30137 

-.22815 

-.19244 

-.00045 


mm 


ft/sec® 

in 

.98054 

.97047 

.92561 

.89687 

-.87675 

.87843 


1.00945 


ft/sec^ 

~Tn 

.14250 

.11217 

.09751 

.07301 

.07008 

.08800 

.11042 

.12421 

VLx 


-.00996 

-.00799 

-.00448 

-.00371 

-.00289 

-.01035 

-.02571 

-.02927 

S'Lx 

rad/ sec® 
rad/ sec 

-.63036 

-.71659 

-.86876 

-.87455 

- . 79734 

-.69428 

-.58319 

-.40661 

VLx 

rad/ sec ^ 
rad/ sec ' 

-.05234 

-.02797 

-.10977 

-.07677 

-.05688 

.02671 

.04799 

.13095, 


rad/ sec® 
in 

.45903 

.45705 

.44170 

.43230 

.42630 

.42645 

.44179 

.47353 

^6 /I 

r XX 

rad/ sec* 
in 

-.12454 

-.13516 

-.13185 

-.13554 

-.13256 

-.12602 

-.12594 

-.13779 

N„/I 
v' zz 

rad/ sec® 
ft/sec 

.00122 

-.00016 

-.00236 

-.00456 

-.00706 

-.00827 

-.00668 

-.00351 

"p'Lx 

rad/ sec® 
rad/, sec 

-.02085 

1 

-.01101 

-.02424 

-.04001 

-,05921 

-.07441 

-.08324 

-.08048 


rad/ sec® 
Yad/slic 

-.04986 

-.02575 

-.03287 

-.03310 

-.03704 

-.06380 

-.06944 

-.10919 


rad/ sec® 
in 

.02925 

.02744 

.02580 

,02385 

.02313 

.02376 

.02537 

.02650 

"s/Lx 1 

rad/ sec® 
Lg_ — ■ 

.17285 1 

. 17149 

.16373 , 

i 

.15906 

,15540 

.15545 1 

,16332 

.17935 


*Dat:a extrapolated from values at 80 knots. 
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Table D-6, YHC-IA Longitudinal Stability Derivatives (Rate of Descent = -1 500 ft/min) 


0 

•< 

0 

m 

S 

0 

• 

1 

n 

• 

• 

0 

z 

in 

< 

Z 

in 

ID 

8 7 

l t. 

X j 

• 

00 

c 

a 

r 

z 

0 

H 

0 

z 

> 

(D 

0) 

> 

n 

I 

c 

09 

3 

OD 

□ 

O 

0 

u 

• 

I'omoAui) 

\'i:i,OCT.TV 

KTS 

0 

20 

40 

60 

80 

lOO* 

).20" 

140’'- 


‘Ic-l’. . 

9.09f)i:i 

8.24561 

6.51804 

4.34680 

1.74707 

.85800 

-2.0445C 

■ -7.20897 

O 

In 

.566U 

.45244 

. 38166 

.63183 

1.02115 

.18179 

.40690 

,58596 

6 '■ 

in 

6.51301 

6.30131 

5. 94 108 

5.94845 

6.38886 

7,21746 

8.58952 

10.55994 

f 

Ib/hr 

1102.055 

1076.035 

1023.741 

1012.939 

1055.591 

1128.295 

1170.255 

1531.148 , 

Xu/m 

ft/ SCC2 
I'k/'sec" 

-.02649 

-.01935 

-.03260 

-.04491 

-.05277 

-.06214 

-.07091 

-.07641 


t't/s'ecs 

I'C7soc 

.07505 

' .08122 

.08209 

.07521 

.06568 

,08363 

.06543 

.03686 

x,/» 

racf/^'ec 

.15633 

.20875 

.26690 

.. . 

.17058 

.09717 

,61510 

.14239 

-.38179 

X5 /m 
e 

L't/scc^ 

in 

.17842 

.16228 

. 14484 

.16227 

.18105 

-,03466 

-.02077 

-.00255 

^6^, /m 

ft/ sees 
"**'“11^ 

1.16992 

1.04000 

.90570 

.78861 

.62563 

,82428 

.66268 

.42569 


TF/ see” 
f't/sec 

To^TT^ 

-.02454 

-.00971 

.00737 

.02520 

.07201 

||||H 

.06104 


Et/ sec'^ 
‘f t/Fec. 


-.60121 

-.72925 

-.85332 

Ha 

I -1,02179 

-1.06668 

-1.12632 

BBHI 

ft/secs 
rad/ see 

-1.07268 

-1.20010 

-1.25010 

-1.38662 

-1.47426 

-1,87479 

-1.91395 

-2,13375 

^6e/ni 

‘Tt/secs 

Tn 

.05444 

.16010 

.36218 

.41296 

.41472 

,35286 

.30425 

.25688 


ft/ sec® 
in 

-7.5865 

-7.6A041 

-8.2477 

~~-8.95^ 

-9.7090 

-10.48495 

-11.30635 

-11.93003 


raci/ sec® 
ft/ see 

.00251 

||H 

.00321 

-.00458 

-.00467 

-.00070 

-.00005 

.00025 


1 rad/ see® 
j t 


.00327 

.00726 

.00765 

.00730 

.00532 

.00350 

-1.49278 

. 00267 
-1.41/S.l 


! rncl/sec® 
Jl'AUrZsrnc__ 


-1.03020 

-1.24040 

-1.37428 

-1.4243 

-1.52214 

"‘o"yy 

rad/ SCO® 
in 

.37472 

.37930 

.40093 

.43815 

.46558 

.49995 

.5175S 

.529SS 

■•'6,/X.yy 

rad/ see® 

-.04570 

-.03118 

.01909 

.03637 

.04212 

.02194 

.01435 

.00992 


b'ltii cxl:i;nj)olriI:ecl from values at 80 knots , 
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Table D-7. YHC-IA Lateral Stability Derivatives (Rate of Descent = -1500 ft/min). 


FORWARD 

VELOCITY 

KTS 

0 

20 

40 

60 

80 

100* 

120* 

140* 

o 

in 

.16567 

.15081 

.14181 

.15535 

.20065 

.31888 

.46288 

.51919 


in 

-.23383 

-.21772 

-.28008 

-.42584 

-.61869 

-.75973 

-1,05289 

-1.14699 


ft/ sec® 
ttT^ec 

-.02060 

-.08487^ 

-.10455 

-.12256 

-.14842 

-.16761 

. 18487 

-.22163 

Yp/m 

ft/sec“ 
rad/ sec 

-.40249 

-.48197 

-.63817 

-.68143 

-.63180 

-.42134 

-.16458 

.30789 

Yp/iti 

ft/ sec® 
r^/sec 

-.04204 

.04052 

-.06921 

-.04086 

.01514 

.20713 

.19847 

.29458 


ft/ sec® 
in 

1.03950 

1.04228 

1.04066 

1.04266 

1.06287 

1.06455 

. 1 . 10945 

1.19557 


ft/ sec^ 
in 

.14982 

.14802 

.12894 

.11875 

.12938 

.14730 

.16972 

.18351 


rad/ sec® 
ttZsec • 

-.02127 

.01846 

-.01978 

.02294 

-.02866 

-.03612 

-.05148 

-.05504 

>-p"xx 

rad/ sec® 
ra^T/ sec 

-.37496 

-.40345 

-.45470 

-.46318 

-.42905 

-.32599 

-.21490 

-.03832 

L,/I 
r' XX 

rad/ sec. 
rad/ sec 

.02907 

.06798 1 

.02113 

.03740 

.07088 

.15447 

.17575 

.25871 


rad/ sec® 
in 

.48113 

.48228 1 

.48271 

.48417 

.49124 

.49139 

.50673 

.53847 

‘•‘r'^xx 

rad/ sec® 

-.13312 

- . 13406 

-.14131 

- . 14556 

-.14455 

-.13801 

-.13793 

-.14978 

V'xx 

rad/ sec® 
ft/ sec 

.00108 

.00098 

.00044 

.00052 

.00041 

-.00080 

.00079 

.00396 

N 7l 
p zz 

rad/ sec® 

-.02133 

-.02122 

-.02478 

-.02754 

-.03643 

-.05163 

-.06046 

-.05770 

N^/I 

r zz 

ra^ sec® 
fau/ sec 

-.07653 

-.06348 

-.07218 

-.07450 

-.08262 

-.10902 

-.11466 

-.15441 

Na /I 

rad/ sec® 

.03076 

.03053 

.02424 

.02856 

.02918 

.02981 

.03142 

.03255 

"6^/'xz 

rad/ sec® 
in 

. 18338 

.18355 

.18366 

.18425 

.18730 

.18735 

.19522 

.21125 


*Data extrapolated from values at 80 knots. 



APPENDIX E 

STABILITY DERIVATIVES CH-47 HELICOPTO? 


This appendix presents the stability derivatives provided by NASA LaRC for 
the CH-47 VALT research helicopter, which will be used to evaluate the advanced 
display concepts for commercial VTOL aircraft. The column labeled 0.00000 was 
derived by fitting a third order curvp through '■40, -20, +20 and +40 values. 
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